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Prefacio 


¿Cómo es que llegamos a que en el XXXII Colóquio Brasileiro de Matemática haya 
un curso (y un texto) como éste?. 

Interesante pregunta cuya respuesta pasamos a detallar. Esta respuesta puede, cómo 
no, ser tomada como un ejemplo del llamado efecto mariposa en lo social. 

Por allá por el año 2015 dos estudiantes de pedagogía en matemáticas que estaban 
iniciando sus estudios de magister en matemáticas (Pablo Díaz y Juan Vivanco) se pusie- 
ron de acuerdo con un profesor (Rafael Labarca) para estudiar, semanalmente, algunas 
materias básicas para la formación de un matemático. 

Las asignaturas básicas que se debería cursar en dicha formación incluyen: álgebra 
lineal; álgebra exterior; teoría de grupos; estructuras algebraicas; cálculo de una y varias 
variables; variable compleja; análisis real; espacios métricos; topología; análisis abstracto; 
teoría de integración; ecuaciones diferenciales ordinarias y parciales; elementos de proba- 
bilidad e inferencia estadística; Cálculo numérico; Fundamento de la matemática; Geome- 
tría diferencial; análisis geométrico; otras geometrías; cálculo de variaciones; elementos 
de sistemas dinámicos, entre otras. 

Habiendo tanta posibilidad debieron elegir. Se pusieron de acuerdo e iniciaron con 
un repaso de cálculo diferencial e integral de una variable real; siguieron con métodos 
de cálculo para resolver ecuaciones diferenciales ordinarias y parciales. En este estudio 
aparecieron, naturalmente, las funciones Gamma y Beta. En este punto los tres percibieron 
que eran muy pocas las demostraciones asociadas a las propiedades de la función Gamma y 
surgió, naturalmente, la idea de realizar las demostraciones de varias de estas propiedades. 

Luego, siguieron con elementos de sistemas dinámicos y, por ahí, nació la idea de 
estudiar la dinámica de la función Gamma. Para esa fecha Pablo Díaz estaba completando 
los cursos básicos del programa de magister en ciencia en la especialidad de matemática 
en la Universidad de Santiago de Chile y, cómo no, la propuesta del profesor fue empezar 
a entender la dinámica de la función Gamma. A fin de saber si este tema había ya sido 
estudiado contactó a colegas del área de sistemas dinámicos en Brasil, Francia, Estados 


Unidos y México; para saber si esto ya había sido estudiado y, sorpresa, no había sido 
estudiado y, por ello, empezaron dicho estudio. 

Al Iniciar el trabajo para la tesis (alumno Pablo Díaz; guía Rafael Labarca) quedó 
claro que lo que se quería era escribir un texto que tuviera dos partes: en una de ellas se 
desarrollarían diversas propiedades clásicas de la función Gamma y en la segunda parte 
se estudiaría la dinámica de dicha función. 

Es claro que todo el material que se quería desarrollar excede lo que se pretende con 
una tesis de magister pero, igual se hizo. Parte de ese material constituyó la tesis de magis- 
ter de Pablo y la base de un trabajo original, firmado por Díaz y Labarca, que se pretende 
publicar en una revista matemática. 

Este material se completó el segundo semestre de 2018 y, para esa fecha, la orga- 
nización del Coloquio Brasilero de Matemática hizo un llamado público para presentar 
propuestas de cursos introductorios y avanzados para el XXXII Coloquio. Se hizo la pro- 
puesta de curso introductorio y, nueva sorpresa, fue aceptada. 

Entendemos que el comité científico aceptó la propuesta porque revisita un tema anti- 
guo (la función Gama) utilizando una visión moderna (desde la perspectiva de los sistemas 
dinámicos). 

Eso; humm, y,..... ¿dónde entra aquí esto del efecto mariposa?. El efecto mariposa se 
explica aquí cuando vemos que un genuino interés por saber matemática de tres personas 
en algún lugar lleva, a lo largo de un par de años, a concluir un libro que, esperamos, resulte 
entretenido para muchos otros aprendices de matemáticos y/o estudiantes de ciencia y 
tecnología en toda la América Latina. 

También, estamos seguros que en el texto hay más de algún error y, con seguridad, 
algunas demostraciones podrían haber sido más cortas o más elegantes. Por eso, y para 
quienes estudien el texto, les dejamos nuestros correos electrónicos (pablo.diaz(Wusach.cl 
y rafael.labarca(Wusach.cl) para que tengan la gentileza de permitirnos progresar en una 
mejor versión del mismo ya que, como no, podemos soñar con una segunda versión. 

Para finalizar, y como una manera de recordar a uno de los pioneros en la investigación 
en matemática en el Brasil y fundador del Coloquio Brasilero de Matemática, dedicamos 
este libro al Profesor Dr. Mauricio Matos Peixoto, recientemente fallecido en Río de Ja- 
neiro. 


Pablo Díaz y Rafael Labarca. 
San Carlos- Brasil y Santiago de Chile, Mayo de 2019. 


Introducción 


Todo indica, ver las referencias Davis (1959) y Fuss (1843), que la función Gamma 
tendría su origen en un intercambio epistolar entre dos destacados matemáticos del siglo 
XVIII, los que se correspondían con la idea de resolver el problema de extender la función 
factorial, definida para números naturales, a números no naturales. 

Introducimos, de esta forma, el estudio de la función Gamma que se desarrolla, a lo 
largo de tiempo, con contribuciones de una amplia gama de matemáticos que han trabaja- 
do, diversas veces, de manera colaborativa. En los inicios, muchas de sus propiedades se 
obtienen usando una variada gama de técnicas e ideas matemáticas como series, límites, 
derivadas e integrales; que eran la frontera del conocimiento en la época. 

Si el lector piensa, por el hecho de que el estudio de la función Gamma tiene ya casi 
300 años, que no hay nada nuevo que ver, está equivocado pues la dinámica de la función 
Gamma no ha sido objeto de estudio. 

Este trabajo contiene, básicamente, dos capítulos que se refieren a la función Gamma: 


FT:R-—(0,-1-2,...) >R-—(0) definida por F(x) =/ ÓSEA 
0 


En el primer capítulo establecemos diversas propiedades y resultados de esta función 
como: sus diferentes representaciones, la aproximación de Stirling-De Moivre, la fórmula 
de reflexión y algunas funciones que se derivan de la función Gamma (en especial la 
función digamma). 

En el segundo capítulo, nos centramos en el estudio de la dinámica de la función Gam- 
ma. Para ello, lo primero será verificar que para x > O hay dos puntos fijos, a saber x = 1 
y x = 3,5623828215736; y un punto de mínimo; a saber xy = 1,4613361518823 con 
TP (xo) = 0.8855226070359. Luego, probaremos que el punto fijo x = 1 es de tipo atrac- 
tor y que el punto fijo x = 3, 5623828215736 es de tipo repulsor. La dinámica para x > 0 


resulta ser bastante simple: o la órbita es atraída para el punto atractor o es atraída por el 
punto al infinito. 

La riqueza de la dinámica de la función Gamma se encuentra en x < 0: En particular 
nos interesa el conjunto de valores negativos cuya imagen I'(x) < 0. Estos valores se 
encuentran en los intervalos de la forma: ]-2k — 1,—2k[,k € N. En el interior de es- 
tos intervalos consideraremos los subintervalos Jz = [-x4, =xf] e lx = [—yé, =yFl, 
donde los extremos satisfacen: 


res Ms A=2<=1, Fez = MS A=22 


y para los que podemos probar que | '(z)| > 4 > l;parak € N yz € Jg¿U Iz. Con 
esto concluiremos que el conjunto maximal invariante An = [Zo (Uj=1 Li), con L; = 
FJ; U [; para cada i = 1,2,...,n;es hiperbólico y tiene una dinámica topológicamente 
equivalente a la dinámica del shift de 2n símbolos. De hecho, la biyección que hace la 
conjugación topológica se define desde A, hacia el conjunto de sucesiones: Xy = (0 : 


No > (1,2,3,...,2nj. Esto es, asociado con cada z € A, tenemos la sucesión: 
0(=) € Lan definida por O(21(1)=j => T'(2)eL; 


esto define la aplicación z > 0(z) que resulta ser una biyección de A, en 2»,. Esta 
biyección nos permite definir una topología en 23, por :A C 2», es abierto si y sólo si 
971(A) C A, es un abierto. 

En 2, consideramos el shift o : 2, > 22 definido por 0(0)(1) = 0;+,1, esto 
es 0(00,01,...) = (01,02,...), cuya entropía topológica es igual a log(2n), es decir, 
hroplo) = log(2n) Labarca (2011). Probaremos que O : A, > 2, es una conjugación 
topológica entre P'| y, : An > An y 0 : Lan > 2, (estoesoo0= 00 T|z,), con 
lo que concluíremos que A;op (T'| 1,) = log(2n). Cómo consecuencia tendremos que la 
entropía topológica de 1” es infinita. 


Función 
Gamma 


1.1 Antecedentes Históricos 


De acuerdo a lo que señala P. H. Fuss, en Fuss (1843), el surgimiento de la función 
Gamma, se remontería a un intercambio de cartas entre dos grandes matemáticos del siglo 
XVIII. El primero de ellos fue nada menos que Leonhard Euler (1707-1783), que a la 
sazón tenía 22 años de edad y ya se había convertido en un gran matemático. El segundo, 
menos conocido y mayor en edad, era Christian Goldbach (1690-1764) que a la fecha era 
secretario de la Academia de Ciencias de San Petersburgo. 

La función Gamma habría surgido de la intersección de dos desarrollos matemáticos. 
El primero era el problema de interpolación, que se había hecho muy popular durante el 
siglo XVIII. El segundo era el cálculo integral. 

Goldbach considero el problema de la interpolación de los números factoriales 1! = 
1,2! = 2,3! = 6,4! = 24,........ Es decir, definir (o encontrar) una función que diera 


A . . S a Ñ 
sentido, por ejemplo, al factorial de 7 o sea definir la cantidad (5) !, como un valor 


numérico. 
En una carta de Daniel Bernoulli (1700-1782) a Goldbach, fechada el 06 de Octubre de 
1729 y escrita en San Petersburgo, Fuss (ibíd.) tomo II, le propuso la siguiente expresión: 


xy 1 2 3 4 A 
(4+ ) . ñ aa 
2 l+x 2+x 3+x A-=1+x 
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como una interpolación del factorial de un número, donde la precisión del resultado au- 
menta mientras A sea más grande. Por ejemplo para x = 4 y A = 6 tenemos que la 
expresión nos da el vaor 24, 380952 para el factorial de 4. (En la carta Bernoulli calcula 
el valor de la expresión con A=8yx= 2). 


Por esa fecha Euler desarrollaba independientemente su propia representación del fac- 
torial de un número, y animado por la expresión de Bernoulli, le escribió una carta a 
Goldbach (Davis (1959) tomo I) en la cual representaba los n—términos de la sucesión 
1,2,6,24,... por el producto: 


AO AO 


Euler afirmó que su producto convergía en el infinito, es decir, mientras más términos 
había, entonces el valor de esta expresión se aproximaba a n!. Aún cuando los productos 
de Euler y Bernoulli son diferentes, ambos convergen al mismo número en el infinito, y 
se verificó que la fórmula dada por Bernoulli converge más rápido. 


Entretanto, y en la carta enviada por Euler a Goldbach, no entregaba mayores detalles 
sobre la interpolación del factorial de un número. Sin embargo, y un año más tarde en 
1730, Euler entregó más detalles en su articulo Euler (1738), titulado “De progressioni- 
bus transcendentibus seu quarum termini generales algebraice dari nequeunt”. En este 
articulo, Euler afirmó que este producto era “maravillosamente adecuado para interpolar 
el factorial de un término, cuyos índices son fraccionarios”, pero él tenia la intención de 
entregar métodos más convincentes que le proporcionaran valores exactos. 


Este método alternativo daría como resultado una segunda interpolación, entregada en 
forma de una integral definida. El impulso para esta nueva dirección parece haber sido 


Ñ ; . 1 A 
dada por el factorial de 3 Euler, en su articulo, sustituyo n = 2 en su producto, es decir, 


considero la expresión 


DADA: 


que podemos escribir de la forma: 


1.9” gl=n ,3n 3gl=n, qn qlon,5n 
l+n  2+n  3+n  44+n 


> 


1 
entonces al reemplazar n = 3 obtuvo: 


242 242-443 243.44 244: 45 
ss "CC TT 7 
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y ordenando cada término dentro de una raíz obtenemos la siguiente forma: 


y SE 4-6 6-8 8-10 
3.3 5:5 7:17 9:9 
y esta última relación ya había sido estudiada de forma indirecta por John Wallis (1616- 
1703). En su libro ”Arithmetica Infinitorum” de 1656 Davis (1959), Wallis había estable- 


cido que el área de una circunferencia con diámetro unitario podía ser expresado por el 


producto infinito 
2.4 4-6 6:8 8-10 


33-35 71 DO 


1 
Por ello Euler concluyó (En el punto 2 de su articulo) que 5% era igual a la raíz de esta 


sE 
-. 


1 
área, estableciendo el resultado central: 7) = 


Ahora bien, 7 aparece en el cálculo del área del círculo, y esta área significaba a su 
vez integrales, pero Euler ya estaba familiarizado con la integrales que exhibían el mismo 
fenómeno. Por lo tanto se le ocurrió buscar una transformación que le permitiera expresar 
su producto infinito en forma de integral. Para hacer esto usó una integral de la forma: 

1 


/ x(1—x)dx. 


Mencionamos que casos particulares de ésta integral habían sido considerados por 
Newton, Stirling y Wallis. Tomando n como un número entero y e un valor arbitrario (no 
es la constante e que conocemos hoy en día), Euler expandió (1—x)” en su forma binomial, 
y sin dificultad encontró que: 


1-21 
(e+ D(e+2)---(e+n+1) 


1 
/ x*(1—x)"dx = 
0 
La idea de Euler era aislar 1 -2---n del denominador para obtener una expresión para 


n! en forma de integral. 


Ahora seguiremos con la misma nomenclatura e ideas de Euler, marcando con un (*) 


las fórmulas que entrega en el articulo original. Él escribió la integral ' como / y 
0 


sustituyó e por f y obtuvo: 
g 


1 £ gra 1.-2--n 
E(L—x"dx = : 
[> CSS o OPA 
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entonces: 


1.2---n F+(n+1e £ E 

= 2] xs dx(1-x)”. 
(F+ES +28) +18) yo 

Él observó que podía aislar 1 - 2--=1 si hacia f = 1 y g = 0 en la parte izquierda 


de esta igualdad, pero en ese caso, obtuvo en la parte derecha de la igualdad una forma 
indeterminada, la cual escribe extrañamente como: 


(a*) 


/ x0dx(l—x)" (0%) 


p9r+1 


Luego procedió a encontrar un valor para (b*). Primero hizo la sustitución x = z/+s en 
lugar de x y esto le dio : 


£ 
dx = z 3+f dz (c*) 
FE 
reemplazando estos valores en la parte derecha de la igualdad en (a*) se convierte en: 
Ff+(M+De 
gnel cr (127%) dz. (d*) 


esta expresión podemos reescribirla de la forma: 


A ira Y 
FF Jo A e/F +28) 


Una vez más, Euler probó con los valores f = 1 y g = 0 y obtuvo la forma indeterminada: 


_»0yn 
JE 7 dz. (e*) 


(12 


Ahora él considero la relación G(t) = , y para valores cercanos a £ = 0, derivo 


el denominador y el numerador, aplicando la regla de 1”Hópital y obtuvo: 


=2*dtlz 
A f* 
di (1) 
con lz = In(z), entonces para 1 = O se obtiene —/z, por lo que: 
1 — 0 
a 1 (2) 


0 
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CA ay, (1) 


Por lo tanto, Euler concluyó que: 


1 
nt=/ ElnzYdz. 
0 


Con esto Euler obtuvo una expresión para n!, como una integral que se puede evaluar 
para valores distintos a números enteros. Observamos que la integral de la derecha es una 
integral recursiva que se calcula usando integración por partes. 

Notamos que, en ese tiempo, ésta integral de Euler no era vista como una función en si 
misma (de la variable x), sino como una fórmula para obtener los valores de los factoriales 
de números no-naturales. 

Fue Adrien-Marie Legendre (1752-1833) quien estableció el valor de la integral cómo 
la definición de una función. En efecto, en su “Exercices de Calcul Intégral”, Vol. 1 (1811) 
Legendre (2010), ( páginas 339 a 343) Legendre realizó diversos cambios de variables 


para llegar a la integral: 
1 1 x—1 
ru)=/ (m2) dt, x>0 
0 


que definió cómo (x — 1)!; así y con la notación actual, la función Gamma había nacido. 
Haciendo el cambio de variable u = In (2, obtenemos la representación tradicional y 
más conocida de la función Gamma, también conocida como la integral de Euler: 


[,0) 
T(x) =/ ute dy x>0. 
0 


Legendre, utilizando la definición de la integral para la función Gamma probó que 
T"(1) = 1 y la relación de recurrencia "(x + 1) = xI'(x) = x! para x > 0. 

De esta forma la función Gamma extiende el factorial de números naturales a todos 
los números reales positivos. Sin embargo, y durante casi 200 años de investigación de 
esta función surgieron las siguientes preguntas: ¿la función Gamma es única con estas 
propiedades?, ¿es la única forma de interpolar el factorial de un número no-entero?, las 
respuestas a estas preguntas llegaron bajo el concepto de convexidad. 


Observacáo. Una función real f(x) se dice convexa en el intervalo (a, b), si para todo 
c € (a,b) la recta tangente en c, al gráfico de f, esta por debajo del gráfico. 


En, 1922, Harald Bohr (1887-1951) y Johannes Mollerup (1872-1937) demostraron 
que de todas las funciones, log-convexas, que extienden el factorial de un número real 
positivo, la función Gamma es la única que cumple con (1) = 1, (x +1) = xI'(x). 


6 1. Función Gamma 


Observacáo. Una función f(x) definida y positiva en un cierto intervalo se dice log- 
convexa si log( f(x)) es una función convexa. 


1.2 Propiedades de la Función Gamma 


Definición 1.2.1. Para x > O se define la función Gamma como: 
[0,0] 
To) = / in 
0 


De está definición de la función Gamma, podemos verificar que está bien definida para 
x>0. 


Proposición 1.2.1. I'(x) esta bien definida para todo x > 0. 
Prueba. Esto quiere decir que I'(x) € R. 


00 1 o) 
T(x) za pd El pa ear] Poetas, 
0 0 1 
A AAA 


1 L 


fi)  — Six => 1, entonces f¡ es convergente puesto que su integrando es una función 
continua (y entonces integrable) en [0, 1] (integral de Riemamn). 


— Si0 < x < l entonces, y para todo £ > 0, tenemos la siguiente desigualdad: 
0 < 171871 < 471. Por el criterio de comparación, si ocurre que la integral 


1 
/ podi converge, entonces /¡ converge. Ahora: 
0 


+ 
1 1 [e 1 1 
/ ¿=1dt= lim Ptdt= lim =|.='3, 
o+ a>0+ Ja a>0+ X |y Xx 
por lo tanto, /; converge para cada x > 0. 
LS A pe y F() 
I,) Usando el criterio de comparación al límite con g(t) = 17“ tenemos 20 = 
g 
(Lt 
= pe + let entonces 
(2 
. t (+1 
lim F6) = lim (**e*= lim = 
t>+00 g(t) t>+00 t>o+oo el 


1 
y como g(t) = Po converge, la integral de f (1) también converge, por lo tanto, 12 
converge. 


1.2. Propiedades de la Función Gamma ñ 


Proposición 1.2.2. La función Gamma satisface: (1) = 1. 


Prueba. Fácil, solo hacer x = 1 en la definición 1.2.1. 


Proposición 1.2.3. Six > 1, entonces I'(x) = (x — 1) '(x — 1). 
[0,2] 
Prueba. Dado I'(x) = / 147227 dt, integrando por partes, con u = 1%! y du = 
0 
e”* dt, tenemos: 


a 


lim Pelada 
a-—00 0 


Il 


a [0,0] 
+/ =D SRrd 
0 0 


lim (=1471e7*) 
a—>00 


Il 


(«—1) IN PRrtdt= (=D (r= 1D. 


por lo tanto, P(x) = (x— DD) (x—-1), Vx > 1. 
Proposición 1.2.4. Six € N, entonces T(x) = (x — 1)! 


Prueba. Usando la proposición anterior tenemos que I'(x) = (x— 1)" (x—1); aplicando 
nuevamente la proposición pero a (x—1)six—1 > 1, tenemos I'(x—1) = (x-2)I (x— 
2). Sucesivamente, IP (x) = (x — 1)(x — 2)(x — 3)---2-1-1'(1) = (x— 1)! 


Proposición 1.2.5. Six > 1 y x = [x] + r dónde [x] indica la parte entera de x y 
r € [0, 1], entonces: T(x) = (x— D(x—2)---(x— [xD (1). 


Prueba. Sabemos que TP (x) = (x — 1)P'(x — 1). Además, como x > [x] con 0 < 
x — [x] < 1, entonces 
A 


r 


TG) = (2 Dr 2) ++ (x — [xP G— [x) = (4 DG — 2) ++: (x — [xD (1) 


como señalado. 


Observacáo. El resultado anterior nos dice que basta con conocer los valores de la fun- 
ción Gamma en ]0, 1[ y multiplicarla por la expresión polinomial (x— 1)(x—2) ---(x—n) 
para conocer sus valores para x tal que x e]n,n + 11. 


[0,2] 
Proposición 1.2.6. Se cumple que: / e” dx = VI. (Integral de Gauss) 


—00 
[oe] 


2 
Prueba. Sea J =/ e * dx, entonces 


0 
2 2 cd 2 
ps] tax] e dy 
0 0 
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aplicando el Teorema de Fubini a la integral doble 


[o,0) [,0) 2 2 
/ / e 0RID ddy 
o Jo 
[,0) [,2) 2 2 [,0) 2 [,0) 2 
l 0 ey dedy =/ e* dx- | e Y dy =3?. 
o Jo 0 0 


Ahora, usando el cambio de coordenadas x = rcos(0) e y = rsin(0), para calcular la 
integral doble, tenemos: 


tenemos 


2 Eje 
J =f pos rdrd0 = di re Par d0 = 
¿2 1 _,2 Tr. 
= 3 im Eze ¡l , = 5 dim, [37 + | Sa 
luego, J? = 7) entonces: J = sE . Por tanto: [%, e dx= yr. 


Observacáo. Mediante el cambio de variable + = y? en la definición 1.2.1 se tiene que: 


2 PA 
Tx) = 2/ y Tady 
0 


1 
Proposición 1.2.7. (5) = yr. 


Prueba. De la observación obtenemos que 
E da o 2 ym 
de (3) = 2/ e”? dy y usando la proposición 1.2.6 tenemos que | e * dx= 2 
a 0 
1 
y, por tanto: f” (5) =D: E == 


Como consecuencia tenemos: 


3-.-(2n—1) 
9n 


1 
Proposición 1.2.8. T' (n + 3) = V/Tr, para todo n € Ñ. 


Prueba. Usando la proposición 1.2.3 tenemos que TP (x + 1) = xP (x), entonces podemos 
reescribir T” (n + 3) de la forma 


P(n 14341) =(n 1+3)0 (1 1 + 3), pero a su vez 
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dE (r 1+ 3) =P ( 2+ + 1) y así usamos esto n—veces tenemos que: 


1 1 1 
r(a+z)=(2-1+3)1 (n-1+3) 
=(n-1+3)-(2-2+3)r (n-2+3) 
2 2 2 
E 1 1 1 r 
=(n-1+3)-(2-2+3)--(5) 
2n- 1 2n-3 3 1 
= RE ME 


por lo tanto: 


1 1.3...On—1 
r(a+3)- El IE 


Generalizando esta última proposición se obtiene la siguiente: 

2) _ [0 Dk + pl: [(n =2)Kk 7 
k kn 
Prueba. Para demostrar esto, usaremos la proposición 1.2.3, es decir, podemos reescribir 
rP (rn + 7) dela forma (rn =-1+ z + 1), entonces T (r =1+ z + 1) = (rn -1+ -) r ( 
pero asu vez ” (rn =1+ 2) =I (r -2+ A + 1) y así usamos esto n—veces tenemos 
que: 


Proposición 1.2.9. Para todo n, p,k € N, secumple T ( + 


np pre 


cl a ra 


por lo tanto: 


r(n+2)= A E E 


para todo n, p,k € ÑN. 


Veamos ahora cómo extender la función Gamma a los números negativos. Para ello, 
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T 1 
aplicamos la proposición 1.2.3 y obtenemos la relación: (x) = ERE 
x 
Ejemplo 1.2.1. Para x = =) tenemos: 
N T( 1 
Pp (-3) = (2) =-2T (5) =-2Vx 
2 => 2 


Iterando este proceso para x < 0, x 4 —1,—2,-—3,... tenemos la siguiente: 
Definición 1.2.2. Para cada x tal que —n <x <—n+1,n e N se define 


T(x+mnm) 
x(x+1)---(x+n-=1) 


PQ) = 


De la proposición 1.2.3 tenemos la siguiente: 


Proposición 1.2.10. lim P(x) = +00. 
x>0 


TMx+1 
Prueba. Tenemos que I'(x) = ED entonces: 
x 
lim T(x+1) 
Tr 1 tE 
ln Pip PRL = +00. 
x>0+ x>0+ x lim x 
x>0+ 
Obtenemos ahora el siguiente resultado: 
Proposición 1.2.11. lim T(x) = —00 
x>0 
ad T(x+1) 
Prueba. De la definición 1.2.2, con n = 1, tenemos que I'(x) = ————, entonces: 
XxX 
| o PEAD A tD 
lim (x)= lim = - = —00. 
x>0— x>0— x lim x 
x>0 


Otra de las consecuencias de la definición 1.2.2 es que I'(n) > too para valores 
enteros negativos y que el signo de /"(x) en el intervalo —n < x <—n+1,n € N será 
(-1)” y por tanto la gráfica de la función Gamma para x € (R/(0,—1,—2,-—3,...)) será 
de la forma: 
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| Tr) 


Figura 1.1: Gráfico de la Función Gamma 


1.2.1 Funciones Log-Convexas 


En esta sección probaremos la unicidad de la función Gamma en el conjunto de funcio- 
nes log-convexas. Esto es, veremos que si una función satisface las propiedades señaladas 
en las proposiciones 1.2.2 y 1.2.3 y además es log-convexa, entonces coincide con la fun- 
ción Gamma. Para ello, definiremos este concepto y presentaremos algúnos resultados 
relacionados (ver Andrews, Askey y Roy (1999), Artin (1964) y Carlson (1977)). 


Definición 1.2.3. Una función f : (a,b) > R es convexa si para cualquier x, y € (a,b) 
se cumple f(Ax + (1—A)y) <AF(6) + (0 —A)FO), VA € (0, 1). 


Teorema 1.2.12. Una función f : (a,b) > R es convexa si y solo si para cualquier 
a<x< y <z<b, se tiene: 


0)= 00 MH) TO) 
yor o z=y 


Prueba. Sean los valores x, y y z tales que: a < x < y < z < b. En este caso tenemos 

=- =- z= =x ' 

que EE ad ette 22% Y Ahora como y = dee z, aplicando 
Z=Xx Z=Xx Z=x Z=x 


f tenemos: 


A 


entonces (2 —x) f(y) < (2— y) f() + (y —x) F(2), por lo que (z— y UF (y) — FG) < 
(y —x)(F(2) — F()), por lo tanto: 


HN 0-0) 
PEL O 7 Ey 
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El reciproco se prueba haciendo las manipulaciones inversas realizadas en la primera 


HI) - FE)=F0) 


parte de la demostración, es decir, si , entonces f(y) < 


PAN zZ=Y 
z— —x z— — X zZ= 
TY a+ 2 (2) y como 1= A a a 2), con A = e] 
Z—=xX LX Z—X Z—=xX Z—=X 


vi=i=? 
z 


X 
, entonces 
=X 


FAx+(-24)2) <SAf60)+(1-2A)f() Vx,z e (a,b). 


Ahora, variando y podemos obtener que el valor de A recorre todo el intervalo ]0, 1[; así 
f(x) satisface la definición de función convexa. 


Teorema 1.2.13. Si una función f : (a,b) > R es dos veces diferenciable, entonces es 
convexa si y solo si f(x) > 0 para cada x € (a,b). 


Prueba. Dado que f es convexa y es dos veces diferenciable en (a,b), elegimos a < 
x< y <w<z< b y por el teorema 1.2.12 tenemos que: 


FO)=469 - OF) 


y=x ZW 


consideremos los límites como y > x* y w > z7,entonces f'(x) < f'(z) con x,z 
dos puntos arbitrarios tales que x < z, entonces f” es creciente en (a, b), por lo tanto, 
f"(x) > 0 para todo x e (a, b). 

El recíproco se obtiene recorriendo el camino inverso de la prueba de la primera parte. 


Veamos ahora que la función Gamma es log-convexa, para ello necesitamos probar las 
desigualdades de Young y de Holder. 


Inicialmente probemos la siguiente proposición: 


Proposición 1.2.14. (Desigualdad de Young) Sea y : [0, co[> R una función continua, y 
estrictamente creciente y tal que lim p(u) = 00, p(0) = 0. Sean y = p”! y definamos 
UT—>00 
Xx 


p(x) := / pl(u)du, W(x) := 1 v(u)du, entonces para todo a,b € [0,oo| vale: 
0 0 

ab < p(a) + W(b). La igualdad vale si y solo si b = p(a). 

Prueba. Es claro que H(a) + Wíp(a)) = f, p(u)du + Jue y (u)du = ap(a), donde a 


y g(a) son lados del rectángulo (ver figura 1.2). De esta forma, p(a) + W(b) = ap(a) + 
Y (b) - W(p(a)). 
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Figura 1.2: Interpretación geométrica de la desigualdad de Young 


Los siguientes gráficos representan las áreas ab, p(a) y W(b) 


Figura 1.3: Interpretación geométrica de la desigualdad de Young 


En ambos se verifica que ab < p(a) + W(b) y que la igualdad vale si b = g(a). 


Sean ahora p, q números tales que p > 1 y 7 + 7 = 1, por ejemplo, q = 47. 


Corolario 1.2.15. Sean p > 1,a >0yb > 0, entonces 
a? p1 
+ SEACIAr 
Pp q 
y son iguales si y solo si a? = b%. 


Prueba. Sea y : [0,oo[> [0,oo[ la función p(x) = x””!. En este caso su inversa es 
y : [0, oo[> [0, col, y(x) = "Y x 


o = [Cocoa [ras La 


0 


b 1 
vo) =/ vz du = — 
6 PP 
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De acuerdo a la desigualdad de Young tenemos: 


ap  p1 
ab == 
p q 
como queríamos probar. También, la igualdad ocurre si y solo si b = p(a) = a?”!, esto 


p 
es bPT= a, luego b% = DPI = (07) =P. 


Ahora estamos en condiciones de probar: 


Teorema 1.2.16. (Desigualdad de Hólder) Sean p y q números positivos tal que > + m7 0 


l, entonces para cualquier par de funciones integrables f,g : (a,b) > R tenemos la 
desigualdad: 


b b a b a 
/ sreoscoras < (/ ES (Y colas] 


1 l 

Prueba. Sea a = VU SG01? ax) "yB= (e leo)? dx) * sia = 0 quiere decir 
que | f(x)]? = 0 para casi todo x Ela, b| y eso implica que f(x) = 0 para casi todo 
x Ela, b[. Así, para toda función integrable g(x)se cumplirá f(x)g(x) = 0 para casi todo 
x Ela, b| y, por tanto e |f60)gG0)1| dx = 0. Asi, la desigualdad pasa a ser una igualdad 
trivial, 0 = 0, lo mismo ocurre si $ = 0. 

Si ocurre que 4 = +00 0 f = +00 entonces la desigualdad se cumple de manera 
trivial. 

Consideremos el caso en que « 4% 0 y 6 X 0 y denotemos por u(x) y v(x) a las 


funciones u(x) = LEGAL y v(x) = Lera Tenemos que: 
b 1 p? 
/ urd= | FI dx =1 (.1) 
a a? Ja 
b Lp? 
vldx = 7] lgl2dx =1 (1.2) 
f, 00 qa 


ahora aplicamos el corolario 1.2.15 a u(x) y v(x),Vx € (a,b) y tenemos: 


GA" y UE 
Pp q 


u(x)u(x) < 


integrando con respecto a x en (a, b) tenemos: 


b b b 
/ u(uobodx< | rf CA, 
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y aplicando las igualdades (1.1) y (1.2) llegamos a: 


Ñ Lo 1 
u(x)v(x)dx <=+->==1 
a p q 
esto es: A 
Sa 100300) ldx ¡ 
0% IG01? dx)? Va lg ()/1 dx)" 
luego: 


b b a b nn 
/ sensores < (| ES (/ lsco¡t4s) 


como señala la desigualdad. 


Tenemos ahora el siguiente resultado: 
Teorema 1.2.17. I” : (0,00) > (0, 00) es log-convexa. 


Prueba. Seal < p < oo y q tales que: 5 + 5 = 1. De la definición 1.2.1 de la función 
Gamma, tenemos que: 


X % E y 
r(2+2)-/ (2 +p7197t gy 
Pp q 0 
00 , 
=/ RC AL 
0 


sl 1 1 
=/ (PA E E 
0 


1 1 
y usando la desigualdad de Hólder con f(x) = (1*71e71)? y g(y) = (17le7*) % tene- 
mos que 
00 1 1 00 l 00 +] 
/ (41870) P (1971277 (/ ptetd) (/ peta) 
0 0 0 
por lo que: 
x y 1 1 
r(2+2) sir tro. 
p q 
Ahora si aplicamos log a ambos lados de la desigualdad y además hacemos 4 = 5 enton- 
ces 1—-1= a se cumple que: 


n[PAx+(0A4)y»)] <A [PF 60)+ 0-4 Ihn[P60)]. 
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Este resultado es para todo x, y € (0, 00) y por lo tanto la función Gamma es log- convexa, 
de acuerdo a la definición 1.2.3. 


El siguiente resultado prueba que la función Gamma es la única función log-convexa 
que satisface (1) =1 y I(x +1) =x1I (x): 
Teorema 1.2.18. (Bohr-Mollerup) Existe una única función f : R* => R tal que log( f(x)) 
es convexa y cumple con f(1) =1y f(x+1) =xf(x). 


Prueba. Primero consideremos x € (0, 1) y n entero positivo. Usando el teorema 1.2.12, 
con log( f(x)) que es convexa en los intervalos [n,n +1], n+1,n+x+1] y n+1,n+2] 
obtenemos: 


log(f(n + 1D) —log(f0) _ los(fn +x + 1)=log(fM+D) 
n+1-n e M+x+1)-(n+1) de 
PEL + 2)) —log(S (a + 1)) 
MmM+2)-(n+1) 


así: 


(1.3) 


(E >) d log(f(n + x + 1)) — log(f(n + 1)) E (h=»): 
Fu) Xx fum+1) 


Sabemos por hipótesis que f(x + 1) = xf(x) y (1) = 1 de lo que deducimos que: 
fM+1)=n! y además fM+x+1)=(M+x)fM+x)=(M+x)M+x—1)f(n+ 


x=l=...=(M+x)mM+x-=1)-xf (o), entonces en (1.3) tenemos: 
os n! ) < log((n + x)(n + x= 1)---xf(c)) — log(n!) us 
(n —1)! xl 
e ÓR € 22) 


multiplicando esta última desigualdad por x y usando las propiedades de factorial y loga- 
ritmo, obtenemos: 


N+xM+x—1)---xf(x) 
n! 


x log(n) < log < xlogín + 1) 
8 


si restamos x log(n) y aplicando el límite para n > oo, tenemos: 


0< lim [os (E z e: > 23) Ae los(60) | <o0, 


por lo tanto: 


Ji 08 (A) e og) = 0 
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y dado que el logaritmo es continua, se obtiene que: 


fa) = qe (14) 
E M+x)l(m+x-=—1)-- 


Ahora veamos que la identidad (1.4) es valida para x € (0,00). Sea m € Z* tal que 
x Elm,m-+1|[. Sabemos que f(x+1) = xf (x), entonces f(x) = (x—1)(x—2)---(x— 
m) f(x — m) con (x —m) € (0, 1), entonces aplicando (1.4) a f(x — m) tenemos: 


nin*> 
HR E M+x—m(n—1+x=m)---(x—m)' 
entonces 
is (x— Dx —2)---(x—mjnin*o" 


os N+x—m(m-1+x=m)--"x(x—1)---(x—m) 


(1+x)M—1+x)--(n—(m—1) 


si multiplicamos POr 0 m-1+0-0-(m-D7 


= y por propiedades del límite nos da: 


nin* . (n+x)---(n—(m-=1)+x) 
nin? . x x—(m-=— 1) 
por lo tanto: 
Ax) = dat Vx € (0, 00). (1.5) 


Ma M+x)MmM—1+x)-- 


Ahora veamos que 


nin* 


To) = 2171 di = lim 
57 / Ñ n>0 (n+xMm=1+x)--x 


para0<x=<l. 


n 

, Í iñ ; 

Sabemos que e * = lim ( - s) , entonces la función Gamma la podemos reescri- 
n 


n—>00 


[0,0] t n 
bir de la forma I'(x) = / 21. lim ( = s) dt, entonces tenemos que ”(x) = 
0 n 


limines fa (1= py” 71 dt, usando el cambio de variable u = Mos da que: 


1 1 
T(x)= lim d (1—u)y" nuy indu = lim " (1—uu* "dy 
n>00 Jo n—>00 0 
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integrando por partes tenemos: 
(1 —u)u* d 


X 


1 
lim "of (1-—w"u* dy = lim w| 


1 _any1,,x 
>+ 1 E EA du] 
0 Xx 


1 
= lim n” - z/ Au) du 
Xx Jo 


n—>00 


e integrando por partes (n — 1)—veces está última expresión, tenemos que: 


e AS | 1 
T(x)= lim n?- dt y 1 gy 
n—>00 x(x+1)---(x+n-=1) Jo 
n(n—1)---1 1 


= lim n?* 


naco” ED Uen—D) x+n 


entonces 
nin* 


[0,0] 
T(x) = Ple tdt= lim , 1.6 
da le n>00 x(x+1)---(x+nm) (156) 
Por lo tanto f(x) = F(x) para0 < x < 1. Cómo £(1) = F'(1) = 1, tenemos 
la igualdad de las dos funciones para O < x < 1. Las relaciones f(x + 1) = xf(x) 
y Mx +1) = xP (x) para0 < x < 1, implicarán f(x) = P(x) para0 < x < 2 y, 
sucesivamente, tendremos f(x) = T'(x) para cada x > 0, como señalado. 


Nota. Notamos que si —1 < x < 0, entonces de I'(x + 1) = xI'(x) obtenemos: 


TM(x+1) nin*+! 
T(x) = ———— = lim 
n>0 x(x+1)---(x+n)(x+n+1) 
nin? 

= lim . 

n>00 x(x+1)---(x+nm) 
De aquí conclui 1 ión li cali 
e aquí concluimos que la expresión lim — CENA vale para todo x % 


0) LL si 


1.2.2 Representaciones de la Función Gamma 


Desde que Euler empezó a trabajar en la interpolación de los números factoriales, di- 
ferentes matemáticos han representado la función Gamma de diferentes maneras. En esta 
sección exhibiremos algunas de las representaciones más comunes de la función Gamma 
y demostraremos su equivalencia. 
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Observacáo. (Límite de Gauss) De la demostración del Teorema 1.2.18 concluimos que 
para todo x 4 0,—1,-—2,... se verifica: 


] nin* 
TM(x)= lim á 
n>00 x(x + 1)M(x+2)---(x+m) 


Observacáo. Enla sección 1.1, vimos algúnos antecedentes históricos de la función Gam- 
ma y, allí observamos que, al interpolar el factorial de un número, Euler representó esta 
interpolación como un producto que convergía en el infinito, por ello tenemos la siguiente 
representación, que es equivalente a la encontrada en el teorema de Bohr-Mollerup, la que 
dice que para x > 0: 


a n!(n + 1)* 
Ae pes x(x + 1)(x +2)---(x+m) 


que resulta ser una expresión equivalente, en el límite, a la expresión que usaba Carl Frie- 
drich Andrews, Askey y Roy (1999): 
La Constante de Euler-Mascheroni (y) 

Consideremos la siguiente sucesión: 


1 1 1 
Sa =1+2+>5+...+--—In(n). 
23 n 


Sin > oo tenemos lim S, = oo — oo y sabemos que "oo — oo” es una forma indefi- 
n—00 


nida, pudiendo ser convergente o divergente. Calculemos los primeros términos de esta 
sucesión: 


S¡ =1—In(1)=1 


1 
S2 =1+ E In(2) = 0, 8068528194 


11 
$3 =1+ 5) + 3 — In(3) = 0, 7347210447 


$10 =0, 6263831610 
Ssp 0, 5871823229 


S100 0, 5822073317 
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estos cálculos parecen indicar que la sucesión (Sy nen converge. Euler fue el primero en 
e r . 14 n—>00 
probar que la sucesión converge a un límite que llamó y. Esto es: (Sa nen ——> y. Hoy 


este límite es conocido como la constante de Euler-Mascheroni (y = 0,5772156649 ), es 


n 
1 
decir tenemos: y = lim O a In(7)). 
k=1 


Antes de probar que esté límite existe, probaremos algúnos resultados previos: 


1 1 
Lema 1.1. La sucesión ad, = 1+ > +...+ — —In(n + 1) es creciente para todo n € N 
n 


Prueba. Tenemos que: 


1 1 1 
A A) 
esto es: 


1 
An+1 — Aa = —— + In(n + 1) — In(n + 2) a.7) 
n+l 


aplicando el teorema del valor medio a la función y = In(x), entre los puntos (n + 1) y 
by — 
(n + 2), tenemos que LO HO = f'(c) con c € (a, b), entonces: 
=a 
In(n+2)—In(n+1) 1 


M+2=(M+1) 


ce(nm+1,n+2) 


A 


n+0 
1 
In(n +1) —Innm+2)=-==>— ,0€ (1,2) (1.8) 
n+0 
1 1 
de (1.7) y (1.8) tenemos que 4/+1 — Ap = pal + In(n + 1) — In(n + 2) = sm 


a ES 0 € (1,2), ent bt 
= como , , entonces obtenemos que a — ln = 
0 = 


> 0, por lo tanto, 4y+1 — An > 0, entonces 441 > An, por lo tanto, 


+ Da+0) 
(An)nen es creciente, como queríamos probar. 


1 1 
Lema 1.2. La sucesión Sy = 1 + 3 +...+ — — In(11) es decreciente para todo n € N 
n 


Prueba. Tenemos que: 


1 1 1 1 1 
- =1 = ee = —- —] 1)-([1 = Cote =-] 
4150 = 14h AD) (1+3+ +. n(0)) 
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esto es: y 
Sa+1 — 5n = —— + 1n(2) — In(n + 1) (1.9) 
n+1 
aplicando el teorema del valor medio a la función y = In(x), entre los puntos (1) y (n+1), 
F(b) — fa) 
b=a 


tenemos que = f'(c) con c € (a, b), entonces: 


Inn +1)—In(1) — 1 


] 1 
(+ D=() iii 
1 
l 1)— Inn) == —— ,0€(0, 1 
nn +1) In(1) => +9.6(0,1) 
1 
l —1 1)=-— ,0 1 1.1 
n(n) — In(n + 1) e R € (0, 1) (1.10) 
de (1.9) y (1.10)t e +In(n)-In(n+1) = E E 
e (1.9) y (1.10) tenemos que S»+1 we n(m)—In(n STE 


0-1 0-1 
——————— y como 0 e (0, 1), entonces tenemos que S441—Sn = == < 
a+4DO +0)? (0,1 AAA EN 


O, por lo tanto. sh+1 — Sp < 0, entonces Sa+1 < Sn, por lo tanto, (Sn)nen es decreciente, 
como señalado. 


Lema 1.3. Para las sucesiones [an nen y [Sninen, definidas en los lemas 1.1 y 1.2, se 
cumple que Sn > a para todo n € N 


1 1 1 1 
Prueba. Sy > An y aolo st des 1601) + a A 
n n 


1 1 
solo si In(n + 1) — In(n) > Osiysolosiln (14 7) >0ycomo1+-—>1VnenN, 
n n 


tenemos lo anunciado. 


Teorema 1.2.19. (S1)nen y (An Jen, definidas en los lemas 1.1 y 1.2, convergen al mismo 
límite. 


Prueba. Del lema 1.3 tenemos que s, > 4, VYn € N y por los lemas 1.1 y 1.2 sabemos 
que s, es decreciente y a, es creciente, por lo que: 


$1 >8$2>...54 > lp >...d2 > 01] 


con ai = 1—1In(2) = 0,3068528194 y s; = 1 —In(1) = 1 por tanto, s, está acotado 
superiormente por 1 e inferiormente por 0, 3068528194 (0, 3068528194 < sy < 1) .Con 
lo que (Sn)nen es decreciente y acotada, así converge. Sea y = lim s,. 

n—00 


Falta ver que lim a, = y. Para ello notemos qu:e lim (s, — 4/) = lim (In(n + 
n—>00 n—>00 n—>00 


1 
1) — In(n)) = io (6 ( + 3) = In(1) = 0, por lo tanto, ¿HS 0 AN An = Y, 
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como señalado. 


El valor límite y se conoce como la constante de Euler-Mascheroni. usando esta cons- 
tante, Weiertrass obtuvo una expresión para el recíproco de la función Gamma Havil 
(2003): 


Proposición 1.2.20. (Producto de Weiertrass) 
Para x 4 0,—1,-2,-** y y =lima->00 (Dx-1 + — In(n)) se verifica: 


Prueba. Representemos la función Gamma, como el límite de Gauss, dado en la observa- 
ción : 


y nin? ; n?* 
T(x)= lim = lim 4 = 
n>00 x(x + 1)---(x+n) n—>00 Axe 1) (FA) 
n* n* 
= lim = lim 
n>0 y. tn a Qrr2) SÓ ata n>00 xX(1 + 0 + 2. (+ 2) 
podemos escribir n* = e*(MMM-123==55) . ¿+ E+=+5 entonces 
Az O 
lim = lim = 
n>0 x(1+x)- "(14 7) n>00 x(1+x)--(1+ 2) 
= lim ¿A Un) 13... 5) E l a e” E ez oa en = 
n—>00 xx lx d+) (+7) 


= lim e 04-450). tim pare (14 he] = 
n->00 n>00 A 


n 
E l ly . XNY 71 y 
= y 1 limo (1+3+..+3=—In)) . tim TT (+23) pe 
n—>00 ? j 1 
l= 


1 1 
y como vimos lim (1+ 3 +...+——In(1)) = y, con y la constante de Euler-Mascheroni, 
n—>00 n 


por lo tanto: 
=yx o ex 


Xx 
n=1 a sE ») 


(o) = 


y el recíproco de está expresión es conocida como la fórmula de Weiertrass, a saber: 
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Usando el producto de Weiertrass tenemos que la función Gamma se representa por la 


expresión: 
00 ex 
II ara (1.10) 


Representaciones Integrales 


La definición 1.2.1 nos entrega una representación integral de la función Gamma, co- 
nocida como la integral de Euler. Esta no es la única representación en forma de integral 
de esta función. Por ejemplo, la que entrega la observación y que está dada por: 


E 2 
Tx) = 2/ yA" gy. 
0 


En lo que sigue, veremos otras representaciones. 


Proposición 1.2.21. Para a > O tenemos: 
[,0) 
T(x) = af y 179 dy, 
0 


Prueba. Usando la definición 1.2.1 y el cambio de variables £ = ay tenemos el resultado. 


Proposición 1.2.22. 
[0,2] 
T(x) = / eve? dy. 


[o,0] 


Prueba. Es sólo usar la definición 1.2.1 y el cambio de variables £ = e?. 


Proposición 1.2.23. 
ES e 
T(x+1)= Í e? dy. 
0 


Prueba. Nuevamente usamos la definición 1.2.1 y el cambio de variable y = £* se obtiene 
el resultado. 
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1.2.3 Función Beta de Euler 


La función beta se origina de la integral de Euler / x*(1 — x)"dx y que Legendre 
llamó primera integral de Euler. La forma moderna de la función beta es: 
Definición 1.2.4. Para x, y > 0, la función beta se define por 
1 


B(x, y) =/ pa=0 "dt 


Al igual que la función Gamma, es fácil ver la función beta es convergente para x, y > 
O. 


1 
Teorema 1.2.24. La integral | Po X(1—1)7* dt converge para x, y > 0 
0 


Prueba. Como en la función Gamma, separaremos la integral como la suma de dos inte- 
grales y veremos su convergencia por separado: 


1 


É 2 1 
/ pay dt Se pL1— e dt +/ PALA 
0 0 


bl 
=Í + L. 
l x—-1 y—-1 x-1 
Para f¡,0<t< 7 tenemos que -(1—1£) < 21%”, entonces: 
1 1. 
3 , (72 ! 97 _ gx 97x+1 
ES 2/ ldt=2 lim |—| = lim = . 
0 a>0+ | X |g a>0+ Xx Xx 


1 
Para la segunda integral, cuando 3 < t < 1, tenemos: 1471. (1—£)71 < 2(1 —1)7! y 
entonces: 


Ñ 1-1)? 11 
152 /f (1-0 dt =2 lim USO? 
1 a>1> =y 1 
2, Za 
2 (1 -ayY 
=2 lim 22-(1-a” 
a>1—7 y 
97y+1 
A 


Esto establece la convergencia de la función beta para x, y > 0. 
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Como sucede con la función Gamma, la función beta tiene diversas representaciones 
y propiedades: 
5 
Proposición 1.2.25. B(x, y) = 2 Jl sin?*1(0) cos?71(0)d0. 
0 
Prueba. Usando el cambio de variable £ = sin?(9) y reemplazando en la definición 1.2.4 
tenemos: 


T 


B(x, y) = 2/" sin?*71 (9) cos?"1(9)d0. 


pal 


[0,0] 
Proposición 1.2.26. B(x, y) = 0) aya? 


1 

Prueba. Usando el cambio de variable + = —, y reemplazando en la definición 1.2.4 
7) 

tenemos: 


1 x-1 -1 00 
-/ (5) ( _ y du =f ld (1.12) 
oo NU U 1 


y ahora usando otro cambio de variables, u— 1 = z y reemplazando en la ecuación (1.12), 


tenemos: 
00 gy 1 
am 
o (1+2z)4+> 


Proposición 1.2.27. La función beta es simétrica. Esto es: B(x.y) = B(y,x) para todo 
x,y >0. 


1 
Prueba. Dada la función beta B(x, y) = / P=X(1 —£)7* dt y usando el cambio de 
0 


variable u = 1 —f, entonces 
0 
Boy == [109097 du 
1 


1 
=/ uu L(1— uy du 
0 
=B0,x), 


por lo tanto, 


B(x, y) = BO, x). 
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La siguiente proposición expresa una relación fundamental entre las funciones beta y 
Gamma. 


PE)rO) 


Proposición 1.2.28. B(x, y) = En) 
A 


para todo x, y > 0. 
[0,0] [0,0] 
Prueba. Sea I(x)T (y) = / prendo | u* 17 dy, ahora usando los cambios 
0 0 
de variables £ = a? y u = b? y el teorema de Fubini, tenemos que: 
[0,0] [0,0] 2 > 
TOTO) = a/ 0 AAA dad (1.13) 
o Jo 


ahora, haciendo el cambio de variables a = r cos(0) y b = rsin(0), obtenemos: 
Es qa =r? 2x1 0 0 2y-1 
P(O)PO) =4 e" (rcos(0)) (r sin(0)) rd0dr 
o Jo 


5 00 
A cos?%7*(09) sin??=*(0)d0 -2 ) A 
2 0 


[0,0] 
y dado que B(x, y) = B(y, x) y la observación , tenemos que TF (x+ y) = 2/ 271? p2x429=1 gy 
0 
y, en consecuencia, obtenemos: P(x)P (y) = B(x, y) (x + y). Por tanto: 


POr) 


B(x, y) = TE +" 


1.2.4 Fórmula de Reflexión 


Una propiedad interesante de la función Gamma viene dada por la fórmula de reflexión. 
Antes de establecer dicha fórmula, demostraremos la fórmula de Euler-Wallis: 


00 2 
Teorema 1.2.29. sin(x) = x II ( —- Ss ,) 
mn 


n=1 


Prueba. El resultado es cierto si x es múltiplo de 7, dado que ambos miembros de la 
identidad son cero. Supongamos, por tanto, que x no es múltiplo de 7. 

Notemos que si n es un entero positivo impar, entonces sin(nx) es un polinomio en 
sin(x) de grado n y cos(nx) es un polinomio en cos(x) de grado n. Ambos polinomios 
sólo tienen potencias impares de sin(x) y cos(x), respectivamente. 


1.2. Propiedades de la Función Gamma 27 


En efecto: Para n = 3 tenemos: 
sin(3x) =sin(2x + x) 
= sin(2x) cos(x) + cos(2x) sin(x) 
=2sin(x) cos*(x) + (cos?(x) — sin?(x)) sin(x) 
=2sin(x)(1 — sin?(x)) + (1 — 2 sin? (x)) sin(x) 


=3sin(x) — 4sin(x)?; 


cos(3x) =cos(2x + x) 
=cos(2x) cos(x) — sin(2x) sin(x) 


=(co0s? (x) — sin? (x)) cos(x) — 2cos(x) sin? (x) 


=2 cos*(x) — cos(x) — 2 cos(x) + 2 cos*(x) 


=4c0s*(x) — 3 cos(x); 


como anunciado. 


Suponemos, inductivamente, que el resultado es cierto para n = 3,5,...,n—2. En 
particular sin((n—2)x) = Pn-2(sin(x)) y cos((n—2)x) = On-—2(cos(x)) dónde P,-2(x) 
y On-2(x) son polinomios de grado n — 2 en x y con sólo potencias impares de x (sin 
término constante). 


El resultado para n se deduce ahora de la expresión 
sin(nx) = sin((n — 2)x + 2x) = sin((n — 2)x) cos(2x) + cos((n — 2)x) sin(2x) 


a la cual aplicamos las hipótesis inductivas sobre sin((n — 2)x) y cos((n — 2)x). 
Notamos que la función sin(nx), se anula cuando x es un múltiplo de % de modo que, 
cuando n es impar, podemos escribir: 


Hnm—1) 


sin(nx) = K¡ -sin(x) Il (sin?) — sin? (5) : (1.14) 


Esta expresión es así ya que, cómo vimos, sin(n.x) es un polinomio en potencias impares 
de sin(x) que se anula en los múltiplos de F.. La ecuación (1.14) también puede escribirse: 


sin (45) 


5(n—1) inZGo) 
sin(nx) == Ka - sin(x) I ( Da ] ; (1.15) 
r=1 
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En estas expresiones las constantes K¡, K2 son independientes de x, pero pueden depender 
de n. 
Reemplazando en la ecuación (1.15) el valor de nx por x, tenemos que: 


Hn-1) 00 x 
sin(x) = Ka sin (5) TI (1-2) 


r=1 


Así, para los valores de x considerados se tiene 


3(n=1) 
_fsiaco ] FED) 
K2= E al j 1 ( a) : 


r=1 


En esta última identidad el miembro de la izquierda es independiente de x, y por con- 

siguiente, lo mismo debe acontecer al de la derecha. Para determinar su valor común, 

hagamos x > 0, entonces 

sin(x 

K», = = lim ( ) = 
x>0 sin e ) 


La igualdad de la derecha se obtiene aplicando la regla de 1"Hópital. 


Luego, para todos los valores de x considerados y para todos los enteros positivos n, 
se tendrá: 


20D sin? 
sin(x) =n + sin (2) yl f an 


r=1 


] =m- sin (5) fr (14 MA). (1.16) 


r=1 


Aquí, las funciones f; (1) (x) están definidas por: f,(n)(x) = 0, sir > > y fr (mx) = 


_ sin? (%) 


caso r < == 


sin (2%) ' 
Ahora, cómo para 0 < 9 < 5 se cumple que  < sin(9) < 6, obtenemos para n > 211 
que: 
AM dE 


IAE) S 12.2 = 452" 


n Ar27/2 4y2 


j A TIÓX 
Observacáo. De este último resultado, vemos que sin(xx) = xx II ( = 5) en- 


m2n2 
mn) - | ( Ml 5). 
n 


n=1 


n=1 
tonces 
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Ahora podemos probar el siguiente teorema. 


Teorema 1.2.30. Fórmula de reflexión de Euler: 


TGO)PU-=x) = EZ 


RT 
sin(7rx)' 
Prueba. Dado que la función Gamma tiene varias representaciones, usaremos '(x) = 
*p! 
n*n! 


lim , entonces 
n>00 x(x+1)---(x+n) 


Xp! n Xp) 


POyra A a men pe pee e pene er la 


n!-n! 


n 
A AAC AAA 
e [n1]2 a (1-2-3---ny e 
n>00 x(1—x2)1(4—x2)---(n2—x2)  n>00 x(1—x2)-»-(n2 — x2) 
= lim . = lim a 0 
a 
id E 1 de 1 e TT 


n>00 E A (1 E 5) MEE (1 E 5) peón sin(7x) 


Otra manera de probar la fórmula de reflexión de Euler, es utilizando una integración 
de contorno en el plano complejo. 


Prueba. Hagamos 0 < x < 1 y usemos el hecho que B(x, y) = HE y Blx, y) = 


00 pl 
——————dt. Hacemos (x, y) = (x,1 — x) y tenemos que 
il ES Cy) =( AS q 


pl 


(1++) 


PQO)FU-—x) 
TMx+1-x) 


B(x,l—x) = =rra=o= dt. 


Ahora para calcular este integral, lo hacemos mediante la integral de contorno: 1 
c1i=2Z 


30 1. Función Gamma 


Figura 1.4: Contorno para integral compleja. 


Donde el contorno C = C¡ UC» U C3x UC, consiste en dos circunferencias centradas en 
el origen y de radios R y e respectivamente, como muestra la figura 1.4. Por el teorema 


de residuos tenemos: 
¿xl 
dz = 211, 
C l-z 
entonces 


x—-1 x—1 x—1 x-1 
— 21 i =f - de+ | ? d+ | ' do+ [ ] dz (1.17) 
ci 1-2 c, 12 c31=z Ca 1-2 


donde C es la circunferencia exterior de radio R, C2 denota el segmento que va desde 
—R a —e,C3 denota el radio interior e y Ca es el segmento que va de —e a—R. 

Entonces para el circunferencia exterior hacemos el siguiente cambio de variable z = 
Re?? y tenemos: 


x—1 T ex=1 i(x-1)0 Ñ T ¡px ix0 
a y aa 
c1-z ar 1-—Re' 7 1-Re' 


y dado que O < x < 1 tenemos: 


7 ¡px ix0 
L'= lin | 4, 
R>o J_y 1— Rel? 
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entonces 
Tr ¡RXgix0 Tr ¡pr 
[11] = lim / ——5d0|< lim Je d0 
R>00 |)_7 1— Rei? R>o)- 7 R=1 
R*g [7 
= lim 
Roo R=1|_z 
Rx 
=27% lim 
R>o R=—1 
_ 0 
o 


aplicamos 1*Hópital y tenemos: 


11] < 27 lim xR**=0, 
R>o00 


por lo tanto, f1 = 0. 
Para la circunferencia interior, /3, hacemos el siguiente cambio de variable z = ee*? y 
tenemos: 


qe=l xi ¿XL pix 1)0 . —xr ¡ex gix0 
L = / dz = / ——_ AUee?) == il pd, 
c31=z Fi l — ee! nn  l1-—ee! 


entonces 0 
—T XylX 
lee 
3 = lim —d0, 
e>0 Jr 1 — cet 
por lo que: 
—xr ¡EX eix0 —A [gx 
|13| = lim / ——— d0| < lim le d0 = 
e>0 Jr  1—eei0 e>0Jxr  e€-1 
ep” 
= lim 
e>0€—=1 |, 
Xx 
€ 
=-—2x lim 
e>0E€—= 


por lo tanto, f3 = 0. 

Para el segmento que va desde la circunferencia exterior a la circunferencia interior y para 
el segmento que va desde la circunferencia interior a la circunferencia exterior, Ca y Ca, 
escribimos z =—1 = te"! y z =-—t= te” respectivamente. Entonces reemplazando 
en (1.17) tenemos: 


0 ¿x—1,ixx 00 ¿x=1,—ixTT 
Í e Í e 
1 | A 
o 1+t 0 1++1 
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entonces 
00 pl . 
1) dt [pa Ed poa] , 
o 1+t 
por lo que: 
Na 27 1 
0 ++ e xx ecixa 
a xi_p—xi 
y dado que sin(x) = “5 *” entonces 


[9,2] (pel 
/ d= E, 
o 1+1% sin(7x) 


POO)FÚ-x)= 


por lo tanto: 


sin(rx)' 


Proposición 1.2.31. (Fórmula de duplicación para la función Gamma): 


221 Tx) TP (x + 3) 
WE 


Prueba. Dada la definición de la función beta 1.2.4 y la proposición 1.2.28 tenemos la 
siguiente relación: 


rOx) == 


Mac A 
BP) = “ox > 


1 
/ td=0'*d (1.18) 
0 
y E u+l 
y usando el cambio de variable £ = uE tenemos: 


1 1 
| a e y du == aii | a NN u2) du 
-1 0 
ahora, usando el cambio de variable u = y/u tenemos que: 


1 
1 
B(x,x) =a va - vw) dv = IBA) 


a O) _ AYTO) 
T(x +3) T(x +3) 


(1.19) 
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y dado que las ecuaciones (1.18) y (1.19) son iguales, entonces: 


2 ES. TE 
Ta+b  TQOx)” 


por lo tanto: 


221 Tx) TP (x + 3) 


PQx) = 7 


1.2.5 Aproximación de Stirling-De Moivre 

Para el cálculo aproximado del valor de la función Gamma podemos usar la aproxima- 
ción de Stirling-De Moivre. 
Proposición 1.2.32. Sea x > 0, entonces: 


1 139 571 
288x2  51840x3  2488320x* 


1 
TQ) = Ze + a + + 0(1 7) 
6 


Prueba. Sea x > 0, entonces usando la definición de la función Gamma 1.2.1, tenemos 
00 


Tx +1) = / t*e* dt. Haciendo el cambio de variable 1 = x(1 + 4), tenemos: 
0 


ad [,0) 
A =/ (1 + ye dy = / x*(1+uye *e Y xdu 
=i 5 
e [,0) 
=/ yA LX ¿0 px In(1 48) q, — di 3 | ¿e u+In(1+u)) q. 


Expandiendo la función In(1 + u) en serie de Taylor centrada en 0, tenemos que In(1 + 
00 k+1,,k 
DY” u 


u)= y e Y luego: 


k=1 


(¡=p + k u? 00 (=D*+1,k 
—u + In(1 +4) =-u + Y ra =-4>+u e 


k 
k=1 k=3 
k+1 uk 
y? 1 
e 
para |u| < 1. 
Como queremos que la aprostmación final, de la serie, contenga términos de hasta el orden 


1 k+1 uk . 
75 truncaremos la serie . ED hasta el orden 10, esto es garantía suficiente 


k=3 
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para obtener los primeros 5 términos de la serie final. Por lo tanto, truncaremos como: 
(- 1 an k 
0(u) = y 
entonces 


00 2 oo 2 
Mx+1)= aa] AG) y a Pl EE pro gy, (1.20) 
-1 -1 


Por otra parte también podemos expandir e**(% como serie de Taylor para x > 1, es 
decir, 


iú _ > (xo (u))> Ñ > xo 10 (=D 


1 ] 
ao v! A eee k 


e 


truncaremos esta serie manteniendo los primeros 8 términos para asegurar así los 5 prime- 
ros términos de nuestra aproximación, es decir, 


Para expresar estos términos, al desarrollar estas sumas finitas, usaremos la fórmula de 
Leibniz, o sea: 


n! 
(ar +42+...+ap)” = Y —_—_—_—— a a3o ae. 
01 !07!- 
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Así, tenemos que: 


v 


8 10 A 10 k+1,k 
¿Ao (u) y En de y ED A $ ED 
v! 
k=3 


+ 
v=0 k k=3 k 
a 10 (14H 1 ya fp 10 (1 +1yk 
2! k 8! k 
k=3 k=3 
>> y 10 a ¿10? 
=1 o 
E qe TY A E 10) * 
3 4 
a % ¿10 si A E y 10 a 
6 3 4 10 24 3 10 
5 6 
a ¿10 E ea ¿10 A 
1201 3 10 720 3 10 
7 8 
Si e ¿10 A a Joy 
5040 Y 3 sl 10 40320 Y 3 ób 10 
Sa 8 de 15 Ns ¿10 , 
= Xx 
4 5 6 7 9 10 
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a e e, cis A mi 
192 144 225 180 216” 


ñ y0 u 18 es ee ul? A 20 A df 20 a ñ 
7201729 *** 972 1215 1296 


xa” u2l u22 x8 u24 
ll rn ia dd esmas EE aos rl da 
+5040 2187 ON 2916 Ñ + + 40320 | 6561 y 


reordenando los términos anteriores tenemos: 


ES 472 x 87, x? 19x2 de x 9 x? 15317 
== ->]|u U 
480 6 162 180 “9 72 1400 
ES a) 10, e El PT E 4 49317 
10)” 1440 280 )” 1944 — 17280 
A a IAN 77x% 
+ U u + E 
30240 648 50400 25920 
4049: 1, a > TA e y 
q u u 
100800 29160 155520 7776 


19026Íx* Y 16. 23% 49 ed 503%" 1.38 
PRO Por u v) 
2903400 77760 524880 933120 


E 19 107x6 , 20 e ¡2 dl 27 
116640 4665600 11022480 2099520 
8 
A: 
264539520 
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1 


dicci 2 y? A 
ahora usando la sustitución u = (5) w, en la ecuación (1.20) tenemos: 
X 


00 ; 
Tx+1) a | e" DP ¿ol gy = 
-1 


por lo tanto, 


[0,0] 
T(x+1= Viptdes | ¿70? ¿HA() q yy, 
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(1.21) 


aquí hemos usado el cambio de variables o (u) > A(w) y, dado que la serie o (u) converge 


A 


para lu] < 1 y como u = ( w, entonces la expansión en serie será más precisa si 
x 


x > 1. Además, si sustituimos en el límite inferior de la integral el valor —oo por el valor 


1 
Í=-— (5) E tenemos, reemplazando en la ecuación (1.21) que: 


[0,0] 
T(x+1) Vatios | ¿70? ¿XA(W) qyy = 
—00 


[o/o] 4 
E | Pci fa + de wé+ 
66 x 9x 3x2 
(A 2 a (4 2, 16 o 
30x2  x3 9x2  175x3 5x4 


(3 493, 6698 Y 12, (_154 
243x2  270x3 * 945x1 405x3 


8098 Y 14 8 190261 Y 16, 
1575x4 243x3  113400x* 


32 1006 Y 1 428 
A ——_—_——— Y pa 
3280513  3645x* 18225x4 
32 99 32 24 
320 eri5t 


40 39 
+ y M;(x)wY + E MAA dw. 
j=13 j=1 


(1.22) 
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Veamos ahora como son las integrales de la forma: 


G,(0) := / "TU yy = / f(t,n,a)dt V(n,a) € ZE xR?. 


[o,0) 
% at? . 
Para n = 0 tenemos que la ecuación Go(a) = / e 7% dí que es la integral de Gauss, 
0 


TT 
entonces Gp(a) = 7" 


3 
[0,0] 
Para n = 1 tenemos que G; (a) = / te? gy e integrando por partes concluimos que: 
0 
G¡(a) = — 
Ahora, derivando la expresión de G, (0) tenemos: 


E (e) 


[o,0) 

SN £(t,n,a)dt = -/ (1+22242 41 = —Gn+o(a) 
0 

y cambiando índices: n > n — 2, tenemos: 


G, (01) = —— Gn-a(0) (1.23) 


1 1 
Iterando (1.23) y usando que Go(a) = 3 dE Gi¡(a) = — y Fesulta: 


» Paran par: 


1 _1 VIT 
G2(a) q Go(a) da (5 TO 2) E 
d d (1 _3 1.3.7 
Ga(a) = da R(2) d (qua 1) > 
d (1-3 1.3.-5- 
Gola) = ———Ga(a) = 5 Sy = 5 va 
d q2.23 a -.2% 


d 1-3-5---(n=1)- yr 
Gnla) == 7 Gn-2(0) = > = 23 
a. 
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1 1.3.--(n—1 
y dado el hecho que /” € > ) = z ) eN por lo tanto, 
2, 


1 n+l1 
G, (a) = rl > ) 


» Paran = 2k + l impar: 


d eta Í 4 1 
Gs(a) = Gta) == (2) a (53) de 


d d 1 1 2 2 
A A (5) de E =45) = 205 


d d 2 1 6 6 
AS (55) a E el) = 24 


1:-2-3-=k _ kl. 
Da Kk+1 2 1? 


d 
Gor+1(0) = = yy ao aa) = 


dado quen = 2k + 1, entonces k = | _ ml y puesto que k! = P(k+1) 
tenemos: 


n+1 
G, (0%) = TA) 


n+1 
aq 2 


Notamos que esta última ecuación es la misma para n par o impar. Como en la ecuación 
(1.22) el intervalo de integración es todo R y dado que la función f(t,n, a) es par o impar 
(va a depender del termino 1”), resulta, por simetría, que 


1 
SS r(15 ) . 
J p”gat? 21” SInes par 
a 2 
e 0 si n es impar. 


Se observa claramente que los términos en la integración de 


50 40 
1) eu? y M; (x)w2 dw entregarán un valor bastante pequeño. Notamos, además 
7% j=13 
P EZ 
que | en E Nj(x)w2+1dw = 0, Vj, dado el hecho que 2j + 1 es impar. Por 
ca Pra 
tanto, en la ecuación (1.22), tenemos: 
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1 [0,0 [0 
T(x+1) 421% +20 * / e” dw— =/ wit dy+ 
—00 


—00 


4 4 ss 32 En 
ell Dreta: " ¿10? guy, + om) 24? guy 
9 3x2). 2066715x* Jo, 


1 3 4 4 15 
= Verte | a o E ( ) e dE 


4 9x 3x2 8 
2/7 A |1-5+ 


EEN 4096 4x 
4 4115 47 2 1105 4 
+ + 
(5 2) 8 (5 25) 16 (5 
612 16 Y 945 8 493 6698 
5 5=) 32 (2 0x3 + 22.7): 
10395 154 8098 Y 135135 8 
64 (a 511) 128 Cn 
190261 y 2027025 32 1006 Y 34459425 
m0) 265 (osa 55.1) sa. * 
428 654729075 32 13749310575 
+Te225x*' 1024 328054 2048 
32 316234143225 
+ 2066715x* 4096 | 


y simplificando algebraicamente obtenemos: 


' 1 4 139 571 
T 1) = V2ax*+2e Y [1 . 
o peSeS ( + Tax * 288x2 — 5184013 aa) 


Y usando la relación de recurrencia de la función Gamma PP (x + 1) = xP (x), tenemos 
que: 


1 1 1 139 571 
T(x) = V21x*72e*(1 O(x 5), 
e) te O at aegx2  51840x3  2asea20t IA Y 
siendo O(x7*) el error estimado. 
De esta aproximación concluimos que: 
1 1 1 139 571 
T(x) = V21x*72e*(1 O(x 3). 
E e ra s1840x3 7 2a8g320xt + 007) 


2x  288x2 


1.2. Propiedades de la Función Gamma 41 


Observacáo. Usando la última expresión con x =n € N grande tenemos: 


Tm) = V2an"> 


1.2.6 Fórmula de Multiplicación de Gauss 


Otra propiedad de la función Gamma, que es de utilidad, se resume en la fórmula de 
multiplicación de Gauss: 


Teorema 1.2.33. 
E k cálo m 
Tr (++2) = (21) 2 n2 ""Tí(nx) Vxé¿ men). 
Se n n 
Prueba. Usando la fórmula de recurrencia de la función Gamma tenemos que: 
k k k k 
Plx+-)])=P|x+--1+1]=[(x4+--—1]P|[x4+--—1)]. 
n n n n 
Por otro lado, usando el límite de Gauss sabemos que 


mim? 
T == lim , 
qUe m->00 x(x+1)---(x+m) 


entonces: 


y usando la aproximación de Stirling-De Moivre, tenemos: 


=T(m+1) =V 27 (m + 14227 (M+D 
1 1 
=V2m"+2 ( + 7) ( + -) en 
m 


m 


por lo que: 


1 A 
Vlam"+2g Mp X+571 


r (+2) = lim, ER m = ++ 
n ES ES 2). .(uetks nm— 2) 
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de donde: 
k y Vlamre"mx+R=3 pm 
Plíx+=)|= lim a 
n m>00 (nx +k)Mnx+k+n)---(nx+k+nm-—n) 


Ahora, si multiplicamos estas expresiones desde k = 0 hasta k = n — 1 tenemos: 


n—1 n 
k Mm" ¿29M y nx > 
Tr (++) = a, Eme ton 


(Mx)nx +1)---(nx + nm_-—n) ) 


n n=1 k 
m"*=3m2k=0 hp"? ] 
, 


(ax+D)--(Mx+n—1+nm—n) 


esto es: 


a k A (v21)" (mayor ¿27 yr 4 Elo 
II Plx+-=)]= li 
k=0 de 


m ] 
m>00 (nx)(nx + 1)---.-- (mx+n—1+nm-=n) 
n—-1 
Además, dado que yo — es una suma de términos en progresión aritmética, tenemos 
n 
k=0 
-1 
A 
7 = = —— Y, entonces: 
n 2 
k=0 
1-5 py 5 


jane k , (V27)"(mny"*¿ 292 py? 

II Plx+=)]= lim , 

PA n m>00 (nx)(nx +1)---(nx+n-—1+nm-=n) 
por lo que: 


1 


TI P (s + s) =,li (Y mie pz 
k=0 ás 


im ¿ 
m>00 (nx)(nx + 1)---(nx+n-—1+nm-=n) 


Ahora, haciendo la sustitución mn —> m, obtenemos: 


¡e k ] (VJ mtm > 


n-1 11 

k 2 n m__m nXx—=> ¿nx 
T] (s s) hi (21 m"e"m n 
k=0 n 


— m>00 Mxjnx +1)---(nx +m-= 1) 
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n-1 k . ( lazy MNirmn? 2é M yy x az nx 
II Plx+-=)]= lim 
le n m—>00 Mx)(nx +1)---(nx+m-— 1) 


j del 
como lim V21m”*2e7" = ml, entonces: 
m-—>00 ? 


n-1 k ( 2 mim" 1 y 30x 
yl Plx+ = lim 
Pao n m>0o0 (nx)(nx + 1)---(nx+m-=1) 


mim"*=1 


an llax y 
>= (27) 2 n2 lim 
m>00 nx(nx + 1)---(nx+m-= 1) 


>= Cm n20n im decai : de 
m>0o0 nx(nx +1)---(nx+m-=1) (nx—1) 
mim"x=1 


>= COn? (mx —1) lim ; 
m>00 (nx — 1nx(nx + 1)---(nx +m-= 1) 


Usando, nuevamente, el límite de Gauss, tenemos: 


cit k a n-=1 1_ 
II PTlx+=)]= lim (lr) 2 n2"(nx— 1) (nx—D), 
n m>—>>00 
k=0 
por lo tanto, 
n—-1 k A 
II Tr 6 + s) = (27) 2 n2 "I(nx) 
n 
k=0 


m ; 
para todo x É (—; m € N), como anunciado. 
n 


1.2.7 Función Digamma 
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Una herramienta para calcular la derivada de la función Gamma es la función digamma 


o también llamada función psi Artin (1964). 


Definición 1.2.5. Se llama función digamma, denotada por W(x), a la función W(x) = 


TP'(x) 
T() 


, parax 40,-1,-2,.... 


De la fórmula de Weiertrass: Fa 
natural en ambos lados tenemos: 


[o,0] 


a y [Im (143) -<]; 


1 x . A 
) == Es (1 + =) e», aplicando logaritmo 
y n 
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derivando esta expresión, con respecto a x, tenemos: 


To 1 


1 1 2 1 1 
"TO her El +5 A 


entonces 


E ES PERO O E OO 
a XxX MES +) 


w(x) = El Eos a 
X AAA ¡canas 
Y n A x+n-1 n+x 


n 


Y (x) 4 Y +, lim Y, : 

= =y--— == —_——— im ———- — 

dl Y A NA x+n-1 PAE l x+n-1 n+x 
En 7 Y EE RSS 


telescópica 


. 1 1 
Ea os ma) li) 


por lo tanto, 


[oe] 


x—1 
200 ++ El E 


En el cálculo anterior hemos supuesto que I'(x) > 0 (para tomar el logaritmo). Es 
fácil verificar que la misma expresión es valida para el caso en que I'(x) < 0. En efecto, 
en esté caso usamos la expresión: 


Lee] (+ er 
== = (—X)e = Je: 
MES) n 
n=1 
Por tanto, el cálculo vale para x 4 0,—1,-2,... 
Al igual que la función Gamma, la función digamma tiene varias propiedades intere- 


santes, como su propia fórmula de recurrencia, diferentes representaciones integrales, fór- 
mula de reflexión, etc. 


Proposición 1.2.34. (Fórmula de recurrencia para la función digamma): Para x A 
0,—1,-—2,... tenemos que : 


W(x +1) = W(x) + a 
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Prueba. Ocupamos la fórmula de recurrencia de la función Gamma: 


T(x+1)=x1(x), 


aplicando logaritmo natural en ambos lados de la igualdad y luego derivando con respecto 
a x, tenemos que: 

C++ Dd 

FEED 


por lo tanto, 
1 
W(x +1) = W(x) + 7 


Proposición 1.2.35. (Fórmula de reflexión para la función digamma): Para x € 7 tene- 
mos que : 
W(l—x) —- W(x) =  cot(7x). 


Prueba. Ocupamos la fórmula de reflexión de Euler para la función Gamma: 


TOTÚ-x)= a 


aplicando logaritmo natural en ambos lados de la igualdad y luego derivando con respecto 


a x, tenemos que: 
Po) TU(l=x) xcos(x) 


TG) TU-=x) sin(rx)” 


por lo tanto, 
W(l—x) — W(x) = rm cot(7x). 


Proposición 1.2.36. (Fórmula de duplicación para la función digamma): Para x 
0,—1,-2,... tenemos que : 


2W(2x) = 21In(2) + W(x) + Y (s + 5) y 


Prueba. Ocupamos la fórmula de duplicación para la función Gamma: 


2% PT (+ 3) 
Ja j 


aplicando logaritmo natural en ambos lados de la igualdad y luego derivando con respecto 


rOx)= 
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a x, tenemos que: 


2 


P'Ox) M6), +3) 
2 O PEN 


por lo tanto, 


24(2x) = 21n(2) + W(x) + Ww (s du 5) . 


Proposición 1.2.37. (Fórmula de multiplicación para la función digamma): Para x 
0, —1,-—2,... se cumple que 


n—-1 


YV(nx) = In(n) + > Dv (+ 5). 


k=0 


Prueba. De la fórmula de multiplicación de Gauss para la función Gamma tenemos: 


n—-1 
—n k 
T(nx) = Or) pro Il Tr (s + s) , 
k=0 de 


aplicando logaritmo natural en ambos lados de la igualdad y luego derivando con respecto 
a x, obtenemos que: 


n-1 [" pk 
na - =nIn(n) + nta) el 


a 
+ 
3| 
a 


por lo tanto, 
n—-1 


Y(nx)= In(n) + : Dv (+ 5). 


k=0 


Otras Representaciones de la Función Digamma 


La función digamma, al igual que la función Gamma, tiene otras representaciones, 
como las siguientes: 


Teorema 1.2.38. Para x > 0, tenemos: 


a) W(x) = e z (e - +») dt (Dirichlet) 
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b) W(x) = IN (F- ya (Gauss) 


Prueba. a) Evaluamos la integral doble: 


pj] rdidx= | | oraras. 
1 Jo o Ji 


Evaluando la integral del lado izquierdo tenemos: 
Ss [o,0) Ss w 
18) eddi= lim, / e *dtdx 
1 Jo 1 W>% Jo 
s tx w 
[Jen 
1 wW>00 Xx 0 
s ex] (1.24) 
=/ | lim (- + ) dx 
1 wW—>00 Xx XxX 


=f Aito 
1 X 


1 
Evaluando la integral del lado derecho: 
00 ps 00 e s 
/ / e *dtdx a (- Jas 
o 41 0 1 
[o,0) —SX A 
a 
0 Xx Xx 


Usando la igualdad (1.24), en la integral anterior, tenemos: 


/ j Erardi= | ia A (1.25) 
o ¿1 0 id 


Ahora, usando (1.25), calculamos la integral doble: 


- es (+1) 
IN O , dl sdt. 


Para ello, integramos primero con respecto a £: 


- es(+1) 
P a —————————dtds= 
—H _ ¿ost 
a E mies (E —— ) aras 
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Ss 


= In(s) 
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y por (1.25), tenemos que: 


—t —st 00 
de ES =1 6 E ) aras sl 1078 In(s)ds =P (x). (1.26) 
t 0 


Ahora integrando la integral doble con respecto a s tenemos: 


00 00 ¿x=1,=5 
I E] / AS (e* —e *Ddsdt = 
o Jo Í 
se 1 ers est 
le a El pa? a Mo f"* ————— —dsdt. 
A 


[,0) 
En /; tenemos que: / s 1078 ds = T'(x); y para evaluar 1, usamos el cambio de 


0 
variable s(1 + 1) = u. Por tanto: 


1 


oo pt a U AOS e 
I =r 6) / ES -f ' (> ña 5) e 1+0 E Jn 
[0,0 e! [0,0] T(x) 
ll Er! DE (1.27) 


1 
=r ¿(e (+ E) E 


y como (1.26) y (1.27) son iguales, entonces: 


1 1 
Tx) = re) / £ (e - +") dt, 


ÓN ia 3 APR: 
a =b eq) 


por lo tanto, 
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b) De a) tenemos que: 
291 1 
v)=/ A 
0 Í a + t)* 
| 1 
= lim ll — (e =- a+) dt 
e>0+ Jin(i+e) 1 (+12 


o os e“ du 
esot nato 1 n(1r+e) (e” — De 
para la segunda integral usamos el cambio de variable 1 + £ = e”, entonces 


00 et 00 ee du 
W(x)= lim (/ — dt - pez) 
e>ot Win(1+e) * In(1+e) (e% — 1)eX 
[0,0) e. [0,0) e "YX dy 
= lim il —du— | E 
e>0+ AJin(1+e) In(1+e) (e% — 1) + ¿xr 
00 —u 00 =ux 
= lim (/ el só — du) 
In(14e) 4 In(1+e) 1= e 


e>07 
—uUzx 
= lim yA (E Es )«.] 
e>s0t LÍn(1+eo A Y Lee 


por lo tanto, 


Aproximación de la Función Digamma 


Al igual que con la aproximación de Stirling-De Moivre de la función Gamma, se 
puede encontrar una aproximación asintótica de la función digamma. 

Para esto introduciremos los polinomios de Bernoulli, a saber: la sucesión (B, (x) Zo 
definida por las condiciones: 


* Bolx)=1. 


* B(x) =nB),_¡(x),n > 1y 
1 
5 Br(x)dx =0,n =0,1,2,... 
0 


Notamos que B¡(x) = 1 y entonces Bi(x) = x + c. La condición de Bi(x)dx =0 
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1 
implica Si + c)jdx =0y (5 + E = 0. Esto es, l +c=00c= 1 Por tanto, 


Bi(x) =x-— z. 


Procediendo de manera similar para n = 2,3, 4,5, 6 tenemos: 


1 3 Xx 
BAx) == x + Ba) = 0-0 
1 5 5 Xx 
B pa 4_3 3 Po B = Cs 3 
4(x) =x x+x 30 s(x) = x? qe +34 6 y 
5 1 1 
B = 623 5 o. se 
s(x) =x Xx e ¿A Eder 


Observamos que B27+1(0) = B2n+1(1) = 0, para todo n > O y que B2,(0) = Ban(1) = 
B»,, para todo n > 1. 


Usando estos polinomios calcularemos, dada una función diferenciable f(x), el valor 
1 
S. F6ddx. 


Tenemos e Fooddx = de FG0)Bolx)dx = de FG) B| (x)dx. Integrando por partes 
tenemos: 


1 1 
f reoax= [rosca 
0 0 
1 
= OB [Or 
0 
1 Doe 
=340+F0)=5 [| £ 0804 
1 
= 1040093 | £608200(4= [| 100820041 
0 
=40+F0)=3 [0-10] 8,00)+ 
1 
> +3, F"69-3BsC9dx 


20 ( 


5 540) +3 (0))= FO0-F0)+ 


e E FU)BLODdx 
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1 
=UO0+10)- HL 0-10)+ [108401 
1 
E y a 


1 
=UO0+ 0) (FO F0)= 5h F"C)BsCIdx 


woo Ero- 10-74 / 10800 


= 040) -Z(L0-F0)-3 al - 
1 
= [| 149 coracoar] 
l B» / / Ba mM mM 
= SOHO) E LO-F0)-0"0-F"0)+ 
“a "ab FOOBL(dx 
5 FOO) ZO 0)- 0770) + 
+2 aa -f a 
= UO+IO)- HL) (10 -F"00)- 
Ñ EN aaa 
= UO+IO)- HL O-F0)- (10 -F"00)- 
7 somo) S í romana 
=UO+F0)-HLO-F0)- (10) -$"0)- 


ES j y (vi) 
EOS e NE ()Bo(x)dx 


+72 
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y asi sucesivamente para obtener la relación: 


n 


1 BT pora) Qk=1) 
ll 104: =10+101- E 0120] 


1 
AT 


1 Band. 
0 


Ahora si reemplazamos f(x) por f(x + 1) tenemos he Ax + Ddx = E fu)du. Por 
tanto, podemos hacer el cálculo anterior con g(x) = f(x + 1) para obtener: 


; 1 AB 
dl gjdx = 180 +201- 2 [0-0] + 
k=1 ¿ 


; Qn) 
+ / ¿OD ) Ban(a)dx. 


esto es: 


2 
f roa=¿140+f01- La E [f2D)- 0] + 


Qk)! 
(+ 1)B 

tí PAU + 1) Baníxddx. 

A partir de los cálculos para de dx y fi F(x)dx tenemos 


; dx = d ; d Pe 0 1 . 2 
[rwoax= [| 1604x+ [sa =10+ 10+3/0- 


en 
Esla rmvoj: 


+an a FDW +1) Ban(x)dx. 
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Podemos continuar para calcular [7”_, f(x)dx como E F(x + m)jdx y obtenemos 


[rod 10+10+10+...3£00- 
E [10m — 1240) + 


O e 
+2) BA e Mx + k)Godx. 


Esta última ecuación recibe el nombre de fórmula de suma de Euler-MacLaurin Fer- 


nandez (2008). 


1 
Aplicamos lo anterior a la función: f(x) = a y tenemos que: 
zZ+HX 


Li a 1 


/ sood:= || ir be 


Asi, usando la fórmula de la suma de Euler-MacLaurin, tenemos que: 


1 1 1 1 
=/0 + £0+/0)+ fm D + /0m)- 


Z zZ+E+m 


=— Bor S E 
La (LID (mm) — FAIDO)) + 


1 n—-1 
) Baní) LO (+ dx = 
9 k=0 


1 
Cn)! 


1 És 1 
(2 +m- 1)? 


1 1 
e Ey (2 +m- 1) 


2 (Ck)! ( +myk+1 AREA ) + 


1 n—-1 
[no 16 + bar. 
9 k=0 


+ 
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Considerando n > m y haciendo n > oo tenemos: 
Lo 1 O 2 Bon 
UN 22 (z +n) 2 pe 


Ahora podemos establecer el siguiente resultado: 


Proposición 1.2.39. (Aproximación de la función digamma) Para z 4 0,—1,—1,-—3,+*- 
tenemos: 


1 1 1 
W(z) = In(z) — 8 + 


-6 
1222 + 2407. OE. 


Prueba. De la sección 1.2.7, consideramos la expresión de la función digamma: 


is +2 ==...) 


Derivando, con respecto a z, tenemos que: 


[o,0] 


1 Esa 1 
yw (2) =Y' ——— ent Wi+n= Y — == 
(2) 2 Son=0" entonces Y (z + 1) 2 Era? 


De esta forma, tenemos: 


z 222 (z+n) z2n+1 
n=1 n=1 
[0,0] 


1 1 , Bon 
> =p Pere an Tr? 
n=1 


1 


Bon 
32 +2 2941 e integrando esta última expresión con 


1 
entonces W'(z + 1) = — 
z 


respecto a z, nos da que: 


1 B2n 
y 1)=1 —- Ñ 
(+1) nz) + >, y +C; 
n=1 
Se puede probar (ver Arfken y Weber (2001) ) que C¡ = 0, por lo que: 


1 ES 
W(2+1)=1 o e 
(+D 2) 2z 2 2nz21 
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Así, usando la fórmula de recurrencia de la función digamma: 
1 
Wíz+1)=VW()+ e 


obtenemos que: 


1 Bon 
(2) =In(2)- 72 - y 
n=1 


1 1 
y considerando los números de Bernoulli B2 = 6 y Ba = 30 tenemos: 


A 1 
YERMO=>= 2 zas 


POE), 


que es la aproximación que queríamos encontrar para la función digamma. 


1.2.8 Función Polygamma 


Las diversas derivadas de la función digamma generan funciones conocidas como las 
funciones polygamma. 


Definición 1.2.6. La n-ésima derivada de la función digamma para x 4 0,—1,-2,... 


es: 
n 


y (x) = Y, (x) = Á 


Epa W(x). 


Al igual que la función Gamma, digamma, esta función tiene varias propiedades in- 
teresantes. 


Proposición 1.2.40. Para x 4 0,—1,—2,..., tenemos que: 
00 1 
(n) =(_qy+l 
Y9W(x)= (El) nz Era Ft 


paran =1,2,... 


Prueba. De la definición de la función digamma y derivando con respecto a x tenemos 
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que: 


0 + el - +) 


a a | 1 
! = = 
vol ue) 


wW(x)= 


NE» 
el (k + x — 1)2 


derivando nuevamente, con respecto a x, tenemos: 


[oe] 


” == 2 . 
ed La" 


derivando esta segunda derivada, con respecto a x, obtenemos: 


Ses 6 
y!” == : 
(x) 2 +1 


por lo tanto, tenemos que: 


YM) = Wax) = Dn E a 


paran =1,2,... 


Proposición 1.2.41. (Fórmula de recurrencia para la función polygamma): Para x A 
0,—1,-—2,... tenemos que : 


ED! 
xn+l 


Wi (x +1) = Wa(x) + 
Prueba. Usando la fórmula de recurrencia de la función digamma, tenemos que: 
1 
W(x +1) = V(x) + E 
derivando con respecto a x tenemos que: 


wW(x+1)=W(x)- > 
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derivando nuevamente con respecto a x, nos da que: 
YA ” 2 
W(x+1)=Y 0) +3» 


por lo tanto, 


ED! 
Wa (o + 1) = Wi (x) + ar 


paran =1,2,... 


Proposición 1.2.42. (Fórmula de reflexión para la función polygamma): Parax H4 0,—1, 
—2,... tenemos que : 


n 


Wal) + ED) = (07 E 


qañ (cot(77x)). 


Prueba. Usando la fórmula de reflexión de la función digamma, tenemos que: 


W(l—x)-— W(x) = 7 cot(srx) 


derivando con respecto a x, nos da que: 
, , d 
W(l-x)+W(x) = eE (cot(7r x)) 
x 


derivando nuevamente por x, tenemos que: 
d? 
Yw"(1—x)-—w"(x) = 1 =(cot(xx)), 
dx? 
por lo tanto, 


n 


7 (cot(7x)) 


Yl-3+609"*w, (0) = (ED"7 E 
dx 


paran =1,2,... 


Proposición 1.2.43. (Fórmula de multiplicación para la función poligamma): Para x A 
0,—1,-2,...ym > 1 se cumple que 


n—-1 


1 k 
(m) de 
pin+1 PA (++ 5). 


k=0 


ym (yx) = 
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Prueba. De la fórmula de multiplicación de la función digamma tenemos que: 


n—-1 
YW(nx) =1In(n) + - y (s + >) 
k=0 


derivando con respecto a x, tenemos que: 
—1 
IS k 
, PA , UN 
VW (nx) = 72 ) y (++ 
k=0 
derivando nuevamente con respecto a x, nos da que: 
—1 
LS k 
" A " E 
y (nx) = 3 > Y (++ s 
k=0 


por lo tanto, 


LE k 
ym (nx) SS a $ yn) (s E s) 
k=0 


para m = 1,2,... 


1.2.9 Función Gamma en el Plano Complejo 


Hasta el momento hemos estudiado a la función Gamma considerando, como su do- 
minio de definición, el conjunto de los números reales, al que hemos sacado los enteros 
no positivos. Según Davis (1959) la función Gamma es una función de una variable real 
en el sentido en que muchas de sus importantes propiedades funcionan con esta teoría. 
Sin embargo, y ya en los primeros trabajos sobre esta función, se trato de extenderla a los 
números complejos, pero sin mucho resultado. En esta sección ampliaremos el dominio 
de la función Gamma hacia el plano complejo. 

La historia indica que Euler fue pionero en el análisis complejo pero, parece no haber 
considerado la posibilidad de definir el factorial de un número complejo. De acuerdo con 
Davis (ibíd.), la extensión al plano complejo fue iniciado por Gauss y Legendre entregando 
tablas de la función Gamma para valores reales. Entretanto, no fue hasta 1909 que Jahnke, 
Emde y Lósch (1960), basándose en estas tablas, esbozaron un gráfico en tres dimensiones 
de la función Gamma con dominio en el plano complejo: 

El gráfico anterior representa el gráfico de la función z > | (z)|, paraz 4,0,—1,-2.— 
Fasa 

Fue por mediados de 1930 que la función Gamma tuvo aplicaciones en la física teóri- 
ca y en 1950, ahora usando computadores, se publicaron algunas tablas de valores de la 
función Gamma en el plano complejo. 
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Figura 1.5: Gráfico del módulo de la Función Gamma en el plano complejo por Jahnke y 
Emde en 1948 


Si extendemos la integral de Euler a los números complejos, específicamente el lado 
derecho del plano complejo, tenemos: 


n= rtcta: (Re(z) > 0) 


esta función, con términos complejos, coincide con la definición de la función Gamma 
con dominio real. 


Usando el limite de Euler y el producto de Weiertrass podemos extender la definición 
para el plano complejo, excepto para los enteros no-positivos. Para la extensión de la 
función Gamma, en el dominio complejo, probaremos que tiene como polos simples los 
enteros no positivos. A saber: 


Teorema 1.2.44. Para z € C, la función I'(Z) tiene como polos simples a z = —n, 
-1y 

n=0,1,2,...y con residuo Res(I”, —n) = ( , ] 
n! 


Prueba. Sea P(z) = fi 1?7lertdt = fi1?7lertdr + [0 127 Le7t dt, motemos que si 
[,0) 


. . Ñ (=D 
expandimos la función e”? como serie de Taylor, centrada en O, es e”? = EN y 
n! 


n=0 
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tenemos que: 


T(z) =f y EVE - le enidt 
Jo n! 1 


00 1 jn+z-1 00 
=D d+] de 
n=0 0 ni 1 
1 n pn+z 1 [0/) 
y Ed | | +/ ch 
“ol at Ln+z]o 1 


SI 
o 


[0,0] 
La integral v(z) = / e”* 1271 df define una función analítica de la variable z y, por 


1 
tanto, aplicando el teorema de residuos, tenemos que: 


So 


Res(I”,—n) = Him ¡E = (=npyI(z) = ,n E No. 

Concluimos que /'(z) es una función definida en los números complejos menos los enteros 
no positivos, con polos simples en esos puntos y, dado este hecho, la función Gamma es 
analítica en el plano complejo menos los enteros no positivos. 


Para ver más propiedades de a función Gamma en el plano complejo, recomendamos 
ver los ejercicios. 


1.2.10 La Función Gamma Incompleta 


Para definir la función Gamma usamos una integral evaluada en el intervalo ]0, oo| 
podemos, escogiendo x, con O < x < oo, dividirla en dos integrales, evaluadas en ]0, x[ 
y ]x, oo! y, con ello, definir las funciones y (a, x) y (a, x) por: 


y (a, x) =| le tdt; H(a)>0 
0 


[0,0] 
rasw=/ ttetadt 
Xx 
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Estas dos últimas funciones están relacionadas con la función Gamma por: 
(a) = y(a,x) + I'(a, x). 


Para el caso de que el parámetro a es un número entero positivo, tenemos: 


Proposición 1.2.45. Para n e Z*, se cumple: 
n—1 xk 
= 1)! ZW ga — 
y(n,x)= (n—1)! ( e 2 =) . 
=0 


Prueba. Sea y(n,x) = ae "=Le=tdt, si integramos por partes con u = ("71 y du = 
e”* dt tenemos que: 


Xx 
y(n,x) =/ olas 
0 
Xx 
= A A + (n — Df tetdi 
0 
Xx 
=-— xml 4 (n— Df dE, 
0 


Recursivamente: 
y(m,x) = xP 18% — (mn — 1x7... (mM Dixe * + (n- Mee” +1) 
por lo que que: 


n—1 yk 
y(n,x) =(n—1)! ( a 5) 
k=0 


como señalado. 


Proposición 1.2.46. Para n e Z* tenemos: 
n—1 yk 
ps qe, —X Ea 
T(n,x) = (n— le * Y” ak 
k=0 


Prueba. Integrando por partes, al igual que en el caso anterior, tenemos para: u = 1"71 y 
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dv = e dt, entonces 
Se [0,2] 
T(n,x) = 20 e [7 +(n— Df ro? tdr 
me x 
=x" Y 4 (n— Df Sed 
x 
Recursivamente llegamos a: 
T(n,x) =x"" le * 4 (n= Dx? 4... +(n—1)le* 


por lo que: 


n—1 xk 
T(n,x) = (n— 1)le 2 mí 
=0 


como afirmado. 


1.3 Aplicaciones de la Función Gamma 


En sus diferentes representaciones (límite, producto infinito e integrales), la función 
Gamma puede ser, y a veces es, una importante herramienta a la hora de hacer diferentes 
cálculos. Ya vimos cómo se gestó la función Gamma y señalamos sus principales propie- 
dades, ahora presentaremos algunas aplicaciones. 


1.3.1 Para el Cálculo de Integrales 


Algunas integrales lucen bastantes intimidantes, pero se vuelven más amigables si 
aplicamos de manera asertiva la función Gamma (o la función beta). Daremos algúnos 
ejemplos ilustrativos y muy detallados que recopilamos de Attar (2006) y Farrell y Ross 
(1971). 


[o,0) 
Ejemplo 1.3.1. Calculamos la integral e da 


0 
00 


Usando la proposición 1.2.21 sabemos que: T'(x) = a* / 412741 dí Evaluando para 
0 
a =2y x = 4, tenemos: 


[0,0] 
T (4) = 0 ode 
0 


00 T(4 31.03 
3, =2x 
== ===: 24 
/ A it rt 


entonces 
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' d (=D! 
Ejemplo 1.3.2. Probar que: / x” (Ino) dx = mA paran,m € Ñ. 
En este caso usamos el cambio de variable u = — In(x), entonces reemplazando en la 


integral, tenemos que: 
0 00 
a E E 
00 0 
[0,2] 
=c0" f y” ¿M+DA y 
0 
y usando nuevamente la proposición 1.2.21 parax =n+1ya =m + 1, tenemos: 
[,0) [,0) 
TMm+1)= (m+ pan] Me MD (m+ pa y” ¿M+Dua gy. 
0 0 


Por tanto, y ya que f, u"e(7+D% gy = (-1)" L x”(In(x)) dx, entonces (n+ 1) = 
(m+ 13 + 210" L x" (In(x))"dx y ya que '(n + 1) = n! obtenemos: 


ED! 


1 
/ x” (nto) dx = Da 


como pedido. 


Ejemplo 1.3.3. Calcular [<= 


<= 


Para este cálculo usamos, nuevamente, el cambio de variable u = —In(x) y reempla- 
zando en la integral tenemos que: 


0_ 4 q 00 1 
/ il =/ dy =R (5) = VI, 
00 u2 0 


Ls" 


2 f A 1 
Ejemplo 1.3.4. Calcular ÍA ( 
0 


esto es: 


4 
sinó (x) y 0) eii 


Haciendo la diferencia del paréntesis, tenemos que: 


A A 
/ (a 505) coscodx = [| ( PTS ) cos(x)dx. 
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Usando el cambio de variable u = 1 — sin(x) y reemplazando en la integral tenemos: 


9 4 y La 3 
[(H») —du= | ut(l1—u) +du. 


Usando la definición 1.2.4 de la función beta, B(x, y) = e 21 (1 — £)* dt; con los 
valores x = 5 y y = 2, y la proposición 1.2.28 obtenemos que: 


5 1 1 
pe. 2) = PF) =i uz(1—u) du. 
0 


44. TÉ+D 

Ahora usando la proposición 1.2.3 y la proposición 1.2.7 tenemos rG) = 114) y 
TU 

TG) = Lr(3) = Sn entonces 


e 3 PIQOIQ [! ep 
a — 4 4 4 4 
1 ui (1 —u) *du = T. a = E" 


por lo tanto: 


Tí 1 


Epa Y Ir 
/ (e +05) ts 27 


1.3.2 Volumen y Area de una Bola y una Esfera en Dimensión » 


Una bola n-dimensional de radio R, consiste en el lugar geométrico de los puntos, 
del espacio euclidiano de dimensión n, cuya distancia al origen es menor o igual a R. 
Recordemos que las coordenadas de un punto en el espacio euclidiano de dimensión n, es 


la n-upla (x1,x2,..., xn) € R”. En forma de ecuación, la bola n-dimensional de radio R 
corresponde al conjunto de puntos (x],xX2,..., Xp) € R”, tales que: 
ARA PS RO. (1.28) 


Para calcular su volumen, integramos infinitesimalmente el elemento de volumen dV = 
dxi ...dx, sobre la región de n-dimensiones , indicado en (1.28). Entonces el volumen 
de la bola, V,, (R), será: 


Va (R) = p dxidx2...dxpn = Cp R” (1.29) 
A ES 


donde C,, será una constante a determinar. 
Ahora, si consideramos la bola n-dimensional como un conjunto de capas (cáscaras) 
concéntricas, entonces tenemos que el elemento de volumen d V,, de la capa de radio R y 
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cuyo espesor es dR está relacionado con el área A, (R) de la esfera (n — 1)-dimensional 
por d V,(R) = AndR, entonces 


o a a, ge (1.30) 
AN 7 


Veamos una manera, usando a la función Gamma, de calcular el valor de la constante 
C, . Para ello, consideremos la función: 


2 


24 y2 Ll 
Pda) 0 Qi+x5+h... ex) — e R 


si integramos esta función sobre las n-dimensiones, obtenemos: 


27 2 ES ES 2 2 2 
/ e dv, Sl 0] e Hit Tt +A dadas. din 
0 00 —00 


de la ecuación (1.30), tenemos que: 


oo 2 28 2 
/ e ALGR =p e nC, RAY aR 
0 0 
[0,0] [0,0] 
E Pe / a das. de 
00 00 


por tanto, obtenemos: 


ee 2 22 2 2 2 99 2 
"Cr | ESREAR El e *Tdx;y / ida f et”, (1.31) 
0 =00 =00 —00 


Ahora reemplazando en el lado derecho y lado izquierdo de la ecuación (1.31) por la 
proposición 1.2.6 y la observación (reemplazando x por 3) respectivamente, obtenemos 
que: 
n 
nCP(5) = 53. 
2 

entonces 7 

m2 
 TG+D 


Por lo tanto de la ecuación (1.29), tenemos que: 


Ca 


Asimismo, sabemos de la ecuación (1.30), que el área de una esfera (n — 1) —dimensional, 
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denotada por A, (R), es: 
d 
(Va (R)) = An(R), 
entonces 
d 123R” nR" 177 NR" 2 2712 
dRÍTG+1 TrG+1) 51 (3) TG) 


n—1 
R , 


por lo tanto: 
2717 
RE: 


AR) = Fa 
2 


Un segundo método, para calcular el volumen de una bola n-dimensional, es hacer el 
cálculo de forma directa por medio de un cambio de coordenadas -en el espacio euclidiano 
n-dimensional- de coordenadas rectangulares a coordenadas esféricas. 

Así, para calcular la integral de la ecuación (1.29): usamos el cambio de coordenadas 
esféricas: 


x¡ = Reos(p1) 

x2 = Rsin(p1) cos(p2) 

x3 = Rsin(Q1) sin(p2) cos(p3) 

X4 = Rsin(Q1) sin(p2) sin(p3) cos(pa) 


Xn-1 = Rsin(p1) sin(p>) - -- sin($a-—2) cos(pn—1) 
Xn = Rsin(p1) sin(p2) --- sin(Py-1). 


Es fácil comprobar que xi +x3 +... +32 = R?, 
Al hacer este cambio de coordenadas, la integral en (1.29), queda de la forma: 


vr =f.-/ dxidx>...dXpy 
xxx <R2 


(1.32) 
= ff s1ardórdos---dón-, 


donde J es el correspondiente Jacobiano al usar el cambio de coordenadas. 


Observacáo. Las nuevas coordenadas corresponden al radio R y a las n — 1 coordena- 


das angulares, H1,Q2,...,Pn—1, donde los ángulos P1,42,...,Pn-2 € [0,1 [| y $n-1 € 
[0, 27 | radianes. 
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Ahora, el Jacobiano de la transformación es: 


0X; 
Ide 5) con =1,...,n—=1;i¡=1l,...,1n; 
AR, pj) 


no es difícil verificar que: 


J = RI sin"2(41) sin""3 (47) --sin(ón_»), 


entonces reemplazando en (1.32), tenemos: 


Y, (R) = Jas) sinó -ddRdgrdo2 den 
21 


R T 
= | RAR / del ds 
0 0 0 


que es lo mismo que tener: 


Va (R) =f. fu sin""(41) ---sin(da-2JdRdg do» -- db —1 


R” E 2 
AA (2 sin sn) + (4 dtr). 


Ahora, usando las proposiciones 1.2.25 y 1.2.28 y reemplazando en la ecuación (1.33), 


a) 
e O O 


(1.33) 


n r() ra) TO 
R” 1 1 1 
LO 
n—veces 
R"x2 
31 (3) 
por lo tanto: 
AR = 7 
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Ejemplo 1.3.5. Calculemos el volumen y el área de una esfera de dimensión n = 3. 


3 R? 73 R3 E H3R3 4 , 
VAR) == ET a e ds 
Pl 1) E) 2" 316) ¿7? 
y 
217R2 213R?  473R? 
A2(R) = e = = 4xR?. 


A 


1.3.3 Distribuciones en Probabilidad y Estadística 


En esta parte veremos como la función Gamma puede usarse para definir distintas fun- 
ciones de distribución acumulada, en estadística y probabilidad. Para ello recordaremos 
algunas definiciones de teoría de la medida y de probabilidades. 


Sea X un conjunto no vacío y $2 € P(X) una colección de subconjuntos de X. 


Definición 1.3.1. Diremos que $2 es una álgebra en P(X) si ocurren las siguientes pro- 
piedades 


(1) Si A € $ entonces (X 1 4) € 2 y 
(11) (aditividad) Si 41, A2 € $2 entonces 41 U 4» € S2 y 
(iii) X € $2. 


Cuando la colección £2 satisface solamente las condiciones (1) y (11) se dice que es un 
anillo en P(X). 


Cuando la colección $2 cumple las condiciones (1), (111) y la siguiente condición de 
o-aditividad: Si Aj € S2 Vi € N entonces US ¿A;¡ € $2; decimos que £2 es una 0 -álgebra 
en P(X). 


Ejemplo 1.3.6. Si Ocurre que el par (X, T) es un espacio topológico entonces la menor 
o-álgebra que contiene a los abiertos de la topología Tse la llama la o -álgebra de Borel 
de X. Para obtener esta o -álgebra basta intersectar a todas las o-álgebras que contienen a 
los abiertos de 7. Por ejemplo P(X) es una de ellas. 


Al par (X, £2) se le llama espacio medible si X es un conjunto no vacio y £2 es una 
o-álgebra en P(X). 


Vamos ahora a definir el concepto de medida sobre un espacio medible. 


Definición 1.3.2. Sea $2 una 0-álgebra de subconjuntos de X. Una función u : $2 —> 
[0, 00] se dice una medida sobre $2 si: 
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(a) (9) =0; 


(b) Para toda familia enumerable de conjuntos dos a dos disjuntos 4;;i € N se cumple 
que 1(U;A¡) = 2 u(A;). 


Un espacio de medida es una terna (X, $2, 1) donde X es un conjunto no vacio; £2 es 
una 0 -álgebra de subconjuntos de X y yu es una medida sobre $2. 


Ejemplo 1.3.7. En R” consideramos la o-álgebra de los bolerianos B. En este caso existe 
una única medida 4 : B > [0, oo[ tal que para cada rectángulo A = [a1,b1] x [az, b2] x 
--*+X [4 , b, ] se cumple que A(4) = (b, —a1) x(b3—a2)x:-+Xx (by —an). Esta medida se 
llama la medida de Lebesgue en R”.En lo que sigue consideraremos el espacio de medida 
(R”,B,1.). 


El espacio de medida (X, $2, 1) se dice un espacio de probabilidad si ocurre que 
p(X) =1. 


Definición 1.3.3. Consideremos dos espacios medibles (X, $2) e (Y, B). Una función f : 
X > Y se dice medible si ocurre que VB € B se cumple que f71(B) e 32. 


Ejemplo 1.3.8. Una variable aleatoria a valores reales del espacio medible (X, $2) es 
una función medible £ : X > R, en donde hemos considerado la d-álgebra de Borel en 
R. Si ocurre que la imagen ¿(X) C R es un subconjunto finito o infinito enumerable se 
dice que la variable aleatoria E es discreta. En el caso que existan a,b € R,a < b tales 
que [a, b] € ¿(X) diremos que £ es una variable aleatoria continua. 


Consideremos el espacio medible (X, $2) y £ : X —> R una variable aleatoria continua. 


Para valores a < benR oa = —000b = oo podemos formar los subconjuntos medibles 
de X: 
[la <E<b]= fx € X/a < ¿(x) < Dd), 
O 
[a <E <b] =(x € X/a < ¿(x) < by, 
O 
[a <E <b]=(x € X/a < ¿(x) < by, 
O 


[a <E<b]= ix € X/a < 5(x) < bj. 


Consideremos ahora una medida de probabilidad P : $2 > [0, 1]. Tenemos los valores 
reales P([a < € < b]) definidos para cada variable aleatoria £ y cada par de números 
reales a < b. 
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Definición 1.3.4. Una función integrable, no negativa f : [0,1] —> KR, se dirá una densi- 
dad de la variable aleatoria £, asociada a la probabilidad P, si ocurre que: 


b 
Pasé<bp= f food. 


Esto es el valor P((x € X/a < E(x) < bj) coincide con la integral de la función f(x) 
en el intervalo [a, b]. 


En particular, se cumple que: f FG)dx= 1. 


Observamos que el intervalo de definición de la función de densidad se puede extender 
a todo R, considerando f(x) = O para—o0<x<001<x+<o00. 


Ahora que hemos definido la función de densidad de una variable aleatoria podemos 
definir la media, la varianza y la función de distribución asociadas a dicha variable aleato- 
ria. A saber: 


Definición 1.3.5.  (i) la media o esperanza de la variable aleatoria É es el valor yu defi- 
nido por: 


n=E)=/f x > foddx: 


(11) La varianza de la variable aleatoria É es el valor: 
[0,2] 
var) = | 0? gx: 
—00 


(111) La función de distribución acumulada o simplemente función de distribución de la 
variable aleatoria É es la función real F : R > R definida por: 


Fo) =PE<xD=f_ 10 


Ahora, en general, podemos suponer que cualquier función no negativa f : R > R 
que satisface: : f(Go)dx = 1, es una función de densidad de alguna variable aleatoria 
asociada a alguna medida de probabilidad en algún espacio de eventos X. En este caso la 
función F(x) = ES F(s)ds; es la función de distribución definida por la densidad f(x). 
Observamos que si la función de densidad es continua, entonces la función de distribución 
es derivable y 


FU(x)= f(x). (1.34) 


No es dificil verificar (ver Saavedra Gallardo (2012)), que una función F(x) de distri- 
bución (asociada a una variable aleatoria £) satisface las propiedades: 


a) 0 < F(x) < 1 para cualquier x e R; 
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b) P(a < E < b)] = F(b) — Fla), para a < b; 

c) Sia < b entonces F(a) < F(b) esto es, F es no decreciente; 
d) limy>0 F(x) = 1 y limx>-x F(x) = 0; 

e) lim, + F(x) = F(t), es decir, F es continua por la derecha; 


D P(É =1t) = F(t) — F(t7), donde F(t7) = limx->1- F(x). Es decir, P(É = t) es 
el tamaño del “salto” de F en £. 


Ahora podemos, al revés, dada una función F : |R > R que satisface las propiedades 
[a)] — [F)] anteriores, definir una función de densidad f(x) por la relación (1.34). 

Así, dada una función de densidad tenemos una de distribución y dada una de distri- 
bución tenemos una de densidad. 

A lo largo del tiempo se han generado diversas funciones de distribución más conoci- 
das. Veremos algunas. 


Distribución Gamma 


Esta distribución de probabilidad continua se utiliza para modelar el comportamiento 
de variables aleatorias con asimetría y/o algúnos experimentos en que está involucrado el 
tiempo. 


Definición 1.3.6. Sea X una variable aleatoria continua para la que suponemos que su 
función de densidad es: 


xo 
Fx,0,PB)=3 BT (a) 
0 


para x>0,0,B>0 


para otro caso. 


De esta distribución, podemos mencionar tres casos particulares: 


a) Distribución de Poisson: Distribución de probabilidad discreta que expresa, a partir de 
una frecuencia media, x, la probabilidad que ocurra un determinado número de eventos 
a € X, durante un intervalo de tiempo dado. 


Definición 1.3.7. Sea X una variable aleatoria que representa el número de eventos 
aleatorios independientes que ocurren a una rapidez constante sobre el tiempo o el 
espacio. Se dice, entonces que la variable aleatoria X tiene una distribución de Poisson 
con función densidad de probabilidad: 


xelox 


fo,a,1)=3 Tía) 
0 para otro caso. 


para x > 0,0 € Z+* 
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b) 


1. Función Gamma 


Si hacemos «4 — 1 = A tenemos la función densidad de probabilidad que aparece 
generalmente en la literatura: 


er 


CREO ESTO para x>0,1 € Z¿ 
o para otro caso. 


Xx 


Distribución exponencial: La distribución exponencial es utilizada para determinar la 
probabilidad de que en cierto tiempo suceda un determinado evento. 


Definición 1.3.8. Una variable aleatoria X tiene una distribución exponencial si su 
función de densidad está dada por: 


1 2e7%% para x>0,1>0 
1, = 
SS) 


para otro caso. 


La distribuciones de Poisson y exponencial, aparte de ser casos especiales de la distri- 
bución Gamma, guardan otra relación relevante y simple. 


Supongamos que X sea una variable aleatoria, con distribución de Poisson con pará- 
metro A, que describe el número de ocurrencias de un fenómeno por unidad de tiempo, 
entonces la variable aleatoria que describe el tiempo entre ocurrencias consecutivas 
tiene distribución exponencial con parámetro A y en tal caso este parámetro A es la 
tasa media por unidad de tiempo, y $ = z es el tiempo entre estas ocurrencias. Es 
decir, haciendo uso de la distribución de Poisson, podemos asumir que la probabilidad 
de que no ocurra algún evento, en el periodo hasta el tiempo £, está dada por: 


0,—At 
FAt,0,1) = P(X =0) = AZ a 


Podemos utilizar lo anterior y hacer que X' sea el tiempo para el primer evento Poisson. 
La probabilidad de que la duración del tiempo hasta el primer evento exceda x es la 
misma que la probabilidad de que no ocurra algún evento de Poisson en x. Esto último 
está dado por e74*, es decir: 


P(X>x)=e?%,. 
Asi la función de distribución acumulada para X está dada por: 
P(0<X<x)=1-e*, 


Ahora para reconocer la presencia de la distribución exponencial, podemos derivar la 
función de distribución acumulada anterior para obtener la función densidad: 


M0) =1e 94 
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h ad A 1 
que es la función de distribución exponencial con 4 = >. 


c) Distribución ji cuadrada (o xy-cuadrada o de Pearson): Otro caso importante que se 
desprende de la distribución Gamma, se obtiene al hacer « = 5 y f = 2, con v € 2, 
Este resultado se llama distribución ji cuadrada. El parámetro v, es llamado grado de 
libertad. 


Esta distribución juega un papel importante en inferencia estadística y además tiene 
aplicaciones en meteorología. 


Definición 1.3.9. Sea la variable aleatoria X continua. Diremos que X' tiene una dis- 
tribución ji cuadrada, con v grados de libertad, si su función de densidad está dada 
por: 


1 v Xx 
v —_ xx 37 lez para x>0,veZ+ 
52) =3221'(5) 


2 
0 para otro caso. 


Fa, 


Distribución Beta 


Otra distribución continua de uso frecuente es la distribución beta, que suele utilizarse 
para modelar fracciones tales como la proporción de impurezas en un producto químico 
O para modelar eventos que se definen por valores mínimos y máximos. La escala de la 
distribución beta suele modificarse para modelar el tiempo hasta la culminación de una 
tarea. También se utiliza en la estadística bayesiana. 


Definición 1.3.10. Sea X una variable aleatoria . Decimos que tiene una distribución beta 
si su función de densidad esta dada por: 


Bl. QxP1(1-x31 para O<x<l1,p,q eR* 


Fx, pq) = 0 para otro caso 
T 
con B(p,q) = AE 


Distribución de Weibull 


Como la distribución Gamma y la exponencial, la distribución Weibull también se apli- 
ca a problemas de confiabilidad y de prueba de vida como los tiempos de falla o duración 
de vida de un componente, medido en algún tiempo especifico hasta que falla. 


Definición 1.3.11. Decimos que la variable aleatoria X tiene una distribución de Weibull, 
con parámetros a y $ si su función densidad está dada por: 


B—-1,—axb 
apx e ara x>0 
fra, p) = IP p 


0 para otro caso. 
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A simple vista esta distribución no tiene ninguna relación con nuestra función Gamma, 
pero esta aparece al obtener la esperanza y varianza de la distribución Weibull. 


Teorema 1.3.1. La esperanza y varianza de la distribución Weibull son: 


E[X] =0 ET (5 $ 1) y 


[X] =078 r(7+1) [+] 
var =0 = - 5 . 
13 B 
Prueba. Recordemos que la esperanza de una distribución de probabilidad es: 
Epx=f xfa) dx 


siendo f(x) la función densidad. Dado que en nuestro caso la función densidad es válida 


para x > 0, tendremos: 


E[X] =/ xapxé 120 gy 


usando el cambio de variables «x4 = 1 y reemplazando, tenemos 


A) 
LOGO 


af tá 1 2, MN (AA 
=q $ 1B e*dt=0 8 1Be dt 
0 0 


y de la definición 1.2.1, tenemos que: 
1 1 
E[X] =0 ir (7+1) 
$ 
La varianza de una distribución de probabilidad es: 
var [X] = E[X”?] - (E [X)? 
como ya tenemos E [X'], calculamos 


[,0) 
E[x?] sl Papxi g 
0 
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Usando el mismo cambio de variable 1 = a.xÉ, tenemos 


OM 


[,0) 
=2 Zit=l = 
=0 il Ud 
0 


y usando el definición 1.2.1, nos da 


entonces tenemos que: 
var [X] aa ( ale 1) _ [enár ( n 9 
ra r (5 + 1) = [r (5 + y 


Distribución de Student 


La distribución de Student surge del problema de estimar la media de una población 
normalmente distribuida cuando el tamaño de la muestra es pequeño. 


Definición 1.3.12. Una variable aleatoria se distribuye según el modelo de probabilidad 
de Student con n grados de libertad si su función de densidad es: 


rl) ¿8 (+1) 
2 
t,n) = 1 = <t=<oo. 
f(t,n) To Al + =) 00 00 


+00 
Vamos ahora a probar que ft,n)=1. 
—00 


+00 p(utl p -H(n+1) 
I =f o) ( + ) dt 
-oo T(5)ynx n 


ra 00 pp 2 
= === 2 1+ 2) dt. 
PG)ynr do ( n 


Tenemos: 
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¿ 2 
Si hacemos u = 5, tenemos que: 


EE) e 1 Li ed 
I a | 1+4)20+D. /n- 4 2du 
P5H)vnz 0 ( ) yn 2 
1 
OVA (149400 
00 uy! 
y dado que B(x, y) =/ A = Ly x = 5, por tanto: 
rea 1 
he G -B(5, 
A a 


IG TOTO) 
GV T(G+3) 
=1 


esto es: / = 1, como queríamos probar. 


Distribución de Snedecor 


Esta distribución es muy utilizada a la hora de detectar la existencia o inexistencia de 
diferencias significativas entre muestras diferentes. Especialmente, en aquellos casos en 
los que se quiere investigar la relevancia de un factor en el desarrollo y naturaleza de una 
característica. 


Definición 1.3.13. Sean las variables aleatorias independientes X e Y con n; y n, grados 
de libertad respectivamente, tenemos la variable aleatoria 

n2X 

niy 


F = 


cuyo campo de variación es F' > 0, tiene por función de densidad de probabilidad a: 


n ni-2 

PEE AN E HD 
M(F,n1,n2) = n1 a n2 ( -) ni+noy * 
TEST) wa (A 


F 
(1 + mE) 
+o0 


Vamos ahora a probar que J(F,n1,n2) = 1 y como el campo de variación es 
—00 
+o0 


F > 0, hay que probar que f(Ejniyn2)=1 
0 
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Tenemos: 


pe TR) ON F 
¿al r(2)r(3)Ua si 


PA je 
= n1 n2 ( ) / n3-n2 dF 
TAG) 0 (1+ me (13) 
n2 


] n1 F 
Si hacemos u = ——, tenemos que: 
n2 


m1 


A (EN O E jo. yaa 
13 m1 n2 : / n3+n2 du 
MATE) m1 ni Jo w2 


nu) ma a 1 ¿ 
= n n n U 
MES A A 


uy! 


y dado que B(x, y) = / (+ uysty 


du, usamos y = % y x = %, por tanto: 


r(34) n2 M1 
rara 6) 
ra TT) 
OE rt) 
=1 


IT = 


esto es: / = 1, como queríamos probar. 


1.3.4 La Función Gamma en Soluciones de Ecuaciones Diferenciales 


La función Gamma también aparece en la búsqueda de soluciones de diferentes ecua- 
ciones diferenciales. 
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Ecuación Diferencial de Bessel 


La ecuación diferencial x?y"(x) + xy'() + (4? — a?) y (x) = 0, conocida como de 
Bessel, cuyas soluciones son de la forma: 
y(x) =aJa(x) +bJ_a(x) para QU É Zi v 
y(x) =a1Ja(x) +b1Ya(x) con Y, (x) = lim Yo (x) VWn e Z* 
a—>n 


2k+p 
MES) Jo (x) cos(azr) — J-a(x) 
J, y Y, = 
Ea Mo SENOS +0+ Dn” a(x) sin(o77) 
Para os detalles sobre la ecuación de Bessel recomendamos el texto Gui nez, 
Labarca y Martínez (2004). 


Ecuación Hipergeométrica de Gauss 


Para a, b y c constantes reales, se llama ecuación hipergeométrica de Gauss a la ecua- 
ción x(1 —x)y"(x) + [e —- (a+ b+1)x] y'(<) —aby(x) = 0, que tiene como solución 
a las funciones: 


y(«) =0aF(a,b,c,x) + Blx“F(a—c+1,b=c+1,2—c,x) 


con 
T(a+)T(W+IT() y 
T(a)T (9) (c +Kk) 


Fabeo=t+ Y 


con q, $ constantes, c ¿Z,c %4a+b-+ 1 y la cual está definida para 0 < |x| < 1 
Para mayores detalles sobre la ecuación de hipergeométrica de Gauss recomendamos 
el texto Gui nez, Labarca y Martínez (ibid.). 


1.3.5 Aplicaciones en el Campo de la Física 


La función Gamma aparece ocasionalmente en problemas en física por ejemplo, en 
la racionalización de las funciones de onda de Coulomb y el cálculo de probabilidades 
en la mecánica estadística, entre otros casos. Sus aplicaciones, en la física, son menos 
cuando comparadas con otras funciones, pero su importancia está en el desarrollo de otras 
funciones que tienen una implicancia directa en la física y, por ello, su inclusión en este 
trabajo. 


Estimación de Potencial Eólico usando la Distribución Weibull 


Una de las formas más comunes de analizar los datos de velocidad del viento Oyrazo 
y Mancilla (2006), y de elaborar estudios sobre energías renovables basados en energía 
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eólica es hacerlo mediante una función densidad de probabilidad, en este caso la distribu- 
ción de Weibull. Esta describe con bastante confiabilidad la distribución de la velocidad 
del viento para un intervalo de un mes o hasta para un año. 


La función de densidad de probabilidad de Weibull para la velocidad del viento (v) 
está dada por la siguiente expresión: 


ala 


e: Sd 
qe pe) para k,v>0nc> 1 
0 para otro caso. 


sf 


Esta densidad depende de dos parámetros: c que se denomina parámetro de escala y k 
que es el factor de la forma. Existen varios métodos para calcular estos parámetros, y 
estos claramente van a variar de mes a mes dada la estación del año en que estemos y la 
localidad geográfica. Uno de estos métodos es cuando conocemos la velocidad media (v) 
y la desviación estándar (o) de la velocidad, entonces una buena aproximación del valor 


k es: DEBO 
k=(5) > 


El valor c, podemos calcularlo de la siguiente manera, recordemos que la velocidad media 
se puede obtener mediante la integral v = Le vf(v)dv, si f(v) es nuestra función de 
densidad de probabilidad Weibull, tenemos que: 


00 vk V k-1 k 
v7= / pes (5) Y dy 
0 Cc C 
A j : k 
si hacemos el cambio de variable x = (2) , Obtenemos que 
e A 
v= e] xke *dx 
0 
y ahora haciendo y =1+ : tenemos que 
[0,0] 
v= e] Soda 
0 


lo que es justamente nuestra función Gamma, por lo tanto: 


También, dada esta función de densidad de probabilidad, podemos calcular la potencia 
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promedio del viento, es decir: 


SA 1 ES 

Py = ¿a v? f(v)dv 
2 0 

siendo p la densidad del aire. Nuevamente reemplazando f(v) por la función de densidad 

de probabilidad Weibull, y haciendo los mismos reemplazos obtenemos que: 


Py = pAcT ( + 2) 


y dado que c = , tenemos que: 


P(1+2¿) 
7 _PAWT (145) 
[rar pro 


Función de Luminosidad de Galaxias 


La función de luminosidad de galaxias es una aplicación de interés en astrofísica, ya 
que revela la distribución de los distintos tipos de galaxias según su luminosidad. Este 
parámetro depende de cómo se han formado y cómo han evolucionado durante su vida y, 
por ello, proporciona información sobre ambos procesos. La luz que recibimos desde las 
galaxias es el resultado también de sus interacciones con otras, por ello su importancia 
Coendra (2008), Carroll y Ostlie (2007). En la práctica la función de luminosidad (LF) 
nos da la densidad numérica de galaxias en función de la magnitud, la cual puede ser dife- 
rencial, número de galaxias por unidad de magnitud (o luminosidad) o en forma integrada, 
número de galaxias con magnitud superior a una magnitud dada. Sea v(M, x, y, z) el nú- 
mero de galaxias que yacen en el volumen dV en (x, y, z) con magnitudes entre M y 
M + dM. Si suponemos que las magnitudes de las galaxias no se correlacionan con la 
posición, resulta: 


v(M,x, y, 2) dMdV =9(M)ID(x, y 2d MdV 


donde as 
/ Pp(MdM =1 


y donde p(M) da la fracción de galaxias por unidad de magnitud con magnitudes abso- 
lutas comprendidas entre M y M + dM y recibe el nombre de función de luminosidad. 
D(x, y, z) da el número de galaxias (de todas las magnitudes) por unidad de volumen en 
(x, y, z) y se llama función densidad. A través de las observaciones se encuentra que la 
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LF puede ser ajustada por la función de Schechter (1976): 


100,4M-M*) 


$(M)JdAM =0,41n(10)p* 1079 4M-MD(0-D 4 


o escrita en términos de la luminosidad 


LY” _1d 

Los parámetros en esta última ecuación son la normalización p*, la pendiente logarítmica 
a de la parte débil y la magnitud (o luminosidad L*) característica M*, la cual marca la 
transición entre el comportamiento exponencial de la parte brillante y la ley de potencias 
con exponente a de la parte débil de la LF. Cuando se determina la LF se debe especificar 
de qué tipo de objeto se trata y en qué ambiente se encuentra. 

Ahora el número de galaxias total por volumen, es decir, la densidad de galaxias, usan- 
do la función de Schechter será: 


Pena = UN p* (2) E 


si hacemos el cambio de variable £ = L+ tendremos 


/ P(DdL =] oe dt =9*/ ed, 
0 0 0 


y con fo 1%e7*dt = T (a + 1), por lo tanto, 


Poma =$*T (a +1). 
0 


También podemos obtener la densidad de luminosidad total, es decir: 
[,0) 
pL= [suya 
0 


y nuevamente usando la función de Schechter y el cambio de variable £ = Le tendremos 
que 
pL=9*L*T(a +2). 


Recordemos que la función Gamma se extiende para x < O (donde diverge o es nega- 
tiva o es positiva). Contando con que « toma valores típicos entre —0, 5 y —2.0, entonces 
la densidad de galaxias es finita (y positiva, como debe ser) solo si « > —1 (ya que en 
este caso O < 4 +1 <o0, 5) y la densidad de luminosidad de todas las galaxias que entran 
en una LF tiene un valor válido si « > —2 ( ya que en este caso0 <a +2<1,5). 

Ahora si queremos el número de galaxias en un intervalo de luminosidad determinado 
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se necesita la función Gamma incompleta para un L > £j: 


sa = tw (a ep 2) E 
da Es 


Sistemas 
Dinámicos 


2.1 Introducción 


En este capítulo presentamos resultados originales, respecto de la dinámica de la fun- 
ción Gamma. 

En la sección 2.2 enunciamos definiciones y diversos resultados necesarios para las 
secciones siguientes. 

En la sección 2.3 probamos diversos resultados sobre la dinámica de la función Gam- 
ma. A saber, estimamos el punto de mínimo en x > 0; establecemos y estimamos puntos 
fijos en x > 0; estimamos soluciones de la ecuación 7"(x) = —2 en los intervalos ]—3, —2[ 
y ] — 5, —4[] y, además, probamos que la derivada de J” en estos puntos es, en valor abso- 
luto, mayor que 8. También, calculamos la extensión del conjunto estable del punto fijo 
atractor x = 1. 

En la sección 2.4 damos la definición de la entropía topológica; algunas de sus pro- 
piedades, un algoritmo para calcularla; la extensión de la definición a nuestro contexto; 
calculamos las derivadas T”'(xx), T”'(yx) donde xz e yz son las soluciones de Ty (x) = —2 
para x €] — 2k — 1,—2k[ y k > 3; usamos este hecho para construir una conjugación to- 
pológica entre el shift de 2n simbolos y el conjunto maximal invariante de /” en la unión 
de intervalos U+=1 (1-27 — 1, —2][); finalmente usamos este resultado para probar que 
la entropía de /' es infinitamente grande. Como consecuencias sobre las estimaciones de 
la derivada concluimos que el conjunto maximal invariante en Us=1 (-2j-1,-2j[ es 
hiperbólico Vn > 1. Asimismo, dejamos planteada una pregunta para futuras investiga- 
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ciones en el tema. 


2.2 Previos 


Definición 2.2.1. Decimos que una función f : R > R posee un máximo (respectiva- 
mente mínimo) local en el punto p € R si existe $ > O tal que q € R,con|p — ql < 6, 
implica f(q) < f(p) (respectivamente f(q) > F(p)). 


Teorema 2.2.1. Sea f : [a,b] > R una función derivable en xy €la, b|. 


a) Supongamos que f'(xp) > 0. Entonces existe h > 0 tal que si Xy < y < xp + h se 


tiene que f(x0) < F(y). 


b) Supongamos que f'(xo) < 0. Entonces existe $ > 0 tal que si x, y € [xo — 8, xo + óL, 
con x < y, se tiene que f(x) > f (y). 


Prueba. Por definición de derivada, tenemos que f'(x0) = En oa mA 
5 — X0 
a) Dado que f'(xo) > 0, tenemos que f'(x0) = ee FO F0) > 0, luego existe 


—Xo 
un h > Otal que xo < y < xo +h implica que L0-£00) > 0 y dado el hecho de que 
y — Xo > 0 se tiene que f(y) > £(x<0). 
Análogamente, se prueba que si f'(xo) > 0, entonces existe h > 0 tal que xo —h < 
y < xq implica que y — xo < 0, entonces f(y) < f(x0). 
b) Dado que f'(xo) < 0, tenemos que f'(x0) = lim FM F0) > 0, luego existe 
yY>X0. Y —X0 
un ó > Otal que xo < y < xo +6 implica que A 
y — Xo > O se tiene que f(y) < £ (xo). 
Análogamente, se prueba que si f'(x0) < 0, entonces existe $ > O tal que xo — $ < 
y < xq implica que y — xo < 0, entonces f(y) > f (xo). 


< 0 y dado el hecho de que 


Nota. En el primer caso decimos que la función es creciente y en el segundo decimos que 
la función es decreciente, en una vecindad de xp. 


Teorema 2.2.2. Sea f : [a,b] > R una función derivable en xy € [a,b]. Si f tiene un 
máximo (respectivamente un mínimo) local en xy, entonces f'(xo) = 0. 


rueba. lor, 0 > 0 
Prueba. Del teorema anterior, tenemos que si f'(x 0, entonces f'(xo0) > 06 
f'(xo0) < 0, es decir, la función es creciente ó decreciente, entonces f no puede tener 
un máximo ni un mínimo local en xo. 
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Teorema 2.2.3. (Del valor intermedio) Sea f : [a,b] > R una función continua. Sean 
c=infí f(x) : x € [a,b] y d = supí f(x); x € [a, b]j. Entonces para cada z € [c, d] 
existe x € [a,b] tal que f(x) = 2. 


Prueba. Como f es una función continua sobre el conjunto compacto [a, b] se tiene que 
f' alcanza su máximo y su mínimo sobre el intervalo [a, b], llamamos a estos d y c res- 
pectivamente. Dado z € [c,d] sean Az = [x € [a,b] : f(x) < zp y Bz = [x € 
[a,b] : f(x) > z). Es claro que Az y B¿ son subconjuntos abiertos y disjuntos de [a, b]. 
Si fGo) X z para todo x € [a, b], entonces [a, b] = 4¿U B, esta es una unión de conjun- 
tos abiertos y disjuntos, y esto no es posible dado que [a, b] es conexo. Por lo tanto, existe 
x € [a,b] tal que f(x) = z. 


Teorema 2.2.4. (Del valor intermedio para la derivada) Sea f : [a,b] > R una función 
derivable. Si f'(a) < d < f'(b), entonces existe c Ela, b| tal que f'(c) = d. 


Prueba. Supongamos inicialmente que d = 0. Como f'(a) < 0 existe $ > O tal que si 
a<x=<a-+06, entonces f(x) < f(a). Análogamente, como f'(b) > O existe e > 0 tal 
quesib—e€ <x < b, entonces f(x) < f(b). Ahora, como f posee un mínimo en [a, b] 
existe c Ela, b] tal que f(c) < f(x) para cada x € [a,b], por lo tanto, f'(c) = 0. Si 
d H 0 consideramos la función g(x) = f(x)— dx y se tiene que g'(a) = f'(a)—-d <0 
y g (b) = f'(b)—-d > 0, luego existe c ela, b[ que satisface g'(c) = 0, esto es, 
fUc)=d. 


Nota. Aún cuando no los hemos indicado, hemos usado el siguiente resultado.Si f es 
derivable en xy entonces f es continua en xy. En particular, como hemos supuesto que 
f' es derivable en [a, b] se tiene que f es continua en [a, b]. 


Teorema 2.2.5. (de Rolle) Sea f : [a,b] > R una función continua en [a, b] y derivable 
en la, b|, la cual satisface f(a) = f(b). Entonces existe c Ela, b| tal que f'(c) =0. 


Prueba. Si f es constante sobre [a, b] no hay nada que probar. Por otra parte, si f no es 
constante sobre [a, b], entonces ella asume su máximo y su mínimo en [a, b], uno de ellos 
debe ser alcanzado en el intervalo abierto la, b[, pues f(a) = f(b). Es claro que para 
cada uno de esos puntos la derivada f se anula. 


Teorema 2.2.6. (del Valor Medio) Sea f : [a,b] > R una función continua y diferencia- 
ble en la, b[. Entonces existe c Ela, b| tal que F'(cMb—a) = f(b)— fía). 


Prueba. Definamos la función h por h(x) = (f(b) — fla))x — f(x)(b — a). Se tiene 
que h(a) = h(b) = f(b)a — bf (a), y dado el Teorema de Rolle existe un c ela, b[ tal 
que h'(c) = 0. Por lo tanto, f'(cMb—a) = F(b)— f(a), lo que completa la prueba. 


Teorema 2.2.7. (Taylor) Sea f : [a,b] > R una función de clase C"7! y n-veces deri- 
vable en el intervalo la, b[. Entonces existe c Ela, b| tal que: 


Jem 
(n—1)! 


(n) 
0=/0+106-0...+ oa AL ay. 
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Prueba. Sea h : [a,b] > R la función definida por: 


_ LD) 


h) = $0) £0)- $ (9b=)=->= py 


K 
A 


elegimos la constante K de modo que h(a) = 0. Claro que h es continua en [a, b], deriva- 
ble en la, b[, y h(a) = h(b) = 0. 


Notamos que: 


RC) == 0)= 000-3460 = PG = 0. 


PU) EE y 
K n—-1 
esto es, Pe 
ns A 
h'(x) = 1 (b=xy 71 


Por el teorema de Rolle tenemos que existe c Ela, b[ para el cual h'(c) = 0 esto implica 
que fU(c) = K. Como h(b) = 0 tenemos que: 


Fa) 


F(b) = fla) + FlaNb=a)...+ =D 


Or Eo ay, 
n: 


que es lo que queríamos probar. 


Sea f : [a,b] > R una función. Definimos f” (la composición de f consigo misma 
n-veces) como f% = Id (función identidad), f! = ff? = fof;¡f?= fof? y 
fr=fo0f7? paran > 4. 


Definición 2.2.2. Sea f : [a,b] > R. 
a) Decimos que xo € [a, b] es un punto fijo de f si f(x0) = Xo. 


b) Decimos que xo € [a,b] es un punto periódico de período n de f si los valores 
FG), F2(x0),..., F” (xp) están definidos, f*(x0) 4 xp, paral <i<n-—1ly 
F" (xo) = xo. 


c) La órbita de x € [a,b] por f (o la órbita de x bajo f ) es el conjunto orby(x) = 
LP" (x) : n e€ Nj, cuando tiene sentido la composición de funciones. 


Definición 2.2.3. Sea f : 7 > R una aplicación continua y sea xy un punto fijo de f. 
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a) Decimos que xp es un punto fijo atractor si: 18 > Otal que x €lxo — 8, xo + ó[—Lx0), 
entonces | f(x) — xol < |x — xol; 


b) Decimos que xp es un punto fijo repulsor si: 19 > O tal que x €]xo —8, xo + $[-(x0j, 
entonces | f(x) — xol > lx — xo]. 


c) Decimos que xp es un punto fijo neutro si no es atractor ni repulsor. 


Lema 2.1. Sea f : [a,b] > R una aplicación de clase C* y dos veces derivable. Sea 
xo Ela,b| un punto fijo de f con |f'(xo)| < 1. Entonces existe 8 > 0 tal que x € 
lxo — $, xo + 6| implica | f(x) — xol| < plx — xol, para algún p < 1. 


De este lema notamos que se cumple que | f£”(x) — xol < p”|x — xol, de donde 
deducimos que lim f”(x) = xo para todo |x — Xp] < 8. 
n—>00 


Prueba. Sead = f'(x0). Sabemos que para cada x €]a, b| existe t(x) ela, b[ tal que: 


PU) 


21 (x — 05 


Fx) = fo) + Ax —x0) + 


” 
Í An 
Sea ha) LLO) E 0) 
existe 4 > 0 tal que para todo |x — xp] < $ implica |h(x) — h(x0)| < e, esto es, 


f" (1) 
21 


. Por la definición de la función h vemos que dado e > 0 


tu=a sele =xol. 


Eligiendo e > O de modo que p = 4 + e < 1, obtenemos 


|fG6)—xol < plx—xol y |x—xo|<6G6. 


Proposición 2.2.8. Sea f : [a,b] > R una función diferenciable y xy un punto fijo que 
satisface | f'(xo)| < L entonces xy es un punto fijo atractor. 


Prueba. De la expresión 


109 = 106094 10090=x9 + LL 6x0), 


obtenemos 


O COCO [ y EE) =x0)| 


Fo) 
y luego 
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|fGe) — xo] < 1£'Gro) 11 + 8](< — xo)1, 


como podemos suponer | f'(xo)[1 + 8]| < p < 1; tenemos el resultado. 


Lema 2.2. Sea f : [a,b] > R una aplicación de clase C! y dos veces derivable. Sea 
xo Ela, b[ un punto fijo repulsor de f con | f'(xo)| > 1. Entonces existe 8 > 0 tal que 
x Elxo — $, xo + 6| implica | f(x) — xol > plx — xo], para algún p > 1. 


Prueba. Similar a la demostración del lema 2.1 


Proposición 2.2.9. Sea f : [a,b] > R una función diferenciable y xg un punto fijo que 
satisface | f'(xo)| > 1, entonces xy es un punto fijo repulsor. 


Prueba. Similar a la demostración de la proposición 2.2.8. 


Uno de los conceptos clásicos en dinámica es el de topológicamente conjugados. A 
saber: 


Definición 2.2.4. Sean X e Y espacios métricos, f : X > X y g : Y > Y aplicaciones 
continuas. Diremos que f es topológicamente conjugado a g si existe un homeomorfismo 
h:X > Y tal que goh(x) =ho f(x), Vx € X. 


Esto quiere decir, en particular, que hay un homeomorfismo que lleva la trayectoria de 
x bajo f en la trayectoria de h(x) bajo g. 
Otro concepto interesante en dinámica es el concepto de w-Límite: 


Definición 2.2.5. Sea f : 1! > I, I = [a,b] intervalo real. Para x € 1 el conjunto 
w-límite de x, bajo f, es: 


w(x) =fy € I/1n;> 00 con f(x) > y para i¡ > oo), 
es decir, es el conjunto de puntos donde la órbita positiva de x se acumula. 


Proposición 2.2.10. Si f : I > I yg : I' > I', son topológicamente conjugados, 
entonces: 


h(w f (xo0)) = 07 (M(xo)). 


Prueba. a)]Sea y € h(wf(x0)), entonces 3x € Wf(x0) tal que h(x) = y. Ahora, como 
existe existe sucesión 1; —> oo tal que fi (xo) > x; ya que h(£"(x0)) = g” (h(x0)) 
y h es continua entonces h(f*i(x0)) > h(x) = y. Esto es g”i(h(x0)) > y, por tanto 
y € Og(h(x0). Concluimos que h(wy(x0)) E og (A(xo). 

b)Sea ahora y € Wg(h(xp)); entonces existe sucesión 1; —> 00 tal que g”! (h(xo)) > 
y.Como h”! es continua, tenemos que h7*(g%i (h(x0))) = £* (x0) > h7* (o), por lo tanto, 
y =h7*(x) € 01 (x0)y, luego, x € h(w s (x0)). En conclusión: Wg (h(x0)) C h(w f(x0)). 
De a) y b) concluimos que h(wy(x0)) = 2 (M(x0)). 
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Corolario 2.2.11. Sean f : 1 > I y g : I' => I', topológicamente conjugadas por 
medio de h : 1 > 1'. Si x es punto periódico de f, entonces h(x) es punto periódico de 
2. 

Prueba. De la proposición anterior tenemos que si f : l > Iyg: I' => T' son 


topológicamente conjugadas entonces h(wy(x0)) = 0g(A(x0)), entonces dado que xy es 
punto periódico de f tenemos que 


(xo) = (xo, fxo). 00 Y). 


y h(cw 7 (x)) = (xo), HF (o), ++ MI" xo) = 03 ((x0) luego, (xo) es punto 
periódico de g. 


Observamos, entonces, que la dinámica topológica de dos funciones f : 1 => I y 
g: J > J, que son topológicamente conjugadas es, salvo un cambio de variables topo- 
lógico h : 1 > J que satisface ho f(x) = g o h(x), Vx € I; la misma. 


Para efectos de una conjugación topológica se debe exhibir una aplicación continua, 
con inversa continua que haga la ecuación de conjugación. Esto se puede debilitar solici- 
tando que la aplicación h sea sólo sobreyectiva (y no necesariamente inyectiva) y en este 
caso se llama una semi conjugación. 
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2.3 Resultados importantes de la Función Gamma 


En esta sección señalaremos algúnos resultados que nos permitirán hacer un estudio de 
la dinámica de la función Gamma. Entre los resultados que nos interesan están: encontrar y 
calcular los puntos fijos, determinar el carácter de ellos (atractor, repulsor o neutro) y/o los 
puntos mínimos o máximos de la función Gamma para x > 0 y x < 0. Dado que la función 
Gamma está definida como una integral impropia, no resulta simple encontrar, de manera 
explícita, estos puntos. Para ello, usaremos la aproximación de Stirling-De Moivre, pro- 
posición 1.2.32. Antes de encontrar estos puntos, veamos que ésta aproximación produce 


[0,0] 
valores cercanos a los que se obtienen con la función Gamma: T'(x) = 0) Pd, 
0 
[oe] 
Valor | I(x) = 1 ("1271 dt | Aproximación Stirling-De Moivre 

1 ¡ 0,9994994685336 

2 = 0, 8862269254528 0,8661553843017 

2 1 0,99997898067 

a = 1,3293403881791 1,3293307671298 

3 2 1,9999940099595 

2 == 3, 3233509704478 3,3233462787043 

4 6 5,9999956006056 

5 24 23,9999941451898 

6 120 119,9999881350628 

10 362880 362879,9971745871 


Cuadro 2.1: Tabla comparativa de valores de la función Gamma con la aproximación de 
Stirling- De Moivre de orden 5 usando Microsoft Mathematics. 


Gráficamente también podemos observar lo mismo: 


2.3.1 Punto de Mínimo de la función Gamma para x > 0 


Para calcular el punto mínimo de la función Gamma para x > 0 ocuparemos la apro- 
ximación de Stirling - De Moivre, es decir: 


1 1 139 571 
T(x) = VZAX THE 7* ( + ) 


12x + 28812 — 51840x3 248832011 


y también usaremos el método de Newton-Raphson, el cuál es un algoritmo para encontrar 
aproximaciones de las raíces de una función. 

Sabemos que I'(1) = (2) = 1, por tanto, del teorema 2.2.5 y usando el hecho de 
que ""(x) > 0 para x > 0, tenemos que el punto mínimo de la función Gamma está en 


2.3. Resultados importantes de la Función Gamma 91 


poc 
_-_——_—— r(r)=/ tota 
40 


a AZ E 1 l 139 571 y 
—.. Díx)==yY2rn1r*""1 1+3 ) 


2x 28812 5184014 248832071 ) 


a 


Figura 2.1: Función Gamma y aproximación de Stirling de orden 5 generada por Microsoft 
Mathematics. 


el intervalo [1, 2] y es el único mínimo para x > 0. 


Usamos la aproximación de Stirling-De Moivre y calculamos la primera derivada de 
ella. A saber: 


1 1 139 571 
T(x) = V2Ax* 720 ( + ) ; 


12x * 288x2  51840x2  2488320x* 


3 


3 1 1 139 
PQ) =V2xe* | 17% In(x) — 1 = = 
das ( 10) == ] ( +12 Y 28812 7 51840x3 


571 
2488320x* 


sl 571 
622080x5 ) ' 


y 1 1 139 
2 Xx—=>3 —X pa en 
)+ ARES ( 12x2 7 144x2 Y 17280xt* 
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entonces multiplicando y factorizando, obtenemos: 


3 5 
*=2 In(x x* 2 In(x 
OA OM 


T'(x) =V2re "lx 0721 
(o) Te ( n(x) + 12 288 


139x%=3 In(x)  571x*%7ZIn(x) AvÍ xr E Sara 
TE rc . 576 
973x732 571% 
103680 * 552960 ] 


Para encontrar el mínimo en el semieje positivo, x > 0, hacemos T'(x) = 0 y como 
V21me * 40 para cualquier valor de x, tenemos: 


03 In(x) de x*=3 In(x)  139x%3 In(x) 


¿ul 
AAA 288 51840 
STATE Inf) AT AT 5 E 973x173 
=a88320 2 "8 516 ” 103680 
571x%-%2 

552960 


multiplicando por 4976640x7*+%, nos da: 


4976640x* In(x) + 414720x* In(x) + 17280x* In(x)- 
— 13344x? In(x) — 1142x In(x) — 2488320x* — 1658880x* 
— 43200x? + 46704x + 5139 =0. 


Para encontrar las soluciones de esta última ecuación, utilizaremos el método de Newton- 
Raphson: 
a Fin) 


FU)” 


Xn+1 = *Xn 
donde la función, f(x), será: 


f(x) =4976640x* In(x) + 414720x* In(x) + 17280x? In(x) 
— 13344x? In(x) — 1142x In(x) — 2488320x* — 1658880x* 
— 43200x? + 46704x + 5139, 
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entonces 
f(x) =24883200x* In(x) + 1658880x? In(x) + 51840x? In(x)- 
— 26688x In(x) — 1142 In(x) + 4976640x* — 9538560x*— 
— 4959360x? — 99744x + 45562 


y dado que el mínimo de la función Gamma se encuentra entre 1 y 2 elegiremos xp = 
1,5, y obtenemos: 


== o = 1,4651788301846 
X0 
Fa) 
x2 = 101 ==, y = 1,4613788053233 
=D z rea = 1,4613361572095 
X2 
F(x3) 
xa = 13 = oy = 1,4613361518823. 


Dado el grado de aproximación obtenido tenemos que el mínimo de la función Gamma 
eS Xmin = 1, 4613361518823. Evaluemos este punto en la aproximación de Stirling- De 
Moivre, 

T(Xmin) = 0.8855226070359; 


por tanto, el punto de mínimo del gráfico de la función Gamma para x > 0, es: 


(1, 4613361518823; 0.8855226070359). 


2.3.2 Puntos Fijos de la Función Gamma en el semi eje positivo 


De la misma forma que calculamos el punto de mínimo de la función Gamma, calcu- 
laremos el(los) punto(s) fijo(s) de ésta para x > 0. Es decir, usaremos la aproximación de 
Stirling - De Moivre y luego el método de Newton-Raphson para calcular los puntos fijos 
de la función Gamma, en el semi eje positivo. 

Es claro que un punto fijo de la función Gamma es cuando x = 1, ya que (1) = 1. 
Lo complejo es calcular el otro punto fijo de ésta función. 

Primero, veamos en qué intervalo se encuentra éste punto fijo. Del cuadro 2.1 obser- 
vamos que 1"(3) = 2 y que T"(4) = 6. Esto es [3, 4] € T'([3, 4]) lo que nos indica que el 
segundo punto fijo se encuentra en el intervalo ]3, 4]. 

Ahora, usaremos la aproximación de la función Gamma para encontrar este punto fijo, 
es decir, resolveremos la ecuación: 


1 1 139 571 
V2x tig" ( + ) =xX, 


12x * 28812 — 51840x3 248832011 
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que es lo mismo que: 


1 1 139 571 
v2xr xx) =x 1 == <= = = 0 
a ES +2 *288x2  51840x2 2488320x4) * 


y, dado que es complicado encontrar soluciones de esta forma, usamos nuevamente el 
método de Newton-Raphson, 


8 (Xp) 
Xp41 = Xy == 
E Xn 
donde la función, g(x) es: 
1 1 139 571 
tl 
869 A + ix? 288x2 7 51840x3  2488320x1 


Derivando tenemos: 


1-3 In(x) E 03 In(x) 
12 288 


g (a) =V2xe * (+ In(x) + sá 


139132 In(x) 571x%73 In(x) A El 
NO A BO DE 
STE 973x571 

576 * 103680 ' 552960 )- 


y dado que este punto fijo se encuentra en el intervalo ]3, 4[ elegiremos xo = 3, 5: 


xa =x0- £60 3 5662572852031 
8'(xo0) 


xo =x1 ¿6D _ 3 5623978916007 


xa =x2. ¿6% _ 3 5623838217599 
2; 


xa =x E _ 3 5623828215736. 


por lo tanto, tenemos que el otro punto fijo de la función Gamma es: 
x = 3,5623828215736, en efecto g(x) = 3, 5623828215736. 


A continuación, veremos el carácter repulsor o atractor de los puntos fijos de la función 
Gamma. 
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Proposición 2.3.1. El punto x = 1 es un punto fijo atractor de la función Gamma. 


Prueba. Sea x = 1 y usando la función digamma, es decir, 


Tr 
W(x) = ——— 0,—1,... 
E +» (5- x+n- 1) ** 
ra 1 1 
tenemos que: YW(1) = -0 = -—y+ a ( = >) = —y, entonces tenemos que 
r 
w(1) = PT —y, luego (1) = y - (1) y como F(1) = 1, tenemos (1) = —y, 


esto es:|'(1)] = y < 1. Concluimos que x = 1 es atractor. 


Definición 2.3.1. Sea xp €la,b[ un punto fijo atractor de f :]a,b[> R. El conjunto 
estable de xp, bajo f, es el conjunto: W; (xo) = [y ela, bl; £"(y) Ale xo). 
Proposición 2.3.2. El conjunto estable del punto fijo atractor xy = 1, bajo la función 


Gamma, contiene al intervalo abierto |x2,x1[ donde 0 < x2 < 1 < xy; x1 es el punto 
fijo repulsor de T' situado en 13,4| y T(x2) = x1. 


Prueba. Primero probemos que |1"'(x,)| > 1. Sabemos que x, = 3, 5623828215736 es 
punto fijo de la función Gamma. Usaremos directamente la aproximación de Stirling-De 
Moivre, es decir, 


To) = V2rx* 207 “(1+ 


dE 1 139 571 
x  288x2  51840x3  2488320x* 


Derivando 


13 In(x) m 03 In(x) 


Tx) = V2re* (+79 + E 288 


139x%=3 In(x) 571x*32 In(x) A 
ETT TEST 24 
E 
576 * 103680 * 4976640 12 144 * 
139x%-2 57112 

17280 * 622080 ) 


+ 


por lo tanto, 
T"(x1) = 4,0025450750553, 


entonces tenemos que |"'(x1)| = 4,0025450750553 > 1. Concluimos que x, es un 
punto fijo repulsor de f”. 
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Sea ahora x €lx2,1[U]1, x1[. Primero suponemos que x2 < x < 1. En este caso 
tenemos 1 < TP (x) < x1 = IP (x2). 
Consideramos ahora x tal que l < x < x1; en este caso sea x €] 1, x¡[ el punto tal que 
T'(%) =0. Si y € [1, x], entonces (y) € [1(), 1 y 776) e [1, x]. 
En esta situación, para la diferencia: Ea 2(x) — 7 ?(1)|, tenemos: 


(0-10) = 1 M-0|=|1 Fw) Pm] F-11, 
y como sup(| 1 '(60|, | (PG, |P'1))< 1, entonces tenemos que |" (P(w)) P'(w)| < 
p < 1 y se cumple que |7?*(x) - PAD] < plx— 1] o sea |P4x)-1| < pix — 11, 
por lo tanto, 1 < 12(X) < Y. De esta forma obtenemos que 1 < ... < T(x) < 
LAY im Tk (x) = 1. En efecto, supongamos por absurdo que 1 < 
b= im FR: e ser "2(b) = b y tenemos |b—1| =|P?(b) - TAD)| < 


plb — e pao = maxí|P"(b)|, [PP (y), "(My < 1, lo que no puede ser. Por tanto 
cada y € [1?(x), x] satisface 
lim P"(y) =1 y [12(0),x] C W¿0. 


Observamos, ahora, que el resultado, enunciado en la proposición, sigue de observar que 
para cualquier w € [X, x1] Jk € N tal que 2%) < 17 (w) < Y. 


En efecto, aplicando lo anterior, tenemos: 


a) que six2 < x < x1 entonces lim P”(x)= 1. 
n—>00 


b) Si x €]x1,oo[, entonces lim /"”"(x) = oo. En particular x £ W¿(1). 
n—>00 


En efecto, se deduce de que |"'(x)| > |7''(x1)| para cada x > x¡ y de que x < 
T(O)<TUx<...Sib= lim T”(x), entonces lim |P"*Ux) - P"(x)] =0 
y no puede ser ya que |P"(P(xY)-P*60| =P» (DP G)=x| para z € [x, P(x)). 
Si escribimos p = min(|P'(w)|j; w € [x, (x)] > 4, entonces |" (P(x))— LO E 
DP (DEL) =x)1 > p”|P(<)—x|, por lo tanto, lim [PO -TP*(0|= 


Así, Vx > x1 valex < IT (x)<TU(x)<...y lim T*(x) =00 
n—>00 


2.3.3 Soluciones de la Ecuación I'(x) = -—2 
Para encontrar las soluciones de la ecuación (x) = —2, Wx € R — (0,—1,-2,...), 
lo primero que haremos será observar el gráfico de la función Gamma y la recta y = —2. 
Vemos que la recta y = —2 corta el gráfico de la función Gamma en los intervalos 


]-2n—1,-—2n[ paran e N,n > 1. Primero nos concentraremos en el intervalo ] — 3, —2|. 
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Figura 2.2: (x) = 2. 


Para encontrar estos valores usaremos la fórmula de reflexión de la función Gamma: 


TGO)TU-=x)= enn Vx EZ. 


Sea x < O con I'(x) = —2, tenemos que: 


T 
ASAS sin(7rx)' 


haciendo el cambio de variable y = —x, tenemos: 


21 (y +1) = 7 


obtenemos entonces la ecuación: 
2 (y + 1) sin((ry)) —= 71 =0. Q.1) 


Ahora, como y > 0O,podemos usar la aproximación de Stirling-De Moivre para y. A saber: 


1 1 139 571 
Ty) 2 y 2e7Y (1+ ) 


12y + 288y2 — 51840y3  2488320y* 
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Sustituyendo y + 1 por y en esta expresión, tenemos: 
O +0=yP(v) 


1 1 139 571 
e /2my +27) (1 
A + 12y Y 288y2 — 51840y3  2488320y1 


reemplazando estos valores en la ecuación 2.1, tenemos que resolver la ecuación: 


1 1 139 
22 y +7 si 1 S — 
ES A ( +12y Y 288y2 7 51840y 


571 
te | =0, 
2488320y* 


Para resolver esta ecuación, usamos el método de Newton-Raphson con g(y), dada por : 


1 1 139 
=N Ty YET ET si 1 E e 
280) ny?*Ze sintay) + y + 288y2 — 518403 a 
371 
2488320y* : 


para su derivada tenemos: 


2 (y) =2 21 y +2 e7Y (16) +1+ 5) sin(7 y) ( + a 


12y 
1 139 571 , 
+88y2 7 51840y3 says) sintry)(1 hi 
yd 139 571 ] ES 
12y ' 288y2  51840y3  2488320y* 


O A 139 571 
TT COST 
il 12y ' 288y2  51840y3  2488320y* 


7 e Do 571 
—= SIn( TT ds 
AC 12 y2 7 144y3 * 17280y4 — 622080y5 


, 
] a 2'On) 
Las raíces, que buscamos, las obtenemos usando la recurrencia )»+1 = Yn 


20m) 
Ahora, dado que hay dos soluciones en el intervalo x €] — 3, —2| y que —x = y, entonces 


las soluciones (para y) las encontraremos en y €]2, 3[. Probaremos, en primer lugar, con 
yo tal que —xp está más cerca de la asintota x = —2, esto es, usaremos xy = —2, 25. Esto 
es: yo = 2,25, lo que, luego de algunas iteraciones, nos da que: 


y = 2,2192646352. 


Ahora, para encontrar la segunda solución, hacemos la recurrencia con un x más cercano a 


2.3. Resultados importantes de la Función Gamma 99 


la asíntota x = —3, esto es xy = —2, 75, o sea empezaremos la recurrencia con yo = 2,75, 
lo que, luego de algunas iteraciones, nos da: 


y = 2,9047877691, 


es decir, que las soluciones de la ecuación 21 '(y + 1) sin(xry) — 1 = 0 en el intervalo 
]2, 3[ son: 


y1 =2,2192646352 

y2 2, 9047877691 
como x = —y tenemos que las soluciones de la ecuación T'(x) = —2 para x €] — 3, —2 
son : 


x¡ =—2,2192646352 
x2 =—2,9047877691. 


Ahora, el siguiente paso será calcular | ”'(x)| para los valores de la soluciones P'(x) = —2 
en el intervalo ] —3,—2[ , para ello usaremos la fórmula de reflexión de la función digamma, 
es decir, usaremos la igualdad: 


W(l—x) — W(x) = 7 cot(7x). 
1) Para x = -2,2192646352: 
YW(l—x)-W(x) = mx cot(mx), 
entonces 


w(1+2,2192646352) — W(-2, 2192646352) = 7 cot(-2, 2192646352 - 11) 


Yw(-2,2192646352) = Y (3,2192646352) — 7x1 cot(—2, 2192646352 + 71). 


Para calcular W(3, 2192646352) usamos la aproximación de la función digamma: 


1 1 1 
0 mOl=> =D aer 


por lo tanto, 
Y (3, 2192646352) = 1,0058746602, 
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entonces 


W(-2,2192646352) = 1,0058746602 + 3, 8154458674 
= 4,8213205276 
T"(22,2192646352) 
T (22, 2192646352) 
T'(22,2192646352) == 4,8213205276 - T(-2,2192646352). 


= 4,8213205276 


Ahora para calcular P(-2, 2192646352) usamos la fórmula de recurrencia de la fun- 
ción Gamma varias veces, es decir: 


_PG+D a T(x +4) 


200  GCTEDGHDEEN 


entonces tenemos que: 


T'(1,7807353648) 


T(22, 2192646352) = 70. 4632116915 * 


Ahora para calcular 1"(1, 7807353648) usamos la aproximación de Stirling-De Moivre 
para x > 0: 


í 1 139 571 
TG) =vV21x* 2e *|1 
(o) = V21x* Ze ( Fix? 28812 51840x3 soi) 


cuyo valor es: 
T(1,7807353648) = 0, 9263778450, 


entonces 


T'(1,7807353648) 


—0, 4632116915 
_ 0,9263778450 
— —(,4632116915” 


T(2,2192646352) = 


por lo tanto, 
T(22, 2192646352) = —1, 9999016905 


y con ello tenemos que: 


T'(22,2192646352) =4, 8213205276 - —1, 9999016905 
= — 9,6421670738. 
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2) Para x = —2,9047877691: 
YW(l—x)-W(x) = 71 cot(1x), 
entonces 


w(1 +2, 9047877691) — W(-2, 9047877691) = 77 cot(-2, 9047877691 - 11) 


Y(-2, 9047877691) = W(3, 9047877691) — xr cot(—2, 9047877691 - 71). 


Para calcular Y (3, 9047877691) usamos la aproximación de la función digamma: 


E - 1 
AA” o 


por lo tanto, 
Y (3, 9047877691) = 1,2287259217, 


entonces 


W(-2, 9047877691) 1, 2287259217 — 10, 1877321892 
= — 8, 9590062675 
T'(22, 9047877691) 
T (22, 9047877691) 
T'(2, 9047877691) == — 8, 9590062675 - (2, 9047877691). 


— 8, 9590062675 


Ahora para calcular (2, 9047877691) usamos la fórmula de recurrencia de la fun- 
ción Gamma varias veces, es decir: 


AD) TO+9 


(+ DG +2) +3)" 


P(x) 


entonces tenemos que: 


T(1,0952122309) 
—0, 4766509504 * 


T (22, 9047877691) = 


Ahora para calcular 1"(1, 0952122309) usamos la aproximación de Stirling-De Moivre 
para x > 0: 


Tx) = V2mx* The ( + + 


1 1 139 571 
12x  288x2  51840x3  2488320x* 


cuyo valor es: 
T(1,0952122309) = 0, 9529790219, 
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entonces 


F(1,0952122309) 
—0, 4766509504 
_ 0,9529790219 


——(), 4766509504” 


T (22, 9047877691) = 


por lo tanto, 
T(22, 9047877691) = —1, 9993226093 


y con ello tenemos que: 


T"(22, 9047877691) = —8, 9590062675 - —1, 9993226093 
= 17, 9119437875. 


Para ambos valores: 


x= —2,2192646352 
x= —2, 9047877691 


tenemos que |"'(x)| > 8. 


Para encontrar las soluciones de la ecuación P"(x) = —2, para x €] — 5, —4[, usaremos 
nuevamente las ecuaciones (2.1), (2.2) y (2.3). Las raíces, que buscamos, las obtenemos 
_£0n) 

20) 
tervalo x €] — 5, —4[ y como y = —x, entonces las soluciones están en y €]4, 5[. Primero 
probaremos con yg más cerca de la asintota x = —4, por ejemplo con xp = —4,25, o sea 
en nuestro caso con y = 4,25, nos da, luego de algunas iteraciones, que: 


usamos la recurrencia: Ya+1 = Yn . Ahora, dado que hay 2 soluciones en el in- 


y = 4,0202211103. 


Ahora, considerando un x más cercano a la asíntota x = —5, por ejemplo xy = —4,75, 0 
sea, en nuestro caso con y = 4.75; tenemos, luego de algunas iteraciones, que: 


y = 4,9958032768, 


es decir, las soluciones de la ecuación 21" (y + 1) sin(ry) — xi = O, en el intervalo ]4, 5| 
son: 


=4,0202211103 
4, 9958032768, 


y 
y 


N BP 
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pero como x = —y tenemos que las soluciones para la ecuación I'(x) = —2 son: 


x1 =—4,0202211103 
x2 =—4, 9958032768. 


Ahora el siguiente paso será calcular | "'(x)| para los valores de la soluciones '(x) = -2 
en el intervalo ] —5, —4[. Para ello usaremos la fórmula de reflexión de la función digamma. 
A saber: 
W(l—x) — W(x) = rm cot(7x). 
1) Para x = -4,0202211103: 
W(l—x)-W(x) = mx cot(mx), 


entonces 


w(1 + 4,0202211103) — W(-4, 0202211103) = 77 cot(-4, 0202211103 . 7) 


Yw(-4,0202211103) = Y (5, 020221110) — 77 cot(—4, 0202211103 - 7). 


Para calcular YW (5, 0202211103) usamos la aproximación de la función digamma: 


1 Ñ 1 
120x4” 


1 
W(x) = In(x) — De 


por lo tanto, 
Yw(5,0202211103) = 1,5105833571, 


entonces 


w(-4, 0202211103) =1, 5105833571 + 49, 3867259555 
50, 8973093126 
T'(-4,0202211103) 
T(24,0202211103) 
T'(=4,0202211103) =50, 8973093126 - T(-4,0202211103). 


50, 8973093126 


Ahora para calcular P"(-4, 0202211103) usamos la fórmula de recurrencia de la fun- 
ción Gamma varias veces, es decir: 


Pica en TOO 


ox (a+ Dr +2) + 3) + 4) +5)” 


P(x) 
entonces tenemos que: 


T'(1,9797788897) 
—0, 4958096386 * 


T (24, 0202211103) = 
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Ahora para calcular I"(1, 9797788897) usamos la aproximación de Stirling-De Moivre 
para x > 0: 


1 1 139 571 
TO) e Vaz * (1 
0) e ( +12? 288x2 518403 Sus) 


cuyo valor es: 
T”(1,9797788897) = 0, 9915966971, 


entonces 


T(1,9797788897) 
—0, 4958096386 

_ 0,9915966971 

——(, 4958096386” 


T(-4,0202211103) = 


por lo tanto, 
T (4, 0202211103) = —1, 9999544581 


y con ello tenemos que: 


T'(-4,0202211103) =50.8973093126 - —1, 9999544581 
=— 101,7923006650. 


2) Para x == —4, 9958032768: 
W(l— x)- W(x) = m cot(xx), 
entonces 


w(1 +4, 9958032768) — W(—4, 9958032768) = 77 cot(—4, 9958032768 - 11) 


Y(—4, 9958032768) = W(5, 9958032768) — 7 cot(—4, 9958032768 - 70). 


Para calcular Y (5, 9958032768) usamos la aproximación de la función digamma: 


1 1 1 
AA A 


por lo tanto, 
w(S5, 9958032768) = 1,7053565003, 
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entonces 


Y (—4, 9958032768) =1,7053565003 — 238, 2673359476 
= —236, 5619794473 

T'(-4, 9958032768) _ 

T(=4,9958032768) 

T'(-4, 9958032768) = — 236, 5619794473 - T(—4, 9958032768). 


— 236, 5619794473 


Ahora para calcular P(-4, 9958032768) usamos la fórmula de recurrencia de la fun- 
ción Gamma varias veces, es decir: 


T(x+7) 


COS FDAIDETDET DATA 


entonces tenemos que: 


T'2.0041967232) 


P (54, 9958032768) = 0 o0s908827 * 


Ahora para calcular 1"(2, 0041967232) usamos la aproximación de Stirling-De Moivre 
para x > 0: 


í 1 139 571 
TG) =vV21x* 2e *|1 
(o) = V2x* Ze ( Fix? 2882 — 51840x3 soi) 


cuyo valor es: 
T(2,0041967232) = 1,0017607198, 


entonces 


T(2,0041967232) 


—0, 5008908827 
_ 1,0017607198 
——0,5008908827 


T(-4, 9958032768) = 


por lo tanto, 
T(=4,0202211103) = —1, 9999579837 


y con ello tenemos que: 


T'(-4, 9958032768) = — 236, 5619794473 - —1, 9999579837 
= 473, 1140194355. 


Para ambos valores: x = —4,0202211103 y x = —4, 9958032768 tenemos que |”'(x)| > 


106 2. Sistemas Dinámicos 


En resumen: en los intervalos ] — 3, —2[ y ] — 5, —4[, las soluciones de TP (x) = —2, 
tenemos que: 


x W(x) Tx) 
—2,2192646352 | 4,8213205276 —9, 6421670738 
—2,9047877691 | —8,9590062675 17, 9119437875 
—4,0202211103  50,8973093126 | -—101,7923006650 
—4, 9958032768 | -236,5619794473 | 473, 1140194355 


Cuadro 2.2: Tabla de valores para las soluciones de la ecuación 7'(x) =-—2 en ] — 3, —2| 
y |] — 5,41. 


2.4 La Entropía Topológica 


2.4.1 La definición de Adler-Konheim y Mc Andrews 


La definición de entropía que usaremos se debe a Adler, Konheim y McAndrew (1965). 


En lo que sigue M X 4 y t será una topología en M tal que (M, 7) es un espacio 
compacto. 

Por la compacidad de (M, 7) sabemos que si 4A C 7 es un cubrimiento por abiertos 
de M, entonces existe A” C A, que es un subcubrimiento por abiertos de M y que posee 
un número finito de elementos. Esto es la cardinalidad de 4”, denotada 4+(4”), satisface 
1 <+H(4”) < oo, con H(4”) e N. 

Así, podemos definir el número 


N(A) = minftt(4”) : 4 C A es un cubrimiento finito de My. 
Definición 2.4.1. Diremos que el subcubrimiento 4' C A es minimal, si: 
a) H(4”) < oo. 
b) No existe A” C A subcubrimiento finito tal que H(4”) < H(A”). 


Definición 2.4.2. Se define la entropía del cubrimiento abierto A € 7 como el número 


H(A) = log(N(4)). 


Notamos que N(A) = eH4W > 1, por lo que H(A) e [0, ool. 
En lo que sigue S(M) = [A C T; A es cubrimiento de My denotará el conjunto de 
todos los cubrimientos por abiertos de M. 
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Definición 2.4.3. Para A,B € G(M), se define la junción de 4 y B como AvB = 
[ANB:AEA,BeB). 


Es claro que la junción es una función V : G(M) x S(M) => G(M), v(A,B) = 
Av B. 

Además, AC AV A y vale A v A = A sólo si ocurre que los elementos de 4 son 
dos a dos disjuntos y el cubrimiento A contiene al vacío. Caso contrario, la inclusión será 
estricta. 


Definición 2.4.4. Diremos que un cubrimiento B e 6(M) es un refinamiento del cubri- 
miento A € G(M ), denotado A < B, si WB € BJA € A tal que B C A. 


Definición 2.4.5. Un cubrimiento por abiertos 4 se dice nbd-finito si todo x € X tiene 
una vecindad V,, tal que V, MN A 4 4 sólo para un número finito de elementos A € A. 


Propiedades de v y < 


1. Sea Ap = ([M). Es claro que Ay € G(M) y que VA € G(M) se tiene que 
AV Ay = Ay VA = A, es decir, Ay es el elemento neutro de la junción. Además, 


Ap <A, VA E CM). 


2. Se cumple que AvB=Bv A, VA,B e G(M). O sea la junción es conmutativa 
en 6 (M).Esto sigue de AN B = BN A, para cualquier par de conjuntos A, B. 


3. Se cumple que Av(BvC) =(AvB)vC, VA,B,C € 6(M). Esto sigue del 
hecho que AN(BNC)= (4N B)NC para conjuntos 4, B, C cualesquiera. Así, 
la junción es asociativa en C(M). 


4. VAe G(M) vale A < A.0 sea, la relación de refinamiento es refleja en C(M). 
5. Sean A,B,C € G(M) entonces, si A < B y B < C resulta ser 4 < C. O sea, la 


relación de refinamiento es transitiva en C(M). 


Prueba. En efecto, sea C € C. Como B < C existe B(C) € Btal que C € B(C) = 
B.Como A < B existe A = 4(B) € A tal que B C A, luego C C B C A implica 
C CA. Así, VC EC existe A= A(C) € Atalque C CA A<C. 


Observacáo. a) Como < satisface las propiedades 4 y 5, entonces constituye un 
casi orden reflexivo entre elementos de £(M). No es un orden parcial pues 
no cumple la propiedad de antisimetría. Por ejemplos, los cubrimientos 


A= E (o. >) n € N] U(B(0, 2)) 


0 [o(a2)oen 
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10. 
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de la bola cerrada de centro O e ¡R” y radio 1 satisfacen A <ByB< Ayno 
son iguales entre ellos. 


b) La relación de refinamiento no es simétrica. Esto es claro, ya que VA € C(M) 
vale Ay < A para Apo = (M). 


. SIA<A y B< B', entonces A vB < 4! v Bi. 


Prueba. Sean A' € A', B' € B' y consideramos 4' N B”. Como A < 4! existe 
A = A(4”) € A tal que 4' € 4. Como B < B' existe B = B(B”) € B tal que 
B' C B. Luego 4NB" C ANB.ASiVWA NB' € A! v B' existe ANBEAVB 
tal que AN B" CAN B. Esto implica que 4 vB < 4! v B'. 


. Sean U,B € G(M). Como Ap = [My < B;B < B;U <U se tiene Ay VU < 


U v B, osea, VU, B € G(M) vale! <U vB y B < U v B. De aquí concluimos lo 
siguiente: para cada par U,B € C(M) existe A € G(M) tal que ll <AyB<A. 
En efecto, basta tomar A = U v B. 


. SIU,B € GC(M) yU < B entonces N(U) < N(B) (y consecuentemente H(U) < 


H(B)). 


Prueba. En efecto, sea B' = (B1,..., By) un subcubrimiento minimal de 5. Como 
U < B, para cada ¡ existe U; e U tal que B; C U;. Se sigue que U! = [U;,..., Up) 
es un subcubrimiento de M. Así N(U) < k = N(B). 


. SU < B entonces N(U v B) = N(B) y H(U v B) = H(B). 


Prueba. Como B < U v B, entonces de la propiedad 8, tenemos que: 


N(B) < NU v B). 0) 


La hipótesis 4 < B dice que para cada B € B existe U = U(B) e U tal que 
B C U(B). En particular BCE BNU(B)..WB € B, existe BN U(B) e Bvu 
tal que BE BNU, es decir, U v B < B. Aplicando lo probado en 8: 


N(U v B) < N(B). E 


De (*) y (**) concluimos que N(U v B) = N(B). 


ParatodoU1,B € G(M) vale NU VB) < N(UNIN(B) (y consecuentemente H(U V 
B) < HU) + H(B).. 
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Prueba. Sean los conjuntos 4' = (U1,..., Una), B' = (B1,..., By(s)) subcu- 
brimientos minimales de 4 y B, respectivamente. Se tiene que: 


UB = (0 OB MU) = o MB) 


es un subcubrimiento de 14 v BÉ, por lo tanto, N(U v B) < NUIN(B). 


Consideremos ahora p : M > M una aplicación continua. Para cada UU € G(M) sea 
py (U) = [p"HU); U e U). Es claro que p”U) € C(M). 


11. 


12. 


13, 


14. 


Sil/ < B entonces y UU) < p MB). 


Prueba. Para cada B € B existe U € U tal que B C U. De aquí: 


y UB) UU) => y UU) <p (8). 


SiU,B e G(M), entonces UU v B) = y UU) v y UB). 
(Esto es, la aplicación p7* : E(M) > G(M) es un (€ (M), v)-morfismo). 


Prueba. Sigue del hecho que pH(4AN B) = y”*(4)N p”*(B) para cualquier par 
de conjuntos A, B CM. 


Para todo U € £(M) vale N (pg7*(U)) < NU). 


Prueba. Sea U' = [Uy,...,Un(u)) un subcubrimiento minimal de 14, entonces 
y HU”) = [p7(U1),....p”*(Unq))) es un subcubrimiento de p”* (14). De aquí 
que N (p97*(U)) < N(U). 


Sip: M > M es sobreyectiva, entonces N (9 (U)) = N(U). 


Prueba. Sea u, = (pr (U1l),... pa (Uno! wn); un subcubrimiento minimal 
de y”*(U/). Tenemos : 
Now) 
Mec U eu) 

j=1 
Entonces: 

Nou) 

amcel (UY eu) 


j=1 
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N z N 5 y 
Como y (a To )) C UNETO y (paUJ)) plo" UU) CU) 
y p es sobreyectiva, se EL que: 
NU) NU) 
AMm=Mmc U eletu»y)= U 0. 
j=1 j=1 
Esto es: (0, y... Uy (o! on) es un subcubrimiento por abiertos de 14, por lo tanto: 


NU) SN (gu). 


Combinando esta desigualdad con la desigualdad de la propiedad 13 tenemos el 
resultado. 


15. Para cada UU € G(M) se cumple que: 


n—1 
oa (Y vu) = lim Zu votumv.. ve Dg) 


n>00 n 
k=0 
existe y es finito. 
Para demostrar este resultado, probaremos lo siguiente: 


Lema 2.3. Si (an)nen es una sucesión de números reales no negativos que satisface 


K Ñ An 
An+m < An + Am; Yn,m € N, entonces existe im pd 


Prueba. La condición 44 +m < An + Am implica que Ax. < ka, y así 


1 pa An 
— a —. 
ER kn = > 


Haciendo k > oo concluimos que c = inf (2) < En Vn eN. 


*0 < c+e. Así, para n > ny escribimos 
no 
A 0 < q < np. Luego 
An ía Anop+4 < PAny + Ag < Ano men Ag < 
n nop +4 p (no + 2) no+75  Pio+q4q 
P 


Aq 1 
E A <c+te+-suptaj;0 <i<mnpoj. 
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E An An y 
Así para n grande — < c + 2e. Como c < — < c + 2e concluimos que c = 
n n 
z An 
lim — 
n>00 n 


Prueba. (De la propiedad 15) Sea 


k=0 


n—-1 
dn = H (Y vu) =H (Uva Y) 
tenemos: 
nilo =.5 (U vou)... ve "Du 
VU) VTA) Y go) 
=.5 (U vo HUM vw. vo DU 
vo UV WU) V...v O) 
De acuerdo a la propiedad 10, tenemos 
-1 —(n—1) 
ant+m SH (UV € UV... vo (U)) + 
+H (07 (u vou) V...v DU) l 
Usando la propiedad 13, tenemos 
PS (u vo UU) V...v e "=DU)) ES 
+H (u vou)... v DU) 


es decir, Ay +m € An + Am- 


De la definición de (ar )nen se tiene que 4, > 0. Así, obtenemos que existe 


1 n—1 
lim: 2: dim: 28 (V E) 
n>00 n n>002N 


k=0 
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Definición 2.4.6. Para cada U € G(M) y y : M > M continua, la entropía de y 
respecto de l/ es el número: 


n—-1 
ho = pm, 30 (Veo). 
k=0 


as = H(U); d, < An.1 <N-A, Se tiene La, < a1, esto es: 


16. Para cada UU € G(M) y y: M > M continua, se tiene que: 
h(p,U) < H(U) = log(N(U)). 


Prueba. Como a; = H(U) y An = An.1 < NA], Se tiene La, < a. Por lo tanto, 


An y . 
— < a, tomando, ahora, el límite cunado n > oo, tenemos el resultado. 
n 


17. SIU,B € G(M) satisfacen U < B, entonces h(p,U) < h(p, B). 


Prueba. De la propiedad 11, tenemos que UY < B, entonces (U) < p*(B), 
Vk € N. Aplicando la propiedad 6, obtenemos: 


Uv Mv. vo “DU <BVg UB. vo DB). 
De la propiedad 8, obtenemos que 


n—1 n—1 
a (Voto) <n (Voto). 


k=0 k=0 


Luego, h(p,U) < h(p, B). 


18. SiU € S(M) y y : M > M es un homeomorfismo, entonces h (y, U) = 
h(p9*,U). 


Prueba. De acuerdo a la propiedad 14, 


N(oU) =N (9 (p(U)) = NU). 
Así, 


HU pu) NI) = H(9"(u vo WU) V...v 


da 8 geD(u))) = H(u vVoTUU) VW... DU). 


dividiendo por n y pasando al límite, cuando n —> oo, se tiene el resultado. 
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Definición 2.4.7. Sea y : M —> M una aplicación continua. Se define la entropía topo- 
lógica de y como h(p) = supth(p,U) :U € C(M)). 

Nota. Observamos que 0 < h(p) < oo. 


Definición 2.4.8. Diremos que una sucesión de cubrimientos abiertos (U4,,;n € N) es 
refinante si: 


1. UU, <Un+1, Vn € N. 
2. Para cada B e 6 (M) existe U4,, tal que B < UY. 
Lo interesante de una sucesión de cubrimientos abiertos refinante, es el siguiente lema. 


Lema 2.4. Si (14, ) es una sucesión de cubrimientos abiertos refinantes, entonces h(p) = 
lim A(p, Un). 

n—>00 

Prueba. La propiedad 1 implica h (p,U1,) < h (p,Un+1). Así que: 


lim 2 (p,Un) = sup th (p,Un):n EN. 


La definición de h(p) implica que para cada e > O existe B¿ € G(M) tal que h(p) < 
h(p, B¿) + e. La propiedad 3 dice que existe 14, tal que B¿ < Uy. De aquí h(p, B¿) < 
h(p,U,) y luego tenemos que: 


h(p) < h(p,Un) +es h(p,Un+p) +e, VpenN. 


Ya que h (9, Un+p) <h(p) <h (9, Un-+p) +e, Vp € Ñ, debe ser entonces que h (y) = 
lim + (p0,U,). 
n—>00 


2.4.2 Cálculo de la Entropía en Espacios Métricos 
Suponemos que (M, d) es un espacio métrico compacto. 


Definición 2.4.9. 1. Para A C M, se define el diámetro del conjunto A como d(4) = 
supíd(x, y): x,y € A). 


2. Parall € G(M) se define d(U) = suptd(U); U eu). 


Observacáo. Sil/, B son cubrimientos abiertos de M tales que 4 < B entonces d(B) < 
d(U). En efecto, VB € B existe U(B) € U tal que B C U(B) entonces d(B) < d(U(B)). 
Así d(B) = supíd(B); B € Bj < supíd(U(B)); B € Bj < d(U). 


Probemos, ahora, el siguiente resultado: 
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Lema 2.5. (del cubrimiento de Lebesgue) Sea U € G(M), entonces existe e(U) > 0 tal 
quesix € M y B,(x) = [y € M;d(x, y) < r) satisface r < e(U), entonces JU € U 
tal que B,(x) CU. 


Prueba. Como U € G(M) existe subcubrimiento finito 
fU1,..., Un). Para cada ¡ = 1,2,...,n,sea f£ : M > [0,oo[ la función continua 
f(x) = d(x, MNU;). Como (MA U;) es cerrado, se cumple que f¿(x) =0 => x € 
(MNU,). 

Como para cada x € M existe ip € [1,2,...,nj tal que x € U;,, entonces 


Fi) (x) > 0. Q.4) 


Sea f : M > R la función f(x) = max[ f(x), ..., fh(x)). De acuerdo a (2.4) se 
tiene que f(x) > 0Vx € M.Como f es una función continua, entonces F(M) = [81, 82] 
es un intervalo cerrado con $1 = inf[ f(x); x e Mj > 0. 

Sea x € M cualquiera y 0 < 4 < 6,. Se cumple f(x) > 6. Esto es, hay /, € 
[1,2,...,nj tal que f(x) > 6. Esto implica que d(x, MU) > $. Luego, para cada 
y € Bs(x) debe ser que y € U;, pues d(x, y) < 6. 

Concluimos que Bs(x) € U,,.. 


Corolario 2.4.1. Para cada U € G(M) existe e = e(U) > 0 tal que si A es un conjunto 
abierto con d(A) = supíd(x, y); x, y € A) < €, entonces existe U € U tal que A CU. 


Prueba. Notamos que A C Bg(ay(x) para cada x € A. De aquí basta considerar d(A) < 
$1 como en el lema de cubrimiento de Lebesgue. 


Observamos que si 0 < $ < e(U) entonces $ satisface la propiedad. A saber, si A C M 
es un conjunto tal que d(A) < 8 entonces existe U € U/ tal que A C U. En efecto, como 
$ < e(U) sirve el mismo U del corolario 2.4.1 . Esto nos permite la siguiente definición: 


Definición 2.4.10. Se define el número de Lebesgue de 14 como: 
$(U) = sup fe(U); e(U) satisface el corolario 2.4.1]. 


Corolario 2.4.2. Sea U,B € G(M). Sea d(U) = supíd(U); U € U). Si ocurre que 
d(U) < $(B), entonces B < U. 


Prueba. Para cada U € U vale d(U) < $(B), entonces IB(U) € B tal que U € B(U), 
esto es, B < l4. 


Corolario 2.4.3. Sea (U) € G(M) una sucesión de cubrimientos abiertos tales que: 
LM 
2. (Un) >0,n >00. 
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Entonces SU, y es una sucesión refinante. 


Prueba. La primera propiedad es la primera parte de la definición de una sucesión refi- 
nante. Si B € G(M) es un cubrimiento abierto, entonces $(B) > O y ln € N tal que 
d(U,) < 68(5). De acuerdo al corolario 2.4.2, B < Uy (que es la segunda propiedad de un 
cubrimiento refinante). 


En particular, hropl(p) = lim h(p,U4,) para cualquier y : M —> M continua y 
n—>00 


cualquier sucesión de cubrimientos (14, ) como en el corolario. 
Veamos algúnos ejemplos de sucesiones que satisfacen el corolario anterior. 


1 
Ejemplo 2.4.1. Si consideramos U/, = 4 B,(x);x € M,2r < —¿ el conjunto de todas las 
n 


bolas abiertas de diámetro d = 2r < —, entonces la sucesión (14, ;n € Nj es refinante. 


1 1 
En efecto, si B,(x) € Un+1, entonces 2r < El < — y By(x) € Un, por lo tanto, 
n n 
Un < Un+1- ¡ 
Por otro lado d(U4,) = supíd(U); U € Uy = — y así d(Un) > 0,n > 00. 
n 


Ejemplo 2.4.2. Sea (B,) C G(M) una sucesión de cubrimientos abiertos tales que 
n 


d(B,) > 0,n > 00. Sea Uy = V Bx. Como B,, < UU, entonces d(U,) < d(Br) y 
k=1 

luego, d(Uy) > 0, n > 00. 

Además, como 14, < 14, +1 el corolario 2.4.3 garantiza que (14, ) es una sucesión refinante. 

O sea, a partir de cualquier sucesión de cubrimientos abiertos (B,] C G(M) tales que 

d(Bx) > 0, n > oo, se puede construir una sucesión refinante. 


Ejemplo 2.4.3. Sea o : Xy > Xp, 24 = (0; No >(0,1,...,n — 13) el conjunto de 
sucesiones de n- símbolos. Se tiene (ver Labarca (2011), sección 1.5) que hop([o) = 
log(n). 


2.4.3 Extensión de la definición de entropía a espacios no compactos 


Hasta ahora nuestra definición se aplica a espacios métricos compactos. Es claro que 
aplica, también, a espacios topológicos compactos. Vemos, ahora, la definición que usare- 
mos para nuestro caso y que se refiere a espacios no compactos. 

Antes de nada recordaremos, brevemente, algunas definiciones que vienen de la topo- 
logía (ver Dugundji (1966)) 


Definición 2.4.11. Consideremos un espacio topológico (X, 71). Una familia de conjuntos 
D = (Doa C YlO € T) se denomina vecindad-finito si cada x € X tienen una vecindad 
abierta V, C X tal que V, N Dy 4 Y para, a lo más, una cantidad finita de índices a € Z. 
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Definición 2.4.12. Un espacio topológico (X, 1) se dice: 


* de Hausdorff'si separa puntos( para cada x, y € X existen abiertos U,., U, tales que 
x € Ux y y € U, y se cumple que Uy N Uy, = 8); 


» segundo contable si la topología tiene una base enumerable; 


de Lindelof si todo cubrimiento por abiertos tiene subcubrimiento enumerable; 


* paracompacto sitodo cubrimiento por abiertos tiene un refinamiento abierto vecindad- 
finito; 


» localmente compacto si cada x € X tienen una vecindad V, tal que V, es compacto; 


* O-compacto si es localmente compacto y puede expresarse como una unión enume- 
rable de conjuntos compactos; 


Para nuestros efectos (X, 1) será un espacio de Hausdorff, segundo contable, Lindelóf, 
paracompacto y o -compacto. Por ejemplo: la recta real satisface todas estas propiedades 
y cualquier espacio euclidiano, también. 


Definición 2.4.13. Se define C¿(X, X) como el conjunto de funciones f : X > X tales 
que existe K € X compacto para la cual la restricción f|x : K —> K es una aplicación 
continua. 


Nota. Observamos lo siguiente: 


a.) Para cualquiera f € C¿(X, X) y para cualquier compacto K C X donde la restric- 
ción f|x : K > K es continua, se puede calcular hyop(f, K); 


b.-) La unión finita de conjuntos compactos en un espacio Hausdorff es compacta; 


c.-) Sea f : X > X una función y suponemos Ky y K2 son conjuntos compactos tales 
que £|x, y f|x, son continuas entonces f|x,ux, es continua y, en este caso, 


hrop[f, Ki U K2) = maxi hroplf, K1), hropl f, K2)). 


En lo que sigue, y para efectos de espacios (X, 1) no compactos, consideraremos la 
siguiente definición de entropía 


Definición 2.4.14. Para f € C¿(X, X) sea 
Ky =(K C X| Kes compacto y f | x es continua; . 


Se define 
hropl[f. X) = sup Mhsop(f, K)IKe Ks) : 
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2.4.4 Entropía topológica e hiperbolicidad en la Función Gamma 


En esta sección probaremos que la entropía de la función Gamma es infinita y que 
su dinámica contiene un subconjunto invariante hiperbólico. Para ello, probaremos otros 
resultados primero. 


Proposición 2.4.4. Para cada n > 1 se cumple que: 
T'"(x) <0 para cada x € ]-2n —1,-2n|. 
Prueba. Usando la función digamma tenemos: 


E. Y 
O E 3) 


n=1 


derivando con respecto a x, tenemos: 


POP) Oé y 1 


Too (x+n- 1) ds 


n=1 


entonces > 
69-69 MO 


0 
[Goy 


y como [T'(x)]? > 0 obtenemos que: P“(x) - F(x) — [Go]? > 0, lo que implica 
T"(x) - T(x) > [0 (o)? y, como [TF Go]? > 0, entonces T”(x) - T(x) > 0. Ahora, 
dado que x € ]-2n — 1, —2n[, se cumple PF (x) < 0 y, por lo tanto I"”"(x) < 0; como 
señalado. 


Observacáo. a) Dela proposición anterior concluimos que /"'(x), para x €] — 2n — 
1, —2n[ es decreciente, para cada n > 0. 


b) Dado que W(x) = == derivando tenemos que: 
1" , TJ" 2 
Aya O TO-IFO, 


[Por 


de este hecho podemos concluir que la función digamma es creciente para x > 0 y 
para x € ]-n—1,—n[conn e ÑN. 


Teorema 2.4.5. Para cada n > 1 sean -2n— 1 <-—=x3 <—=x] <-—y3 <—y] < —2n 
puntos que satisfacen: I' (-x3) = TP (-1p) = -2n-1yTP (=x7) E P (-y3) = 
—2. Esto es: ([=x3,-x7]) = [-2n-1,-2] y P ([l-13,-y?)) = [-2n — 1,-2), 
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entonces se cumple que: 


min (min re ell: 


(- risks >> 1 


para algún 4 > 1. 

Prueba. Dividiremos la demostración en dos partes: 
a) [MEyz|>8y 

d) [MEN] > 16. 


Antes de calcular la derivada en estos dos puntos, vemos otras relaciones que nos 
servirán en esta demostración. 

Consideramos, para cada 1 < k < n, los puntos —2k — 1 < —xk < yk < —2k. Ya 
vimos, en la sección 2.3.3, que: —y) + —2, 2192646352 y que Y (—y]) = 4.8213205276. 
De este hecho tenemos que W(—y1) > 4 y dado que P(—y1) = —2 entonces: 

T(= 1 
Y(—y¿) = LO» > 4> Ty!) <-8, 
por lo tanto, 
MD] > 8. (2.5) 


También, vimos que —x! == —2.9047877691 y W(—x!) = —8.9590073177. Luego, tene- 
mos Y(—x1) < —8 y dado que T(—x!) = —2, entonces: 


T"(= 1 
YW(=xb = pe <-—8 > TP (xj) > 16, 
por lo tanto, 
|PExp| > 16. (Q.6) 


Ahora, dada la fórmula de recurrencia de la función digamma tenemos que: YW(y) = 


1 
W(y + 1) — —. Aplicando esta relación a y = —x — 1 tenemos que: V(—x — 1) = 
y 


1 
W(=x) + + ; aplicando la relación a y = —x — 2 tenemos: 
x 


1 1 1 
1 O A, de 
o o il e 


Sucesivamente, obtenemos la relación: 


W(=x —2k) = W(- 0 + 2 + 


A 2.7 
x+1 med eme 6) 
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parax>0yx¿Z”. 
Ahora ya con las relaciones (2.5), (2.6) y (2.7) probemos los puntos a) y b). 


a) Probemos primero que J” (-y —2k) > —2, con esto concluíremos la desigualdad 


—=y1-2k <—y HH para1 < bs < n. Enefecto, tenemos: (—y! pe € ]- o —2k — 2| 


y= y € ]-2k — 3, —2k — 2[; es decir, ambos valores: (—y y! — 2k) y mi I están 
en el mismo intervalo. 


Probamos esta afirmación por inducción: 


* k = 1. En efecto, para obtener el valor: N(—y! — 2) usamos la fórmula de 
recurrencia de la función Gamma, y tenemos que: 
Mer=0 1) 
Tey,-2= == * Z 


2-2 Er-Deyx-0D' 
y dado que 1” (—y]) = —2, entonces 


E 
Ey -DEy2= 0" 


Como —y] == —2.2192646352 tenemos y! < —2 y, entonces que y! =2< 
—4e —y! —1 <-3, por lo tanto (—y! — 2) - (—y1 — 1) > 12, entonces 


T(=y¿-2 = 


=D 
(SA) 


> —2, con lo que, Tí=y =2) > 2. 


Suponemos ahora que: P(—y!—2k) > —2 y probemos que la relación se cumple 
para k +1. 


Usando la fórmula de recurrencia de la función Gamma, tenemos: 


P=y2 2) 


T (=y¿ —2k + D) E (—y! —2k — 2) - (—y1 —2k —1) 


y dado que —y] < —2 entonces: —y—2k-2 < —-2k-4 y —y¿—2k-1 < -2k— 
3; y luego tenemos que: (— y! —2k- )-(Ey¿—2k—1) > Qk + 4)Qk + 3) > 12 
y, por lo tanto: 


T(ey¿ —2k) a 
¡TN 
2 
> 
M2) :(=y1 M1) 


PEy2 Uk +1) = 


> —2; 
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b) 


como queríamos probar. 


De este hecho podemos concluir que, puesto que: P(— y! —2k) > —2, entonces se 


cumple que —y! —2k < —yF* para 1 <k<mn. 


Ahora dado el hecho de que —y!—2k y — qa están en el el mismo intervalo, ] -2k — 3, -2k — 


para 1 < k < n, y usando el hecho de que la función digamma es creciente en 
]-2k — 3, —2k — 2[, ver observación , concluimos que: 


Voy -2k)< Uy) 1<ksm 2.8) 
y como W(—y1) < W(—y! —2k) por (2.7), de esto obtenemos la siguiente desigualdad: 


YW(—y3) < W=y¿ —2k) < W=y HF) => W(—yh) < (y A+) 


/ k+1 
1 k+1 po E 
y como Y(—y,) > 4, entonces Y(—y, **) > 4 lo que implica F ;Fr > 4 y dado 
—=2 


Tie) 
que My +) = o 
por lo tanto, Vk > 1 vale |P' (yk +1) > 8. Esto concluye la demostración de la 
afirmación a). 


—2, obtenemos que > 4, entonces Pisyirs) < —8, 


k> 
Para resolver esta parte, primero probemos que (=xk +2k + 1) —3 0. Para ello 
veamos cuánto es el valor de A Tí=k + 2k +1). 
=>00 


Usando la fórmula de recurrencia de la función Gamma, tenemos que : 
Meda) = aa er 
parak € N y como TE = -—2, entonces 
Fe 4241) =-2> (Ple + 1) (=xf + 2k) 
para k € N, además en ésta última expresión: le + DE +2)--- (xk + 


2k) contiene (2k + 1) —términos y cada uno de ellos es negativo, por lo que el valor 
ME + 2k + 1) será positivo. 


Veamos ahora que el término (—xf + 2k), con —1 < —xf + 2k < 0 no se acerca 2.0. 
De hecho —1 < —xF + 2k < —0, 9. Hacemos esta demostración por inducción. 


+ k = 1. En este caso —x] = —2, 9047877691, por tanto 


=x! + 2 = —0, 9047877691 < —0, 9. 
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+ Suponemos ahora que —1 < —xf + 2k < —0,9, y probemos que se cumple 
para k + 1. Para ello probemos antes que se cumple que e AE —xk . Para 
probar esto último, basta con probar que I'(— e cd +2) <P (=xt ) = —2 ya que 
ambos valores están en el mismo intervalo -2k=1 <—xk,=xFFH1 4-2, < 2k. 
Por inducción: 


— Para k = 1. En efecto Ex +2) = DEA + PAR =x*) = 
(x2) (7 — 1) y como 4 < x? < 5, entonces xi > 1yx?—-1> 1, ee 
tanto, 


Mex? 4+2) < Tex?) = 2. 

— Suponemos ahora que vale para k y probemos que se cumple para k + 
1. Usando la fórmula de recurrencia de la función Gamma, tenemos que 
HEEE SNA A EA) ANA) 
y como xF+2 > 1 yx*+2-—1 > 1, tenemos que T(=xF* + 2) < 
Ptexttl) — = —2, por lo tanto se cumple para k+1.Porlo tanto "(=x AL 

2) < Ma) = —2, entonces A +2< ON 


Ahora: A +H42k+1 = (in +2) +2k <—x! + 2k y como se cumple 
para k, tenemos que ad + 2(k + 1) < —0, 9. También, como —2k — 3 < 
e, entonces —1 < —xF+1 4 2k +2, por tanto —1 < —xF+1 4 2(k +1) < 
—0, 9 y se cumple para k + 1. 


Entonces tenemos que: 
M2 += 2 NA e DEA + 2) + (al + 2h), 
por lo que tomando límites, tenemos que: 


lim T(=xÉ +2k +1) =00 
k=>00 


Ahora, como —2k — 1 < —xf < —2k, entonces O < =xk +2k+1<ly Tx + 
k>00 . k k=>00 

2k + 1) ——> oo, por lo que necesariamente tenemos que: (—x% + 2k + 1) —= 0, 

dadas las propiedades de la función Gamma. 

Puesto que —2k — 1 < —xFF1 4-2 < —x*, entonces 0 < PH Uk E DA 1 

—xf + 2k + 1. Esto es, si €; = —x1 +2j + 1, entonces O < €x+1 < €x, Vk > 1 

Estudiemos ahora la diferencia: Y(=xk) - te 


Para ello usaremos la representación en serie de la función digamma, a saber: 


e 1 
E A O e 


n=1 n=1 
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Tenemos: 


00 
1 1 
w(—x* y il ax E 
( 1) E )= a 2 k+1 


$ xk —y ll =.mn : 


Ahora analicemos la serie (2.9) por partes. 


(2.9) 


Sabemos que 2k < A <2k+1yque2k +2 < E 1 < 2k +3, es decir, que para n, 
k k+1 


desde n = 0 y hasta n = 2k, las diferencias xf —n y xi —n, son positivas. Además 
tenemos que E < ae entonces ze =n< dl — n, lo que significa que: 
a . 0,1,2 2k 
< para n=0,1,2,...,2k. 
o E E 
E y : oa yk k+1 
Para n, desde n = 2k + 4 hasta n = oo, las diferencias xi — 1 y Xx; '” —n son 
negativas, es decir, que n — 5 yn-— a son positivas y como xt > A 
1 1 E 
tenemos n — bj >n- e +1 y entonces Es +1 lo que significa que 
+1 
N=XÍ NX] 
1 1 
sn E para n =2k+4,2k +5,...,00 
On  Xí—N 


Para n = 2k + 1,2k + 2 y 2k + 3, la diferencia EN — n es negativa; y la diferencia 
> —n es positiva paran =2k +1 yn = 2k +2, y es negativa para n = 2k + 3. 


Ahora, agrupando todos estos casos, tenemos que la serie (2.9) se puede descomponer 
de la siguiente forma: 


+ + 
a. A? AS Q.10) 


ES 1 1 Es 1 1 
4 de 
2) ld) 


n=0 n=2k+4 


Tenemos entonces que la suma de n = O hasta n = 2k es positiva y que la serie desde 
n = 2k +4 hasta n = oo también es positiva; luego solo tenemos que ver los términos 
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desde n = 2k + 1 hasta n = 2k + 3. Es decir, tenemos que estimar el valor de: 


1 1 1 
== + == 
xE=2k=1 xFFI-=2k-1  xf-2k-2 
1 1 1 
AH 2 2-3  xbt-2k-3 


(2.11) 


Si el valor en (2.11) es positivo, concluimos que la diferencia W(—xk) — W(=x +1), 
es positiva. 


Sumando los términos en (2.11), tenemos el siguiente denominador: 


GF -M=D=<4 7" =2=138 (7 =9=2% 
((__XE>E>XE- >= -]XAAS<]á, > o o po 
<0 >0 <0 
x (a — 2k — 2) x (x% — 2% — 3) x (01 — 2k — 3). 
— pm) pa — pa) 
>0 <0 <0 


Concluimos que el denominador es mayor que cero. 


Ahora analizamos el numerador de (2.11), el cual es la suma de las siguiente expresio- 


nes: 

GE —2k 1) x (xÍ —2k 2) x (a+ — 2 — 2)x 

l 1 1 Q.12) 
xQ F=2k=3)x (1 F7=2k%=3) | 
OM DAR MDAC ja) 
(4 =2=3) 0" 2% =3) 
E=W=1 0" =D 00% =%=Dx 

; Ñ 1 (2.14) 
x («E —=2k 3) x (a E —2k — 3) : 
= (1 -M=1) (7-2 -1= De 2.15) 
«GF =2=3007 2-3) 

k k+1 k 

xi —=2k-1)x =2k-1)x — 2k — 2)x 

(e xd Ix ) 2.16) 


(ico) <=DE=5) 


124 2. Sistemas Dinámicos 


===) (<= Do -1= 2D) Dia 
2-2) =0=3), 


Esto es, el numerador de (2.11) es igual a la suma de las expresiones (2.12), (2.13), 
(2.14), (2.15), (2.16) y (2.17). Ahora, sabiendo que —xF*! 4 2 < —xf, entonces 


A =2k-3 > + —2k — 1; luego notemos que las expresiones (2.13), (2.14), (2.15) 
y (2.16) son mayores que las expresiones: 


= (1-2 =D0-=2=24-2)0F=2=30-%=D (118 


V-2M=D0 =D <= 2=3)6=2b=1) (2.19) 
=(x=2% = 1) =2% +DG0Í=2k=2GF -2k=3)1G6É=2k=1) (2.20) 
(é=2=D0F=M 4 D0 == F=200=2%=M (021) 
respectivamente. Notamos, además, que al factorizar la suma de las expresiones (2.18), 
(Q.19), 2.20) y C.21) es igual a: 
60 === =1=/D)06=%Mk=14 42) 


de esta forma, tenemos que la suma de las expresiones (2.12), (2.13), (2.14), (.15), 
(Q.16) y L.17) es mayor que: 
(0. = 1) x 0-2 =2) x (77? 2 =2)x 
«== =3)2 0" ==) 60 =0=0% (2.22) 
=%=1= 4D GF=2=M4 3/2): 


Para seguir usaremos ahora las igualdades: ex = —xk +2k + l;€x41 = == + 
Ak+D+1y Edd = e — 2 = €x — €k+1. Así, reemplazando estos valores en la 
expresión (2.22) tenemos: 

(Q — Ex +1 Eex — 1)( — Ex 41 ex — 2)(6x — Ex 41)+ 

+ Geek — V2)N=Ex + 12) = Méx — €r+1)+ 

+ 2(—4ej — SerEx+1 + 3Eí41) — ER — CeñEn+1 + DK ER 1 

A E E TS 
(Ex + €x+1) 


= (€ — Ek+41) (+- 6l€x + €x+1) — 2€x 
(Ek — €Ex+1) 


+ oler — €x+1)) >>0. 
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Notar, en relación al tamaño de los ez, que: como al = 2, 904787769, xl = 4, 995803277, 
xi = 6,999900774 y xf = 8,999998622, entonces : e, = 0,095212231,€2 = 
0, 004196723, €3 = 0,000099226 y e4 = 0, 0000013. 


Entonces, hemos probado que (2.11) es mayor que cero, lo que implica que la expresión 
(2.10) es mayor que cero, por tanto, 


[,0) 
1 1 
k k 
YU) wat) => (5 - a.) > 


n=oX%1 79M Xx] 


De aquí tenemos que: Y(=xk) > rd Vk > 1. 
Aplicando sucesivamente esta desigualdad, obtenemos: 


=> tE>- > ben. Er 
y puesto que Yw(=x!) < —8, sección 2.3.3, concluimos que: 


PEA 
————— < -—48 y, luego 
Mex) aa 
<a: 16, 
De (2.5), (2.6), a) y b), resulta que: 
minímin([T(=x4)] > 16,1" (yDl>8)1<k<n3)>1> 1; 


para algún A > 1. 


Nota. Podemos considerar, para efectos de cálculo, 14 = 8. 


Corolario 2.4.6. Para n > 1 sean Jy = [=2b,—x+] y Ig = [yk]. Se cumple 
121 > A > 1; para cada z € Jg U lx paral <k <n. 


Prueba. Primero notamos que si z e Jz entonces |” (z)| > [PM exbl y que si w € lg 
entonces |” “(w) > |T Ye yk )|. El resultado sigue ahora del Teorema 2.4.5 


[0,0] n 

Ahora, para cada n € N sea Ay = A rai (0 JU 1) Esto es, 
j=0 ¡=1 

n 

xE€ An => TM Uel]sUL.VjENo. 


i=1 
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Denotemos L¡ =Í¡,L2 =J¡,L3= bh, L4= La, .... Los+1 = Ls+1, Los+2 = Js+1 
para cada O < s <n-— 1. Tenemos 2n intervalos (L1, L>,..., Lonj y además 


min (|P'(2)),7 € L¡UL2U...UL2anj>4>1. 


00 2n 
Observacáo. A, = A ra (Ú .) 


j=0 ¡=1 


2n 
ZE An => relt; para cada ¡ =1,2,...,2n. 


i=1 


Asociado a cada z € Ay podemos definir la sucesión 0(z) € X2n = [0 : No > [1,2,3,...,21 
por: E 
ADO =j =T(WeLj. 


Proposición 2.4.7. 0 : A, > 23, satisface a o 0(z) =0 o FP (z) para cada z € Ap. 


Prueba. Sea 0(z) = (00, 01, 02,...), entonces 


[0,0] [0,0] 00 
z€ A PT (Lg) = Lo, N A EA E DAS (A ra,,)) : 


¡=0 ¡=1 i=0 
[0,0] 
Así, tenemos que: z € Lg, y T(z) € MíZo T(Lo,,,).La condición T'(z) € A PHLe, 1) 
¡=0 


implica que 0(T'(z)) = (01,02, 03,...), por lo tanto, 
drob(z)=00T(z) 


como anunciamos. 


Proposición 2.4.8. 9 : A, > 23, es una biyección. 

Prueba. Suponemos z, w € A, satisfacen 0(z) = 0(w). Esto quiere decir que O(z)(1) = 
0(w)() para cada ¡ € No. 

Por tanto (7 *(2), 1 *(w)) C Lema) En particular | *(z) — TF (w)| < Laroa| para 
cada ¡ € No. 

Como |P*(2) — T(w)| = (17 Go)l|z — w| con x € [z, w], entonces 


[0 (2) -T*(w) > Amena) + Aema—olz— wl 
> (minfAi,...,AnP' |z — w] > A |z — wi; 


por tanto sii > 00 y z 4 w vale |I''(z) — T'(w)| > oo cosa que no puede ser pues 
P(z)- 1"(w)| < 1si (1 (2), 1"(w) € Loco) por lo tanto, debe ser z = w y la 
(910) 
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aplicación es inyectiva. 


Vamos a probar que O : A, > 2, es sobreyectiva. Para ello sean 


2n 


PT (2lz € U Ly En = max 


i=1 


2n 


IMorze UL; 


i=1 


Pn = min 


> 


2n 
se cumple €, > |'(Z)| > Pn > 14 >1Vz € UL; 
¡=1 


Sea 0 = (00, 01,02,...) € Xzn y consideremos 


A(00) =Lo,, 
A(00, 91) =L0, A EE 
A(00, 9, 92) =L00 A TAEa) A TP” Lo,) 


k 
A(O0,01,....0x) = () 17 (Lo,)- 
i=0 


Es claro que 4(00) > A(00,01) >... A(00,01,...,0x). 
Probemos que |4(00,01)| < p;'. En efecto, si Ag, = [x(00), y(00)] = Lo,. Sea 
A(00. 01) = [x(90, 91), y (00, 91)] C [x(00), y (00)], entonces: 
|T (x (00, 91)) — T (y (00, 01))| = |T'(2)| lx (00, 01) — y (00, 01)| 
para algún z € [x(00, 01), y(0o, 01)], luego: 
Lo, | =|T')|14(60, 01)1 > pr 14(00, 91), 


y por tanto, p;, "| Lg, | > |4(00, 01)! 
De aquí, concluimos que: p,,* > p, *|Lg,| > |4(00, 01)|, como anunciamos. 


Probemos ahora que 4(00, 01,02) C A(00, 01) € A(00) satisface 15 > |4(00, 01, 02)|. 
Sea A(00, 01,02) = [x(00, 01,02), y (00, 01,02)] C [x(00, 01), y (00, 01)], entonces 
[72(x(00, 01, 02)) — 21000, 01, 02))| = |Lo,| = 


= (17?) (2)| |x(00, 01,02) — y (00, 01, 02)| = 
=|P (DP): TD 14(00, 01,02) > p51 A(00, 01, 02)|. 
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entonces 
1 > |Lg,| =|T?(x(00, 01,02) — Ty (00, 01, 02))| > ph 14(00, 01. 02)! 


por lo tanto, 
Pz" > |A(00, 01, 02)| . 


Sucesivamente para 
A(00,01,...,0x) C A(00,01,...,0x-1) CC... C A(00, 01) C A(00) 


se cumple que |4(00, 01, ...,0x)| < py *. Así tenemos una sucesión de intervalos encaja- 


dos 
A(00,01,...,0x) Cc A(00,O1,...,0x-1) (A A(00) 


k> 
tal que |4(00,91,...,09)] == 0. 
[0,0] 


Así 8 A(00,01,...,0x) = (x(00, 01,02,...)) = un punto. Claro quex = x(0p, 01,... 
¡=0 
satisface 
x € Lo, T(%) € Lo, TU) E Lo,,..., 1 (%) E Lo, 
y, en consecuencia; 0(x) = (00,01,...). 
Ahora resulta claro que: 


07(00.01.02...)=(]T“Lo,). 


i=0 


Con lo que la aplicación resulta ser sobreyectiva. 


Usando Ó podemos dotar a 2, de una topología 
AC YX, es abierto <> O71(A) C An es un abierto. 


Considerando la topología inducida en 2, tenemos que 0 : Ay > 2, es un homeo- 
morfismo que conjuga | 1, : An > An eno : Xan —> 2p y, por tanto Arop ("| 4, ) = 
log(2n). 


Corolario 2.4.9. hsop(I') = oo. 


Prueba. En efecto; es claro que A¡ € A € A3C€...C€ An € An+1---y Sl Ax = 
[0,0] 00 

(0 pa U ]-2; — 1,—2j[ |, entonces cada An E€ Av y Pla, : An > An es parte 
¡=0 j=0 

de la dinámica de Pl1. : Ano —> Aoo. Por tanto, 


hrop[Ula.) > hrop( la, ) =1log(Qn) para cada n > 1, 
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es decir, hrop( | 1.) = 00. 


Observacáo. Recordamos que un conjunto A C R, invariante para la función f : R > 
R se dice hiperbólico si ocurre que existe 4 > 1 tal que |((£”)(x)] > A” para cada x € A. 
El corolario 2.4.9 tiene como consecuencia que el conjunto A,, es hiperbólico para la 
función J/”. 


Observacáo. Hasta aquí nuestro estudio de la función Gamma. Quedan diversas pregun- 
tas abiertas para continuar con esta línea de trabajo. Por ejemplo: Si consideramos 


00 n—-1 
An= YPF U (j-1-27D]: 
k=0 j=0 


entonces: ¿Será que la dinámica de /” restricta a A, contiene uno o más atractores?. Con 
seguridad, este es un problema no trivial e invitamos a los interesados a estudiarlo. 


Notamos que todo nuestro estudio, sobre la dinámica, tiene que ver con la versión real 
de la función Gamma. Es claro que muchas preguntas se pueden hacer para la versión 
compleja pero claro...eso es otro mundo....a ver si por ahí hay alguien que se motiva y nos 
sorprende con resultados novedosos en ese ámbito.... 


Ejercicios 


3.1 Representaciones 


NN ¿70 qy = E 
0 2 a 


2 


1. Verifique que: 


Solución. Usando el cambio de variables 4u* = w, tenemos la integral: 


Foo 1 Foo 1 1 1 1 / 
/ ¿CU gy = al w 2e "dw= T = a 
0 2/0 Jo 2/0 2 2Va 


2. Verifique que: 


F ¿0 dy = e 
0 2 a” 


3. Seaa,c E R* yb € (1, +00), pruebe que: 


+o0 1 
/ pt z- 
0 c(a In(b))+< 


3.1. Representaciones 131 


1 
we"ldw 


c(a In(b))+ 


Solución. Usando el cambio de variable b7%!" = e7%_ entonces dt = 
tenemos que la integral: 


1 
AS pde Ez l pe ae Laa ES P la) 
0 0 


c(aIn(by= claln(by)z 


4. Seaa € (1, +00) y b,c € RF, verifique que: 


ES dto 1 
o b+ctt > ppreiioismz)” 
5. Seaa,c E RF y b € R—(40), tal que 4 = b? — 4ac < 0, pruebe que: 


po e ¿au?+buv+co2) dy dy == 21 ] 
0 0 vV—A 


6. Cuando x € (0, 1), verifique el valor de las siguientes integrales: 


(a) qe cos(1*)dt = cos (2) : ab +) 
0 


2x Xx 
En A 
(b) 1 sin(1*)dt = sin (57) . 7). 
+90 sin (1 + E) _asin(2 43) 
(c) 0 a U= OO 


7. Six € (0, 1), verifique que: 


+oo px=1 3 
dt = —xó cot(rx) - esc(srx). 
LS (rx) <sc(rx) 


8. Six € RAIZ” yn € Z¿, verifique la representación integral de Cauchy-Saalschítz 
valida para x € (—n — 1,—n). 


+o0 HF (yn 
T(x) =/ pal (+ - y ( z d Ja. 


n=0 
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9. Se define el doble factorial de un número entero como: 


n-(n—-2):---5:3-1 sin es impar 
ni=4n-(n—2)---6:-4-2 sin es par 
1 sin =-—1,0. 
Pruebe que: 
2n-1 a 
aras 3) - aa 


1 ED"2" 
Tr == == 

(b) (=n +3) Qn — pia 

1 
(c) r(u+>3) or (+ ¿)=e 1". 
(a) Usando la proposición 1.2.8, tenemos que: 

1 1.3-..-(2n—1) 
Tr — 
( + 5) E yT, 


Usando la definición de doble factorial tenemos: 


1 Qn — 1)! 
(b) Usamos la fórmula de recurrencia de la siguiente forma: 
+1xr Polop ad 
= =>: X = = = = <= 
n+> n+> a 
r + z 
= —Nn = 
2 


Prueba. 


Por lo que tenemos: 
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Usamos el hecho que 


Il 


r(on+)) e) 


Sucesivamente, tenemos que: 


=- 2n+1 q yn 
r(=1+3)= T(=n+ 2%) PM a 


2 (n+Dn4+3) 5) > (2n— DM 
2n 


(c) Se obtiene multiplicando las expresiones del ítem a) con las del ítem b). 


10. Sea a, 6 € R — [0], pruebe que: 


m Tax) BP 
x>0I(Bx) a 
Tr 1 T 1 
Solución. Podemos escribir (ax) = ER mismo que P(Bx) = SS 
ax x 


entonces tenemos que: 


mA T(ax) si Bx de Mar: Po 4D 8 


i = lim li = 5 =“, 
x>0 P(Bx) x>00x x>0IT(Bx+1) a T() «a 
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11. 


12. 


13. 


14. 


3. Ejercicios 


Pruebe que: 
1 


1 
k e 
-/ x*In(x)dx = AFD? 


para k > —1. 


Solución. Usando la sustitución x = e*, tenemos que la integral: 


1 0 
-/ x* In(0)dx = -f e4+D! pp 
0 00 


ahora usando el cambio de variable (k + 1)t = —w, tenemos que la integral anterior 
queda de la forma: 


0 +00 
-/ eE+D! ¿gr = + / we " dw = pe) = : , 
Era (k+ 12 Jo (k+ 12 (k+1) 


Pruebe que: 
+o0 
yá 5 
dx=T | - 
Ja ala 


Solución. Usando el cambio de variable £ = x?*, tenemos que la integral: 


ES 4 ES P (2) 5 
o = -— tf E = 4 =YFPl- 3 
/ e * dx ¿/ e “dt 4 (5) 


Pruebe que: 
Pax) 1 


5-90 Mx) a 


Solución. Del ejercicio 11, tomamos los valores « = a y $6 = 1, entonces: 


tim Te) _B_1 
x>0 IT (Bx) a a 


Pruebe que: 


im (04m) TG) = 


(eb 
nt 
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Solución. Usando la fórmula de recurrencia de la función Gamma sucesivamente, 
tenemos que: 


T 1 
im (e +n) Jl (x) = Jim (x Ey dt ET == —= e +) 
TO 
(ED 
por lo tanto: 
4 Ey 
po a EG) Ñ n! Í 


15. Pruebe que la ecuación '(x + 1) = k, k 4 0 tiene infinitas soluciones reales. 


16. Para a > 0, pruebe las siguientes igualdades: 


(a) ys x2s+1 ¿ax? gy = 67D ya 1): 
0 2as+l 


1 
+00 E T(s+ 7) 


Solución. Para ambos problemas usamos el cambio de variable ax? = u, entonces 
tenemos que: 


uE 1 20 r 1 
a) / x?25+1 ga gy = 2] ue “du = Pa) 
0 0 


2qs +1 2as+l 
+o0 00 1 
to: a PT ls +5 A 
b) ll ed —/ u” 2e "du = at como señalado. 
0 2a:+2 Jo 24 


17. Expresa, en términos del factorial de un número, el n—ésimo coeficiente de la ex- 
pansión de (1 + x)? y(+ 72. 


Solución. Al usar la expansión de Taylor para ambas funciones, tenemos: 


ÍA 1 5 
1 2=1 a A as 
Ue A 
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; 3 5 35 
1 73 =1- e e A 
Mein a a” 


4 o e 1 
entonces tenemos que n-ésimo término de la expansión (1 + x)2 es: 


Qn — 3)! 


20D an 2)" para n=2,3,... 


An = id 


PE de A e _1 
y que el n-ésimo término de la expansión (1 + x)72 es: 


Qn—1)! 


2 yr PRA 1,2,... 


An = (-1)" 


18. Considere la expresión de la función Gamma dada por el producto de Weiertrass: 


verifique que ésta representación cumple con las condiciones del teorema de Bohr- 
Mollerup. 


Solución. Verificamos que ésta representación de la función Gamma cumple las 
condiciones del teorema de Bohr-Mollerup. A saber: 


(a) TO =1 
(b) P(x+1) =x1 (x). 


(c) In(P'(x)) es convexa. 


En efecto, usando la expresión (1.11), tenemos que: 


a) 
i d E 
POD= ki (14 45+...+ $ —In(n)) k 
ES Me 

k=1 
=- mi ¿am , 1.2. ”=1, de 

. 12 n—1 n 

= lim n Cer j 

n>00 2.3 n n+1 


por lo tanto: (1) = 1. 
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ds 1 
e Kk,tenemos: 


b) Usando el hecho que e”? = lim n II 


oro) E: 
er a+1) Eds 
run, Il a 
P(x) e—yx T] em 
e n=1 (1 a 2) 
yl en (1 + =) 
eE Ñ n=1 
ox+1 2 x+1 
I+ ) 
n=1 n 
xe Y Fl (1+2 en 
di A 


; E 
además como e Y = lim n II ek, entonces tenemos; 
n—>00 


A A A AN 
PQ) x+1l n>0o Fa pnl 


Xx E x+2 x+n-1 x+n |- 
pS 


= lim n . id 
x+2 x+2 x+n 


n>0o0 x+1l 


= lim n——— = 
n>0o0 x+n>+1 


T(x+1 
Luego, > = x, por lo tanto, 


+1) =x1(x). 
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c) Aplicando logaritmo natural al producto de Weiertrass tenemos: 


era E eñ 
| A M4) 


Lx Ea x 
== yx —In(0) +97) Im(1+2) 
n=1 n=1 


l l 
y dado el hecho que y = lim (1 + > +...+ a In(1)), entonces tenemos 
n—>00 


In(PG)) = 
= lim DT DN - Eno) 
cas k=1 k=1 k=1 


ahora derivando dos veces a In(1'(x)) tenemos 


E [0,0] 
qx UP ()]= 3 + aX eE > 


y dado teorema 1.2.13 tenemos entonces que In(7(x)) es convexa. 


3.2 Función Beta 


1. Verifique las siguientes identidades: 
(a) B(x, y) = Ba +1, y) + B(x, y +1). 
(b) B(x, y): BH y, 2) = B(y,2) : B(x, y +2). 
(0) B(x. y) = O +1) 


(d) B(x, y) = ol, y=1) 
Trcsc(Ty) 


(e) Para y £ Z se cumple B(x, y) - B(x + y,1— y) = pa 


Solución. Aquí usamos la relación de la función Beta con la función Gamma, es 


decir: 
PCO)PO) 


B(x, y) = CESE 
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a) tenemos que: 


ee + yYNPC6)PO) es AP6)P0)  yPEJPO) 
(x + y) (x + y) Mi+y+D T+y+D 


entonces tenemos que: 


TG+DIO). TOO)TG+10 


di O EE MA E 


= B(x+1, y) + B(x, y + 1). 


b) Tenemos que: 


TIGO) TE+NTE) _ TOO) 
PP TG+N TG+y+D TU+y+em? 


entonces multiplicando este resultado por Pr + », tenemos que 


TO0+DPOTOrEa) 
TO+DIGM+y“1+2) 
IDO) TOTO +2) 
TO+z) I(x+y+2) 
=P, 2): Ple, y +2). 


B(x, y): Pla + y,2) = 


c) Tenemos que: 


TOTO +1) 

Tx+y+1) 

__JPCIPO) 
(«+ y) + y) 


Bl, y +1) = 


e 
ETE 


por lo tanto, 


x+ 
Bas y) = Bo y +0) 
d) Tenemos que: 


_ TE+DIOG-D_AT60-5% 


XxX 
ds T(x+ y) 1) y 


por lo tanto, 


1 
BG,y)= —B(:+1,y 1) 
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e) Tenemos que: 


Bl, y) Ple + y,1— y) TOTO) TE+NTAa—-y 


T(x+ y) TG x+10 
_PE)PO)POG - y) 
AS xT(x) 


y usando la fórmula de reflexión de la función Gamma 
MYNTU- y) = rese(xy) 
para y É Z, entonces 


2. Pruebe que la función beta está dada por el producto infinito: 


da 
Bix, y) = “Ta > 


a) 


Solución. Usando el producto de Weiertrass, tenemos que: 


eYO+y) Tre e 
T(9T 0) 2 (14 DA+ 
B(x, y) = = y 
T(x + y) eYA+N poo ek 
x+y ¡paa +) 


cancelando, obtenemos que: 


= 1 
per y) = 211 HA 


3. Verifique la siguiente identidad para n,k € ZF yk<mn. 


nY 1 
k) M+DBM—=k+1k+1) 


3.2. Función Beta 141 


Solución. Usando la proposición 1.2.28 tenemos que: 


TM-k+DTK+D_TM-k+DT(k+1) 


Bin=k>+1,k+1)= = 
CO pal Tm +2) 


y dado que n,k e Z*, tenemos: 


(1. =KJU!- (1 EJ 
M+D! — (M+Da!” 


por lo tanto: 


nY 1 
k) (M+DBM-k+1,kk+1) 


4. Si-1 <q < $ pruebe que: 


+0 Gima, 1 part pray (E) (és) 
/ Coshy8 7 2 ( a ) > 


Solución. Usando el cambio de variable u = sinh?*(x) y la igualdad cosh*(x) — 
sinh?(x) = 1, tenemos que: 


py ICA = u? du 

o  (cosh(x))8 o (1+ uy? 2u2(1+u)? 
1 /(% uE du 

y ), (+) 


Ahora usando la proposición 1.2.26 y la propiedad simétrica de la función beta, 
tenemos que: 


+00 (sinh(x))* Ml 
/ ao a 


1 
mM 


a+l Pa 
ACE, 
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5. Si 0 € (—x1,70), pruebe que: 


+oo poo ab 
/ ¡) ee +y*+2xy cos(9) q x dy e Ñ 
0 0 2sin(0) 


6. Si f(x) = tan(x) o f(x) = cot(x) y a € (0, 1) verifique que se cumplen las 


siguientes igualdades: 


[acoras = Ñ a 1 UY (+) dx = E) 


raya Ja y 


7. Se llaman integrales de Wallis a los términos de la sucesión de integrales: 


T 


3 2 
W, = / sin” (x)dx = Í/ cos” (x)dx. 
0 0 
Usando estas integrales, pruebe que: 


(a) Se cumple la igualdad, es decir: 


T 


f sin” (x)dx = MA cos” (x)dx. 
0 0 


1 ¡0% Qn)! 
a ( o 5) ETS 
e _ Lo 0 IN ed 10 
Marto a 
1 n+3 1 n+l 
dd) Wer =B , = W,. 
(50 Ha 25 5) n+2 ” 
3 
Qr)? 
(e) Wy == 1112 
(r (4) 
2 
(ra) 
Pili > 
a 2421 
8. Pruebe las siguientes igualdades: 
da n=0 


1 i al 
(a) [aa Sd Qn—- DM. 
> NAT 


3.2. Función Beta 


1 ; TC; n=0 
1x2) 2x2 = 2n — DI 
(b) Fa UR A 
Qm! 
1 92n+1 n!-n! -1 
2n +1)! a 
a : e n=0,1,2,... 
Qn+D!! 
si 2m)!! 
mia RP ( 
/ A E TT 


(n — 5) im! 


1 
22(1_x?Pdx= 2. >>> 
A 


9. Evalué la siguiente integral en términos de la función beta: 


1 
/ a+ 0*41-xdx. 
-1 


Solución. Usando el cambio de variables: x = cos(0), tenemos que: 


9ga+b+1 IN cos?a+1 (5) sin22+1 (5) d0 
. 2 2 


y haciendo el cambio u = E obtenemos que: 


zx 
2 
ga+b+1 d 2/ cos24+1 (u) sin?22+1 (u)du 
0 


y usando la proposición 1.2.25, tenemos: 


1 
1 (1+10*(1=x0 dx = 2 ++ Bla +1,b +0. 
-1 


10. Pruebe, por medio de la función beta, que: 


a dx e TT 
/ (2 — x)-%(x —1)2 — sin(sro) 


para0 <q <l. 
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11. 
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Solución. Aquí usamos el siguiente cambio de variables: u = ==, conx=—t = 
u(z—t)yz=x = (z—1t)(l —u). Observe que ahora los límites de integración son 
de O hasta 1, por lo que: 


1 
eE =] A 
0 


(z — x) 1-9 (x — 1£)2 u%(1 —uyiza 


por lo tanto: 


dx AS 


(2 — xx — 12 r( 4) sin(ma)' 


Pruebe la integral de Dirichlet: 


f frrraxa E Bíp+1,q+1) 
p+qw+2 


donde el límite de integración es el triángulo formado por los ejes positivos x e y y 
la recta x + y = 1. 


Solución. Escribiendo la integral usando los límites de integración, tenemos que: 


l=x 
xP dx yldy = 2 1-05 dx 
lí / e En 


_BO + e + 2) 
q+1 
y como B(p+1,q+2)= E TOS 2 = E A entonces: 
Plo+La4+2)__ Mi+DIG+D _fir+La+D 
q+1 (p+9+2D(p+3+2) p+qa+2 ' 


12. Defina la representación integral de las funciones de Bessel por: 


Tm2(v 


Jy(z) = En E le sin?” (0) cos(z cos(9))d0 
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1 
para R(v) > => 
Pruebe que : 


(=py* 24 
Jy(z) = o ) > 


que corresponde a una representación en serie de las funciones de Bessel. 


3.3 Constante de Euler-Mascheroni 


1. Encuentre la recta tangente en x = 1 para la función Gamma. 


Solución. Sabemos que dada una función diferenciable, f(x), la recta tangente a 
su gráfico en el punto (xo, yo) está dada por y = yo + f'(xo)(x — xo). En nuestro 
caso, tenemos que: xo = 1, yo = 1 y T''(1) = —y. Por lo tanto, la recta tangente 
es: 

y=-y:x+y+1l. 


2. Verifique las siguientes igualdades: 
[,0) 
(a) / e” "In(r)dr = —y. 
0 
[,0) 
(b) Í/ re" In(r)dr =1-y. 
0 


(c) y re" n(r)dr = (n- D!+n de riole=" In(ridr. 
0 0 


Solución. Usando la función Gamma: I'(x) = des 171 dt y derivándola con 
respecto a x, tenemos que: 


T(x) = 1) PL In(Pdt. 


Evaluando esta expresión en x = 1, tenemos: 


/ e” "In(r)dt =— 
0 


ra) 
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Sea ly = qa r" e" In(r)dr, integrando por partes tenemos que: 
[0,0] 
I, ==—r"In(re "|P + / riole=" (n In(r) + Ddr 
0 
[0,0] 
=0 + nI,-1 + / rrieTdr=nlp-1+(1- Dt 
0 


que es la parte c), del ejercicio. 


Para la parte a) tenemos que /y = —y, entonces para la parte b), obtenemos: 


L=01+1f=1-y. 


Observar que 1, = T"(2). 


3. Verifique las siguientes representaciones de y: 


l 00 
(a) y = =] In=In(0)dr = + di 


x—1 
ÓS [ [ CE 


(c) y = E 1 0 ) . Donde ¿(n) es la función zeta de Riemann, definida 


n=2 
como: 


4. Pruebe los siguientes límites: 
E 1 
(a) lim ( - roo) = 
x>0NxX 
; 1 
(b) lim (- - ren) 0 
x>0 Xx 


Solución. a) Tenemos que: 


lim (- = ro) = lim (==) 
x>0x x>0 Xx 


( ra LN-E0) 
= lim | - tt 
x>0 X 
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y esta es la definición de la derivada en el punto x = 1 y como I”'(1) = —y, 
por lo tanto, 


1 
lim (2-70) === 
x>0Mx 
b) Similar a la parte a). 
3.4 Función Digamma y Polygamma 
1. Pruebe que: 
a a + W(x) 1(3S 
XEM) == += 4 + ——Á n= A 

xo x+l x+n-1 En GS 


Solución. Usando la fórmula de recurrencia de la función digamma: tenemos que: 


W(x +1) = W(x) +2 
1 


W(x +2) =Y 1 —— 
+D)= P+D) + 


1 1 
NE 
$ x+1 


1 


x+n->—1 
1 1 


+ 
x+n->—1 x+n->-2 


W(x +n) =W(x+n-1)+ 


= V(x+n-2)+ 


1 1 
A A 
6) x+1 x+n->-1 
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2. Pruebe que para p,q €Zy0< p < q se cumple: 


1 
y (2) =-— y —In(2q)- 7 cot (2) > 
q 2 q 


212 Y om(22) msn (22). 


Solución. Ver Andrews, Askey y Roy (1999), página 13. 


3. Calcule los siguientes valores: 


(a) T'0). 
(b) T'G). 
(c) P*(4). 


Solución. Usando el ejercicio 2 de la sección 3.3 (ejercicios relacionados con la 
constante de Euler- Mascheroni), tenemos que: f,, = nÍ,-1 + (n — 1)! con lp = 
TD, = P'Q),..., por lo tanto, 

a TO)=fh=l1h+01=-y+1. 

b) TB) = BL =21 +1! =-2y +3. 

e) (4) = Bb =3L+2!= —6y +11. 


4. Usando la fórmula de multiplicación de Gauss pruebe que: 


n—1 


090) 1000 + Y w (+ E). 


k=0 


Solución. Aplicando logaritmo natural a la fórmula de multiplicación de Gauss, da- 
da en el teorema 1.2.33, tenemos que: 


n—-1 

k 1 1 

y (s + >) =ó ' nQr) + 6 =nx) Inr) + hh PF (nx), 
n 

k=0 


luego derivando con respecto a x, obtenemos que: 


Abs) (mx) 
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por lo tanto, 
n—-1 


k 
090) 1000 + Y 0 (+ 7). 
k=0 
5. Pruebe que: 

1 | (2m)! yz 

r” =|=|-y-2In(2)+2 . 
(+ 3) (> bi No) 291 
param =0,1,2,.... 


Solución. Consideramos la expresión en serie para W(x) dada en la sección 1.2.7. 
Evaluando Y (3), tenemos que: 


dl 
BÉ pa 2 al) 
qye 
k=1 
=-— y —2!1n(2). 
1 1 1 
Como ro) = HH) tenemos 
PG) = (y —21n(2) YT. 


Con lo que verificamos lo solicitado para m = 0. 
Ahora, para m > 1, tenemos: 


sS) 
m+z)|=-y = 1 
2 E ml 
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Ya que: 


tenemos que: 


[ee] m [0,0] 1 
2 = = 2% k—1 y 
=m+1 k=1 s=1 


y, reemplazando en la igualdad (*), tenemos: 


3. Ejercicios 


1 1 A | 
y jaa 1 o Y oO 
(+3) pe 22 Lit 27 1 


==. MaS. 
y — 21In(2) + 2% 


como W(x) = EE tenemos, usando la proposición 1.2.8, que: 


: 1 e 1 1 
F (+ 3) = (27 -2109+23 57) 0 (n+) 


Ñ A | 1-3---(2m-1) 
= (>7-200+23 5) ( 7 vz) 


0 = 1 1-2-3-4---(2n—1)-2n 
- á 20+23 7 )( 2-4--:2n va) 


= (-»- 2m0)+23> 3 . q) or 


6. Usando la representación en serie de la función digamma, pruebe que: 


2 | 
—1In(2 —— |. 
Diao” | a+ 35) 


param e Z*, 
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7. Verifique la siguiente igualdad: 


SS Oo 1/1,f1 a 
2 ant 16 (y e (5) 


La aproximación de la suma de esta serie se conoce como la constante de Catalan 
Mes 


Solución. Para este ejercicio usaremos el caso n = 1 de la proposición 1.2.40, 
Tenemos la siguiente expresión: 


A E 1171 Ñ 
Cn+ 12 112 32" 52 pea 


n=0 
SÍ (E 16 2 16 ] 
16 NE (4n-39? E (4n—1) 
1 ($ 1 bn 1 ] 
= Ni 37 Y 3|= 
al += al en=1) 
E y! E y 3 
16 4 4 
8. Verifique la representación integral de la constante de Catalan: 


tan”! lin(1 
%.=/ af ao) 
0 [ o LFP 


9. Verifique la representación integral de la función polygamma: 


n_,—xt 


ya) = (-1 n+1 dd re 
(e) =(=1) il 


1] —e* 


dt. 


Solución. Usando el teorema 1.2.38, tenemos que: 


00 ent ext 
Yv(x)= AS 
did / ( t 1 — ) él 
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entonces derivando con respecto a x, obtenemos: 


[0,0) —xt 
w(x) = / te 
0 


1—e* 


y si derivamos nuevamente respecto a x, está última integral, nos da que: 


0/9) (2 —xt 
w"(x) = -/ e 
0 


1] — et 
Inductivamente: 
00 [Me xt 
YM (x) = (19141 / dt. 
0 l— et 


10. Verifique las siguientes igualdades para todo x, y € RAIZ” y conx X y. 


E A) y 
E (+) + y) 


Y +D)-VO+D + 


b +0), 
) (<,y)>(0,0) x— y 6 El ) 


1 
Aquic(s) = Mi —, denota la función zeta de Riemann, definida para valores 


complejos s € C con re(s) > 1. 


Solución. a) Usando la representación en serie de la función digamma, tenemos 
que: 
Wíx+1)-W(y +1) = 


E n+x=n-> y 
2 1 


SI (1+x)nm+ y) 


n=1 
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b) Usando la parte a), tenemos que: 


YVx+1)—v(iy+5n) E 
x-y =D 


n=1 


entonces cuando (x, y) > (0, 0), nos da que 


Wi tl=viy lt. 1 
ES xy =2 Da 
sea | 
2 
n=1 
=((2) 
nm? 


: 1? Es 
Nota Para la igualdad 5 =>» 1 — ver Belanzario (2001). 
n 


11. Pruebe que: 
TOTO 
TTé=1 (0 —k) 


T(x=n)= 


Solución. Usando la fórmula de recurrencia de la función Gamma, tenemos que 
P(<) = (x — 1) (x< — 1), luego derivando con respecto a x nos da: 


TM) =IG-D+(G-D5(x-1). 


De esta forma vemos que la expresión señalada vale para n = 1. Ahora de I'(x — 
1) = M(x-2+ 1) = (x — 2) (x — 2), tenemos: P'(x—1) = T(x-—2) + (x — 
2)" (x —2) y entonces 


Dx) =D) +(— 1) (Tx —2) + (x— 2) (x —2)) = 
=TMx-D+(G-DITG-2+(-— D(-2T(x-2) = 
Tx) T(x) 


=p tay te DE-I-2 
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. De aquí obtenemos: 


i 1 1 
nor rol + 5) 
(x — 1)(x — 2) 


T(x-2)= 


que es la expresión señalada para n = 2. 


Sucesivamente, tenemos que: 


n 1 n 
Tx) =1(9) Pe TM—=m] |). 
k=1 


k=1 
por tanto: 
n 1 
A ES 
[Té=1 (0% —k) 
12. Pruebe que: 
T”(x) 2 
w(x)= — (y z 
(0= 3-40) 
A cis pia Tx) : 
Solución. Usando la definición de la función digamma: Y (x) = To) y derivando 
x 
con respecto a x, tenemos que: 
TOTO -UTE"0y TI" 
Y (x) = 0) —( ¡A 0 _ ue. 


My) MES 


13. Aplicando el ejercicio anterior verifique las siguientes igualdades: 
poa 
(a) 10) = al v=¿0)+y. 
1 
1) PD) =-y? — my 203). 
3 
(0) TUD =y* + 1%y? + 8€(3)y + —r*. 


20 


14. Demuestre que para a,b e R* y a 4 b, se cumplen las siguientes igualdades: 


a STO ¡QS 
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(0) e O dr = 4 ($) 
% cos(bx) —cos(ax) a 
of OO o tm (2). 


15. Verifique las siguientes identidades simétricas: 


() LB) = BG, y) (2) — + y) 
9 
(0) ¿7800 y)= By O) — A y) 


Solución. Para ambos ejercicios usamos la relación B(x, y) = ELOpO) 


Tx+y) * 
a) 
Py O) 
ly) 
Z p QIOTA+YN TOO y) 
(MG + y)? (PG + y)? 


POP) _ TOTO) + y) 
T(x + y) T(x + y (x + y) 

PITO) 

DNCTRCES B(x, y)V(x + y) 

= B(x, y) (x) — B(x, y) (x + y) 

= B(x, y W(x) — W(x + y). 


b) Ídem ítem a), pero derivando con respecto a y. 


3.5 Función Gamma incompleta 


1. Pruebe que: 


xy la=0D! a 
E 


Solución. Integrando por partes la integral definida para y(a, x), tenemos: 


a,¿—X 1 Xx 
y(a,x) = FP tay o E =/ tte “dt 
0 a a Jo 
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156 
e integrando por partes vamos a tener la serie: 
xXx xa+lox 
a,x)= + =p 
y (05%) a ala +1) 
por lo tanto, 
A (a — Da 
Xx 
os » ano 
2. Pruebe que: 
d” 
() A] = (PAT (a + m,x), 


d” *P 
0) qq leva = Y 


3. Pruebe que y(a, x) y I'(a, x) satisfacen las siguientes relaciones de recurrencia 


(a — m, x). 


(a) y(a+1,x) = ay(a,x)-x*e*, 
(b) Pla +1,x) =aT(a,x) + x*e *, 


Solución. (a) Podemos escribir: 


Ñ x d(t%e7t) 
A 1 dt ñ 


o 
=/ (ar 1g7t — xt) dt = 
0 
eS Xx 
Sl area | x“e tdt = 
0 0 


=ay (a, x) — y(a + 1,x), 


por lo tanto, 
a XxX 


y(a + 1,x) = ay(a,x)-x*e 


(b) Ídem ítem a), pero usando la definición de (a, x). 


4. Se define como la integral exponencial a la siguiente expresión 


Ej(=x) = Er (x) = y ar, 
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De igual forma se definen la integral coseno, la integral seno y la integral logaritmo, 
como: 


* Sib) = -f" may 


: ci =- [Da 


t 


' * dt 
> Liqx)= / (0) 


respectivamente. Pruebe que: 


00 pt oo (=D"x 


E 1”x 2n+1 
2n + DOn+ DU 


(b) Si(x) = 5 + 2 


E D”x 2n 
(c) Ci(x) = y +00 +2 Sr ant 


5. Si relacionamos la integral exponencial con un cambio de variables podemos obte- 


ner que: 
[0,0) ext 
En(x) =/ A dt. 
Xx 


Pruebe que £,, (x) satisface la siguiente relación de recurrencia 


xk 

E = e *- —_=— = Lidia 
n+1(x) e ECO para n 

Además pruebe que E, (0) = 27, paran > 1. 


6. Se definen las integrales de error a las siguientes expresiones: 


2 E =P 
» erfíz)= F) e 
* erfelz)= Af 4 


Pruebe que: 
(a) erfc(z) = 1— erf(z). 
(b) erf(z) =117*y (5) 
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¿il 
(c) erfcl(z) = 7 ir (52). 


2 a 2 Re 
Solución. (a) Tenemos: erf(z)+erfc(z) = Fl) Par] Pd 
T Jo TJz 


2 [ee] 
= — (/ Pdo): usando la proposición 1.2.6, vemos que: 
YT Jo 
2 O 2. yx 
erfclz) + erfíz) = + e“nd|==+—-L=1. 
rro) 


Xx 
(b) Sabemos que y (a, x) = / t91671 dt, entonces con a = > Y = 2% 08 


2 
1 z 
da y (5 2) = / 173 e7* dt. Usando ahora el cambio de variable £ = u?, 
0 


tenemos: 


2 


/ ietdt= 2/ du = Vx -erfle). 
o 0 


(c) Ídem ítem b). 


3.6 Función Gamma en los números complejos 


l. Paraze C y Relz) 4 0,—1,-2,... tenemos que: 
T()=T(2) 


nTE)=a 760). 


Solución. Usando el producto de Weierstrass como representación de la función 
Gamma, la proposición 1.2.20, y las propiedades de los conjugados de los números 
complejos, tenemos que: 


porlo tanto, P"(z) = T(Z). De este resultado se deduce inmediatamente que: In(7"(Z)) = 
In(TG)). 
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2. Seax € R,conx É£ Z., tenemos que: 


is 
(691 sinh(7]x)' 


Solución. De la definición de la función Gamma para una variable compleja, pode- 
mos, al igual que en el caso real, probar que vale la fórmula de reflexión de Euler y, 
por ello, tenemos que: 


M2DPO-—z)=-— , 
sin(7z) 
Así, para z = —¡x, con x X O, tenemos: 
. T 
sin—i1x) 
y dado que sin(x) = mit entonces 
PEIOTa+ ix) = 2 E ds 
—ix ix) = == € EOS 
1er) rc 
Dado que sinh(x) = 2, obtenemos: 
PEIÓTa+i)= == 6.1) 
—ix x= ——. . 
22 Gmh(rx) 


Ahora de la fórmula de recurrencia de la función Gamma, '(z + 1) = zI (2), y de 
la definición factorial de la función Gamma, PF (z + 1) = z!, deducimos que 


Ed)! = Peix+DP(Gx+1)=-ixI ix (ix +1) 
y usando la ecuación (3.1), tenemos: 


ir xT 


APIS an 


(8.2) 


Ahora, usando el producto de Weierstrass para la función Gamma, proposición 1.2.20, 
tenemos que: 


160 


entonces siz =u— iv con u,v € R nos da: 
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cid u—iv+l 
E dd sai Y y PO L a 
lu—iv+1| ¿2111 + A) 
grhertl My ne 
EME a VU +u + 1)? + v2 
y que: 
Y EN 


[(u + i0)!] =|(u+iv+ 1) = 


: HE 
upio+1 Ll a E] 
eru+1) E se 
VU+02 +02, y(M+u + 1)? +2 
con lo que |(u — ¿v)!] = |(u + ¿v)!|, por lo tanto, usando la ecuación (3.2), u = 0 
y UV = —x, tenemos que: 
xx xx 
01)! = |———— |, entonces |(¡x)11? = —“—. 
a sinh(7rx) Ad sinh(7x) 


como señalado. 


. Para n y b números enteros positivos, se cumple: 


Tb 


(1 +1ib)!] = (25) Y AREA 


Solución. De la fórmula de recurrencia y la definición factorial de la función Gam- 


ma tenemos que: 


+ibdl=TPM+ib+D=(M+ibM+ib-DT(m+ib=1) 


y dado que n es un entero positivo, podemos escribirla de la siguiente forma: 


M+ib!l=FM+ib+1)= (+ ibNa+ib=1)--. 


y como "(1 + 1b)) = (ib)! tenemos que: 


Pm+ib+D|=|0+ ib) +ib=D|-- 


(A+ ibFPdw+ib) 


[4 + 1b)11(b)!.. 
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1 


o ; . xT 2 
Del ejercicio anterior sabemos que |(¿x)!] = | —¡— ] , por lo tanto: 
sinh(7x) 


1 


br 2 
A CAE IA A MA 
E CAN (5) 


con lo que: 


Tb 
lo io m1 = (3 a llas p2)3, 


como señalado. 


4. Se cumple 


29 p?2 -3 
resi ira] 


2 
pe (a +1) 


para cada número complejo que no es un entero no positivo. 


Solución. De la representación de la función Gamma como producto infinito, tene- 
mos que: 


Ahora, para z = a: +18, tenemos: z +7 = 2a y z-7 = q? + f?, por lo que 
obtenemos: 


1 AS =z a zz =Z 


n=1 n=1 


[0,0] 

pr z 
=z-zeY6+2) ne =) 

Ide 


- (0? + B?2Je2v2 ml ( + 8) 
n=1 


Considerando $ = 0 en (3.3), tenemos que: 


1 [oe] 


o ay -2a 
np AS 
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entonces 


20 1 
er |] = (3.4) 
n=1 al, (1 + =y [PT (a)? 


reemplazando (3.4) en (3.3), nos da que: 


1 o arpa a (1+ 58) (145) 
TOTO [Mar (+2) 


n=1 
a idad y [add 
NC 


Ñ 1 a? + 2 oo p? 
Or. [1 (1 ar): 


n=1 


por lo tanto, 
1 


pe 
TOTO TO E Ú (1 (0 any) 


y como sabemos que |" (a + ¡B)| = | (a: — ¿B)|, tenemos 


[99] 2 -1 
¡Tía +iB)1? = rr TT + em) y, entonces 
Ñ eN p? <a 
resi ira E) 


como señalado. 


5. Paraa: 4 0,—1,-2,...y 6 € R se cumple que: 


[Pla +i8)| < |P(a)l. 


Solución. Del ejercicio anterior tenemos que: 


eS p2 =3 
Tr +i81=|01 || ( y E) 


Bebe (a + ny)? 
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y como: 


5 el 
ral) <P), 


n=0 


por lo tanto, 
[Pa +i8)| < |P(0)l. 


6. Para y € R, se cumple: 


2 


— cosh(ry) 


T 


1 
Pe 
rl+) 


Solución. La fórmula de reflexión de Euler nos señala que 


MDPO 2) = En Usando, en esta expresión, el cambio de variable z = 
Ly; , 
7 + 1y, tenemos: 
1 1 TT 
Pl=+iyPGiyY= 
(5 ») 5 ») sin (m2 + iy)) 
T 


sin(5 +11 y) 
T 


Ñ sin(<) cos(iyz) + cos(5) sin(iy7) 


T 


costiya) sd 


e ya) y griliya) eY7 y eya 


y como cos(iyx) = 3 = 5) = cosh( y), por lo que 
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reemplazando en (3.5), nos da que: 


r 1 de ra iy) TT 
=+1i == 1) = ———, 
2 dl 2 Y cosh(y7) 

por lo tanto, 

r 1 + a TT 
=+i = ———. 
2 dl cosh(71 y) 


7. Para y € R, se cumple: 


0 O O NO 1/2 
p 1 = —l = . 
4 dl 4 iS cosh(x y) + i sinh(7 y) 


Solución. Sea z = - + iy, entonces f” ( + iy)T (3 iy) =PG+DIE2). 
Ahora sea w = z + l, entonces '(z + 1) (=z) = [(w)I'(1 — w) y usando la 
fórmula de reflexión para la función Gamma tenemos que: 


Fwd -—uw) MENE 


T 
sin(r(3 + iy)) 
_ TT 
sin(% +11 y) 
T 
= sin(%) cos(ixy) + cos(%) sin(i7ry) 
T 


AP cos(irry) + sin(imy)). 
y como cos(ixy) = cosh(x y) y sin(imy) = i sinh(i y), entonces: 
27 
V2(cosh(x y) + ¡ sinh(x1 y)) 
2271 
—2(cosh(xr y) + i sinh(71 y) 


Ñ ay/2 
—cosh(1 y) + isinh(7r y)” 


Mw -—w) = 
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por lo tanto: 


( ) E ) Ty? 
Pl-=+iy|P | -—iy)]= a , 
4 4 cosh(x y) + i sinh(7 y) 


8. Calcule el valor numérico (aproximado) de: 


(a) T(. 

(b) FA+I. 
(e) PU —i). 
(d) T(S+3i). 
(e) P6-—3i). 


1. 1. 
1 


Solución. (a) Usar la aproximación de la función Gamma, proposición 1.2.32, y 
evaluar en /" (1), esto es: 


1 139 571 
r VA (14 
io A. y 7 288x2 — 31840x3 ell 


e 1 1 139 57 
TG) = V2ri2e (1 
(1) ile t(1+ 10i * 288i2 — 51840i3 — 24883201 


a 


ñ ALA 0d : Tí Pe 1 
Acá, podemos « escribir la expresión ¿ = e*?, entonces ¡' 3= (ei Fyi=2 = 
— HE_ ¡Tí == j de 
ectit=e? (cos $ Z —i sin %) y la expresión e? = cos(1)—i sin(1), por 
lo tanto, la expresión: 


V2miiE pi =V21 e? (eos E z — i sin 7): (cos(1) — i sin(1)) 
= — 0, 110968 — 0, 5091245. 


Por otro la dos la expresión: 


a 139 571 
¡ "28812 5184013 248832011 
Le Pi 


12' 7 288 * 31840. 2488320 
= 0, 996298 — 0, 0806521, 
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10. 
11. 


12. 
13. 
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entonces 
T (1) = (0, 110968 — 0, 509124¡) - (0, 996298 — 0, 0806521), 
por lo tanto, 
T (1) = —0, 1549 — 0, 49801. 
(b) Aplicar los mismos pasos del ítem a), ítem b) e ítem c). 
(c) Aplicar los mismos pasos del ítem a), ítem b) e ítem c). 
(d) Aplicar los mismos pasos del ítem a), ítem b) e item c). 


(e) Aplicar los mismos pasos del ítem a), ítem b) e ítem c). 


. Pruebe que: 


arg (P(z)) = yW (x) + (=, — aran 2 )) ] 
2 x+n x+mn 


conz=x+iyyzHXH0. 

Pruebe que "(1 + iy) = iyI (iy). 

Pruebe que: E 
|TGy) 1? = y sinh(zr y)" 

Pruebe que si F(x + iy) = u + iv, entonces P(x — iy) =u-—iv. 


Usando las definiciones de las integrales exponencial, coseno y seno, pruebe que: 


Es(ix) = Ci(x) + ¿Si(x). 


Puntos 
maximos y/o 
minimos de la 
función gamma 


En este anexo se muestran métodos para encontrar los puntos de máximos y/o mínimos 
de la función Gamma. 


A.1 Método I 


De la función digamma tenemos que 

PU 1 | 
De 
T(x) x n(n + x) 


n=1 


entonces para encontrar sus valores extremos debemos resolver la ecuación: 


1 [0,8] 
d A 
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claro que para encontrar las raíces de esta ecuación es un problema complejo. Por ello 


consideraremos una aproximación 
M 
1 1 
A, => 
n=1 
que es lo mismo que 
| 
1+yx—x? —_—— = 
1 


La apropiada elección, del índice extremo M, va a depender de la rápida convergencia de 


[0,0] 
1 
la serie —————. Podemos escribir 
2, n(n+x) 
00 Pp 00 1 
== 9 o” al n(n+x) + ai n(n + x) 
a pa) 
Sp 
[0,0) 
1 -1 
=Sp+ YD) 2 (1+27) 
n=p+1 


z 1 2 1 1 
n n+x) 


(AER 
rd ze n(nN+x)  x 
n=1 n=1 
Ahora, para cierto intervalo (a,b) € D (D dominio de la función Gamma), le impon- 


dremos las siguientes condiciones 


l. b-a<l. 


2. x € [a,b] C (a,b). 
3. p= [|x|] + 1, donde [|x|] indica la parte entera de |x| 


Entonces si reescribimos 
[0,0) 
1 Xx 
(143), 


Il 

Il 
eb 
ES 
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A.1. Método I 
entonces 
Y) 7 = Sp (0) + > AS 1) Ñ 
n=1 n(n +x) n= dl E a 
Es 1 
= Sp(x) + Nene Y =P 00+ 
k=0 n=p+1 
dd DL a ya a 
n=p+1 n=p+1 n=p+1 
| 
kk 
SA Y pk+2 
n=p+1 
a el ce E | 
e OI AD il ID 
n=1 n=1 n=1 n=1 
3 o 1 SL m7 Ent > 1 e 
+x 2472 L ni 
| 
e | 1 
luego usamos la función zeta de Riemamn, es decir, ¿p(k) = YE) 0 Y 
n=1 n=1 
tenemos: 
DEXE (Ek +2 k+2 
a” Sp) + DOE (2d +2) Ep de +2) 
n=1 k=0 
o [0,0) [0,0) 
28 AA = E Ó) 
ra = SEDA (Edo — Ep (4) 
n=1 k=2 
entonces 


[0,0] 1 00 
L+yx A AS e E AA) (A.1) 
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luego, eligiendo un M adecuado, tenemos: 


M 
14 yx —x28p (2) +) EDEx (E(k) — Ep(k)) =0. (4.2) 


k=2 


Debido a que éste algoritmo de optimización depende de la función zeta de Riemann y 
de la rápida convergencia de ésta, es decir, mientras sea mas grande M más rápida la 
convergencia de (A.1). Esto implica claramente la minimización de p v/s x, por ello lo 
importante del valor de p, p = [|x|] + 1. 


Con estos cálculos tenemos, ahora, que encontrar las raíces de la ecuación (A.2) va a 
depender del cálculo de los siguientes parámetros: 
+ El número y yy de cifras significativas (precisión). 


» El valor de M(?) < nm. 


+ El número de cifra significativas de 7 y b. 


Ejemplo A.1.1. A modo de ejemplo, consideremos el caso en que x € (a,b) = (-4,-—3), 
entonces según el método tenemos que p = 4. La ecuación A.2 nos da: 


M 
0 =1+ yx — x?Sa(x) — Y Ent) — fa(k)) 
n=1 
2 a 5 1 
=1+x 193 (5- =) Y Ent (sco- 2 5) 
k=2 


n=1 n=1 


Con pocas iteraciones sucesivas, podemos afirmar que el rango de [a, b] = [-3, 6353; —3, 6352]. 
Ahora si asignamos M = ny = 10(?), tenemos que x == —3, 6352933664, 


A.2 Método II 


Otro método para encontrar los máximos y mínimos de la función Gamma para x < 0, 
es usar la fórmula de recurrencia y reflexión de la función digamma y el método de Newton- 
Raphson. 


Para encontrar los máximos y mínimos de cualquier función, primero se calcula la 


derivada de esta y luego se iguala a O, en nuestro caso I"'(x) = 0, pero según la defini- 
/ 


ción de la función digamma W(x) = 


W(x) =0. 


, entonces tenemos que resolver la ecuación 


(x 
MES) 
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Primero usaremos la fórmula de reflexión de la función digamma: 
W(l—x)-W(x) = xmcot(x); Vx E Z. 


Sea x < 0 tal que W(x) = 0, entonces W(1 — x) = x cot(rx), y como x < 0, entonces 
tenemos que 1 — x > 0 y si hacemos el cambio de variables 1 — x = 1 + y, entonces 
=x = y con y > 0 por tanto, tenemos que: 


W(y + 1) = x cot(7 - —y) = —x cot(7 y), 


entonces 
W(y + 1) + reot(xmy) =0 (A.3) 
y como para y > 0, tenemos la aproximación de la función digamma: 
1 ye lo 
12y2 — 240y*” 


1 
YW(y) =1 - = 
(y) = In(y) 3y 
por tanto, usando la fórmula de recurrencia de la función digamma: 
1 
YU DO=P0)+ 7 


tenemos: 
1 de 1 mn 1 
12y2 * 240y* ' y” 


1 
W(y + 1) = In(y) — AN 
y 


esto es: í 
12y2 + 240y*' 


Reemplazando en (A.3), observamos que tenemos que resolver la ecuación: 


1 
YW(y +1) =Iny)+ 5 - 
2y 


1 


1 
12y? + 240y + a cot(my) = 0; 


1 
1 se 
n(y) + 2y 


y para hacerlo usamos el método de Newton-Raphson con : 


1 1 1 
h(y) =1 =- t ; 
O) = In) + 271 + 240y1 + 71 cot(ry) 
para la cual tenemos : 
1 1 1 1 

hy) === + - — 2? esc? (71 y); 

(y) E (xy) 

1) 


como el método de Newton-Raphson, señala que Y»+1 = Yn — , y como sabemos 


h(Yn) 
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que hay un punto de máximo o de mínimo en cada intervalo x e] —n,—n + 1[, entonces 
y €ln — 1,n[ paran e N. 

Para el intervalo ]O, 1[, empezamos con el yo = 0, 5 y nos da que y = 0, 5044319425 
y como —x = y, entonces x = —0, 5044319425. Para el intervalo ]1, 2[, empezamos 
con el yo = 1,5 y nos da que y = 1,573449476 entonces x = —1,573449476. Para 
el intervalo ]2, 3[, empezamos con yy = 2,5 y nos da que y = 2,610713614 entonces 
x = —2,610713614. Y así sucesivamente. 


A.3 Método II 


Otro método, parecido al método II, se obtiene de usar la fórmula de reflexión de la 
función digamma, su aproximación y una expresión, ver Magno (s.f.), por medio de la 
cual Hermite notó de como se comportaban las raíces de Y, 


Xx =-k+0 (2) E (A.4) 


La fórmula de reflexión de la función digamma establece que: 


W(l—x)-—W(x) = xm cot(1x) = EZ 


tan(xx) 


de esta expresión tenemos que: 


Y(l—x)=VW(x)+ 


tan(1x)' 
entonces para x < O tal que W(x) = 0, tenemos: W(l — x) = y usando el cambio 
an(7rx 
de variable 1 — x = 1 + y, tenemos que 
T 
y 1) = =- E 
ad) tan(—Ty) tan( y) 
y, de esto, tenemos: 
tan(zy) z 
an(my) = == 
2= "YO ED 


T 
Ty = arctan | — 
j ( wy + a) 


1 
y= me arctan ( (A.S) 


7070) 
W(y +1) 


A.3. Método III 173 


y usando (A.4) en que xz = —k, entonces como —x = y, tenemos que yz = k, entonces 
Yx — k =0 y usando la aproximación de la función digamma, obtenemos: 


1 1 
12y2 dl 240y*” 


1 
Vly +1) n(y)+= 
2y 


entonces reemplazando en (A.5), tenemos que: 


1 , TT 
Yk—k == ——arctan E 
TT 1 1 1 
In(k - 
OF ma * aaa 
por lo tanto, 
1 TT 
Yk Y k— — arctan 1 1 1 ; 
In(k =- 
MO ma Y 2a0rs 


y con esto tenemos que: 


T 


1 
Xxx =—k + —arctan 
TT 1 dE 1 
12k2 ' 240k* 


1 
In(k 
n( NE 
conk e N. 


De estos tres métodos podemos obtener los puntos de máximo ( o de mínimo) de la 
función Gamma, para cada intervalo ] — k,—k + 1[ parak > 1. 
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Cuadro A.1: Valores de los puntos mínimos y máximos de la función Gamma para x < 0. 


A. Puntos máximos y/o minimos de la función gamma 


X 


T(x) 


—0, 5044319425 


—3, 5446455403 


—1,5734494759 


2,3024072858 


-2,6107136126 


—0, 8881363587 


3, 6352914559 


0, 2451275398 


—4, 6532370735 


—0, 0527796396 


—5, 6671621403 


0, 0093245945 


—6, 6784180625 


—0, 0013973966 


—7, 6877882422 


0,0001818784 


—8, 6957641148 


—0, 0000209253 


—9.7026725251 


0, 0000021574 


Numeros de 
Bernoulli 


Una manera de encontrar propiedades de la funciones Gamma, digamma, polygamma; 
es usando los números de Bernoulli, Oyrazo y Mancilla (2006). 


Definición B.0.1. Se define la sucesión de números (B, )nen, de manera inductiva como 
sigue: By = 1 y 


para cada m € N. 


Nota. Los números B,, reciben el nombre de números de Bernoulli, en virtud de que fue 
Jacob Bernoulli quien los definió y utilizó por primera vez. 


m+1 


Si la definición de Ba, la multiplicamos por (m + 1) = ( 


mos: 


) en ambos lados, tene- 
m-1 m-+1 
(m + 1)Bm ==) , Ja. 


1 Lia 1 
m k 
k=0 
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por lo tanto, 
2 [(m+1 
El k Jn para m > 1. 
k=0 
Si expresamos esta última ecuación para cada m € N, el resultado es el siguiente 
sistema de ecuaciones: 
2B¡+1=0, 
3B2+3B¡+1=0, 
4B3+6B2+4B¡+1=0, 
5B4+10B3+10B2+5B¡+1=0, 


y realizando los cálculos para los primeros 12 números de Bernoulli, tenemos: 


1 1 1 1 

1 025 3 =0;B4 30 5 =0;B6 7) 
A A A 
E »D8g = 30 9= > 10 = 66' 11-= > 12-= 2730 


Í 
Lema B.1. Supongamos que se expande la función 7 en serie de potencias alrededor 


t 
e — 
del origen, es decir, que se busca una expresión de la forma: 


: 00 
et—1 = de Bmt”. 
m=0 


B 
Entonces se cumple que by = E Vm € No. 
m! 
Í sas 
Prueba. Sea AA = O bt"; multiplicando con(e' —1), en ambos lados, tenemos: 
e pa 


1 = (e -1) Y bt”. 
k=0 


n 
y como e* — 1 = ) ¿21 57, entonces tenemos: 
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(Em) (2) 


Usando el producto de Cauchy: $20 4n* izo Dn = Y jp Cn CON Cn = Dogo AkDn—k- 


tenemos: 
o) m br E 
7 (E cm) 


m=0 


e ; Bo. 
Si igualamos los coeficientes en ambos lados tenemos que bo = 1 = or” mientras que 
para m € N nos da: 


E $ ko > kbz 
E m1 101 km 1 Ko! 


y multiplicando por (mm + 1)!, nos da 
se avi 2 [(m+1 
a 1h 


para todo m € N. Comparando esto último con la ecuación 


". [m>+1 
B, = 
k=0 


más la condición inicial 0!by = By = 1, tenemos que: 


AS 


por lo tanto, 
B 
= k!bx y de aquí bx = TE 


Proposición B.0.1. B>2;+1 =0, Vk € N. 


Prueba. Aplicando el lema anterior y tomando en cuenta que B¡ = 2 tenemos que: 
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entonces ES 
Í Í Br k 
A A 
k=2 
tt .e44+1 8 
a E" =1+ PE ¿ke 
2 et-1 k! 
k=2 


t 


t e+1 e 
Veamos que A es una función par: 
e — 


=t e +1 el —te41l tel+1 
2 et-1 el 2 1-et 2 et-1” 


así llegamos a la siguiente relación: 


[oe] 


Ber te+1l — e +1 En 
=2 


No 


Bk 


Si se igualan los coeficientes de 1* en la última expresión, se tiene que 1 = 1 = 
gu p q Y El 


_pe 
€ q EE , para k > 2, es decir: 


Br=(-D*B¿ Vk>2. 
Además si k es un número par, esta última relación no nos entrega ninguna información, 


pero si k es un número impar tenemos que B = —By, por lo tanto B; =0sik>2yk 
impar. 


B.1 Polinomios de Bernoulli 


Definición B.1.1. Para cada m € No, definimos el m-ésimo polinomio de Bernoulli, 
denotado por B,, (x), como: 


Bm) = Y (1) Bpxrok. 


k=0 
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De esta forma tenemos que: 


Bo(x) =1, 
1 
Bi(x) =x-— > 
a 1 
Ba[x) =x"—x+ 6 
3 1 
B3(x) =x? — qe + 3% 
1 
Bab) =* =D EAS 


30” 


1 
Proposición B.1.1. 7) m1 Bm) = Bm(x), para Vm € No con By, , ¡ (x) deriva- 


da del polinomio Bn+1(x). 


2) Bm(0) = Bm(1) = Bm, Vm € No; m H 1. 


Prueba. 1) Primero observemos que (”?*) milk = (1), entonces de la definición de 


polinomios de Bernoulli, tenemos que: 


1 / 1 d ES m-+1 m+1-k 5 
to (E O) = 


k=0 


1 dd m>+l —k 
= —o + 1—k)B¿x""* = 
m>+l yl k Jm ) Bio 


0 
A 
m+l 


2) Sim = 0, entonces el polinomio By(x) es el polinomio constante con valor By = 1. 
Ahora bien cuando m > 2, entonces es claro que Bn (0) = By. Por otra parte, se tiene 
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que: 


m-—2 
m 
=Bmy + MBmp1 + (+)a- 


k 


=0 
¡jaca Ed 
=Bm +mBm-1 —m (- E bi 
m 
k=0 
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=B im + MBm-1 = MBm-1 = Bm. 


En el caso del polinomio B¡(x), se tiene que B¡(0) = 


2 


= Bj, mientras que 


B10) = > = —Bj. Con esto se completa la información acerca de los B,(0) = 


Bm (1). 
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